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MONTHLY NOTICES 
OF THE 
ROYAL ASTRONOMICAL SOCIETY 


Vol. 117 No. 6 


MEETING OF 1957 NOVEMBER 8 
Dr W. H. Steavenson, President, in the Chair 


The election by the Council of the following Fellows was duly confirmed :— 

Matthew Leslie Billis, 5 Easton Road, Bebington, Cheshire (proposed by 
W. F. Bushell); 

Paul Laurence Byard, Beach House Cottage, Golf Place, St. Andrews, 
Scotland (proposed by A. H. Jarrett); 

Leonard Melbourne Dougherty, 63 Hildry Avenue, Liverpool, 14 (proposed 
by C. W. Allen); 

William Joseph Farquharson, 26 ‘Thornhill Road, Mount Merrion, Co. 
Dublin (proposed by H. A. Brick); 

Arthur A. Hoag, U.S. Naval Observatory, Flagstaff, Arizona, U.S.A. 
(proposed by C. H. Smiley); 

Stephen Jackson, Upway, Rode, Bath, Somerset (proposed by H. Thomson); 

Ram Singh Kushwaha, Dominion Astrophysical Observatory, Victoria, 
Canada (proposed by A. McKellar); 

Francis Joseph ‘Terence Maloney, 23 Lawn Crescent, Kew Gardens, Surrey 
(proposed by J. Vetterlein); 

Donald S. Reed, 57 Meredith Drive, Cranston, R.1., U.S.A. (proposed by 
C. H. Smiley); 

Norman Jack Rumsey, Dominion Physical Laboratory, Lower Hutt, New 
Zealand (proposed by J. C. Begg); 

Alexander Thom, Brasenose College, Oxford (proposed by F. J. M. Stratton); 

Stanley Woodall, 11 Martin Road, Springfield, Port Elizabeth, S. Africa 
(proposed by H. Welsh); and 

Harold Surridge Wren, Bradwell, Longfield Drive, Amersham, Bucks. 
(proposed by A. W. Vince). 


The election by the Council of the following Junior Members was duly 
confirmed :— 

Dorothy Joan Crampin, Royal Holloway College, Englefield Green, Surrey 
(proposed by W. H. McCrea); 

‘Terence James Deeming, Physics Department, ‘The University, Birmingham 
(proposed by R. v. d. R. Woolley); 

Adam Ford, The Rectory, Wootton Courtenay, Minehead, Somerset 
(proposed by W. L. Rae); 

Philip Walford Hill, 14 Comberton Road, Kidderminster, Worcs. (proposed 
by R. v. d. R. Woolley); 
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Stuart Robert Charles Malin, 6 Westover Road, High Wycombe, Bucks. 
(proposed by R. v. d. R. Woolley); 

Derek McNally, Royal Holloway College, Englefield Green, Surrey (proposed 
by W. H. McCrea); and 

Geoffrey Austin Mills, 1 Gordon Avenue, Ross-on-Sea, N. Wales (proposed 
by G. E. Taylor). 


Ninety presents were annpunced as having been received since the last 
meeting, including :— 
A. Armitage, Copernicus (presented by the author); 
M. W. Ovenden, Looking at the Stars (presented by Phoenix House); and 
University of Manchester, Astronomical Contributions, Series III, Jodrell 
Bank Reprints, Vol. 3 (presented by the University of Manchester 
Department of Physics). 


MEETING OF 1957 DECEMBER 13 
Dr W. H. Steavenson, President, in the Chair 


The election by the Council of the following Fellows was duly confirmed :— 

John Azgyrakos, National Polytechnic School, Athens, Greece (proposed by 
S. Plakidis) ; 

George Patrick John Barry, 14 Harell Avenue, Colchester, Essex (proposed 
by S. Boulter); 

Kenneth Beames, Linden Observatory, Sydney, Australia (proposed by 
A. York); 

Leonard Robert Hall Beaumont, P.O. Box 3127, Auckland, New Zealand 
(proposed by W. H. Ward); 

Clement S. Brainin, 320 Washington Street, Mount Vernon, N.Y., U.S.A. 
(proposed by F. Benario); 

Peter Chadwick, 11 Alexander Road, Thatcham, Berks. (proposed by 
D. W. N. Stibbs); 

Som Datt Chopra, Punjab University, India (proposed by E. R. Lapwood); 

David Wright Collinson, Physics Department, King’s College, Newcastle 
upon Tyne (proposed by S. K. Runcorn); 

George Contopoulos, University of Thessaloniki, Greece (proposed by 
S. Plakidis) ; 

Paul Fenimore Cooper, 1564 Massachusetts Avenue, Cambridge, Mass., 
U.S.A. (proposed by D. H. Menzel); 

Corban Assid Corban, 62 Remura Road, Auckland, New Zealand (proposed 
by R. A. McIntosh); 

Kenneth Midworth Creer, Physics Department, King’s College, Newcastle 
upon Tyne (proposed by S. K. Runcorn); 

Frederick Vernen Noél Croyle, Stormont House, Mutley Road, Plymouth 
(proposed by A. M. Newbegin); 

David Cumberworth, 4 Bell Street, Workington, Cumberland (proposed by 
T. L. MacDonald); 

John Davis, Jodrell Bank Experimental Station, Lower Withington, Cheshire 
(proposed by R. L. Waterfield); 
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Vishvas Shanker Dhumal, Indian Naval Dockyard, Bombay, India (proposed 
by S. A. Moos); 

Edward Gerard Docherty, 66 Raglan Street, St. Helens, Lancs. (proposed by 
R. M. Baum); 

Wladyslaw Dziewulski, University of Torun, Poland (proposed by J. 
Witkowski); 

Francis Henderson Elliot, Clay Lane, Brundall, Norwich (proposed by 
G. F. West); 

Robert Seeger Evans, 3960 Cedar Hill Road, Victoria, B.C., Canada (proposed 
by G. A. Hole); 

Peter H. Filer, 4 Sutherland Drive, Giffnock, Glasgow (proposed by 
T. R. Tannahill); 

William Arnold Granger, Rose Lea, 58 Priory Road, Peterborough, Northants. 
(proposed by P. Moore); 

Troilus Stanley Gray, 87 Doncaster Road, North Balwyn, Victoria, Australia 
(proposed by W. E. L. Clapham); 

George Francis Hamm, Private Bag, Port Arthur, Tasmania (proposed by 
T. N. Burke—Gaffney); 

David Andrews Hardy, 44 Griffins Brook Lane, Bourneville, Birmingham 
(proposed by M. Guest); 

Mervyn Frank Harrold, The Cottage, Cherry Burton, Beverley, Yorkshire 
(proposed by F. J. Acfield); 

John Hazlehurst, Dept. of Astronomy, Manchester University (proposed by 
Z. Kopal); 

Antony Hewish, Cavendish Laboratory, Cambridge (proposed by D. W. 
Dewhirst) ; 

John William James Hosken, 13 Chase Court Gardens, Enfield, Middlesex 
(proposed by D. G. Ewart); 

Frank Wilsenham Hyde, 27 Carlton Road, Clacton-on-Sea, Essex (proposed 
by J. Heywood); 

Cecil Ismay, 95 Wastwater Avenue, Westfield, Workington, Cumberland 
(proposed by T. L. MacDonald); 

Wilhelmina Iwanowska, University of Torun, Poland (proposed by J. 
Witkowski); 

John Trevor Jefferies, Sacremento Peak Observatory, Sunspot, New Mexico, 
U.S.A. (proposed by R. Giovanelli); 

John David Jones, 57 ‘Talbot Road, Port Talbot, Glam. (proposed by 
J. H. J. Westlake); 

Thomas Reeve Kaiser, Physics Department, University, Sheffield (proposed 
by C. W. Allen); 

Stephen Keeler, 71 St. Johns Grove, Bilton, Harrogate, Yorks. (proposed by 
G. Fielder); 

Branislav Kostic, 19 Myimbar Way, Nollamara, W. Australia (proposed by 
R. W. Boggis); 

Demetrius Kotsakis, University, Athens, Greece (proposed by S. Plakidis); 

Elizabeth Phillips Harding Lees, Narmbool, Barwon Heads, Victoria, 
Australia (proposed by H. Dawson); 

Martin Rastall Coulter McDowell, Royal Holloway College, Englefield 
Green, Surrey (proposed by W. H. McCrea); 
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Walter E. Mitchell, University Observatory, Ann Arbor, Michigan, U.S.A. 
(proposed by C. H. Smiley); 

Peter John Noonan, 32 Parklands Road, Streatham, London, S.W.16 
(proposed by H. J. Lewis); 

Geoffrey Page, 30 James Street, Tong Road, Armley, Leeds (proposed by 
N. H. Silverton); 

Henry P. Palmer, Jodrell Bank Experimental Station, Lower Withington, 
Cheshire (proposed by A. C. B. Lovell); 

William Robert Pearson, 16 Mellalieu Street, Middleton, nr. Manchester 
(proposed by Z. Kopal); 

Stephen Piotrowski, University of Warsaw, Poland (proposed by J. 
Witkowski); 

George V. Plachy, 320 Washington Street, Mount Vernon, N.Y., U.S.A. 
(proposed by F. Benario); 

Thomas William Rackham, 14 Holland Street, Cambridge (proposed by 
W. H. Steavenson); 

Alexander William Rodgers, Mount Stromlo Observatory, Canberra, 
Australia (proposed by S. C. B. Gascoigne); 

Francis Roles, 62 Wimpole Street, London, W.1 (proposed by C. E. R. Bruce); 

Trevor Rounthwaite, 52 Gladstone Road, Northcote, Auckland, New Zealand 
(proposed by R. A. McIntosh); 

Peter John Runge, 897 Canterbury Road, Lakemba, N.S.W., Australia 
(proposed by H. Wood); 

Eugene Rybka, University of Wroclaw, Poland (proposed by J. Witkowski); 

Wallace Leslie William Sargent, Dept. of Astronomy, Manchester University 
(proposed by Z. Kopal); 

Richard A. Stevenson, Drury Hills R. D., Drury, Auckland, New Zealand 
(proposed by R. A. McIntosh); 

Robert Griffin Topping, Ardri, 28 Norwood Avenue, Belfast (proposed by 
A. P. Fitzgerald); 

Benjamin Oswald Wheel, Giggleswick School, Settle, Yorks. (proposed by 
W. H. Julian); . 

Dennis Richard Whittington, 8 Cranmore Way, London, N.10 (proposed by 
C. R. Armstrong); 

John Siddons Woodward, 43 Vera Avenue, Grange Park, London, N.21 
(proposed by V. W. R. C. Lionheart-Richardson) ; 

Paul Young, 14 Denholme Road, London, W.g (proposed by R. W. B. 
Pearse). 


Seventy-seven presents were announced as having been received since the 
last meeting, including :— 
M. Migeotte, L. Neven and J. Swensson, The solar spectrum from 2-8 to 23-7 
microns (presented by the Institut d’Astrophysique de l’Université de 
Liége). 
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THE SECULAR ACCELERATIONS OF SATELLITES 
Harold Feffreys 
(Received 1957 August 26) 


Summary 

Maximum effects of tidal friction are calculated for the principal satellites. 
Corrections for elasticity and for imperfection of elasticity are made for 
satellites of the Earth and Mars, based on results for the Earth. Bodily tidal 
friction accounts for about 1 per cent of the secular acceleration of the Moon, 
and for about 10~‘ of that of Phobos, and 10-5 of that of Deimos. Since the 
secular accelerations of JV and Mimas have the wrong sign to be due to tidal 
friction, it appears that bodily tidal friction is not the explanation of any 
observed secular acceleration. 

Other things being equal, the greatest effect of tidal friction on a satellite 
orbit should be for JI. The fact that no secular effect on it has been detected 
indicates that the tidal lag in Jupiter is very small. The result is enough to 
indicate that tidal friction may have had a considerable effect on the orbit of JI, 
and secondary ones on Mimas, Tethys, and Triton, but the effects on all other 
satellites are of negligible importance in cosmogony. 





The apparent secular acceleration of the Moon’s mean motion is well known. 
Recently Sharpless (1945), van Woerkom (1950) and Kozai (1956) have made 
determinations for Phobos, Deimos, JV, and Mimas. ‘Tidal friction in most 
cases produces a secular retardation, and Darwin considered that this might be 
of cosmogonical importance. With Darwin’s data, however, it was impossible 
to derive a time-scale. We now have at any rate a rough idea of the imperfection 
of elasticity in the Earth (Jeffreys 1952, 1957), which can be plausibly generalized 
to the other inner planets and their satellites, and we also know approximately 
the age of the solar system. Consequently it is worth while to attempt a 
quantitative application. Generalization to the outer planets would not be 
justified, since the physical conditions are very different, but the methods used 
are applicable to their satellites also, and comparison of results with observation 
may lead to some information on their imperfection of elasticity. 

The couple due to tidal friction, for a homogeneous planet, is (Jeffreys 1952, 


p. 220, eqn. (13)) 
— a me 
N= 5 nfm — pH sin 2e 


when f is the constant of gravitation, m the mass of the disturbing body, c its 
distance, p the density of the planet, and H and 2e are the maximum height and 
the lag of the semidiurnal tide. For the conventional equilibrium tide (the 
disturbing potential divided by surface gravity) 
3 ma‘ 
=i MA 
where M is the mass of the planet. The actual amplitude is in general different 
from this; in fact (Jeffreys 1952, p. 365, equ. (12)) 
H 5 3 ma* 
~ 2(1 + 19n/2gpa) 4 MA 





where is the rigidity. 





e 
; 
: 
§ 
i 
i 
j 
= 
: 
‘ 
: 





586 Harold feffreys Vol. 117 
If the mean motion of the satellite is m, and its present value is m), we write 
n=nge* 


so that € gives a non-dimensional expression of the amount of change. ‘Then 


N=mn,¢," Tt 
and 


dé _sfm (‘) psin 2€ 


dt 2" Mn\c) 1+ 19u/2gpa 


gmn/fa\>  sin2e 
~ 4M (<) 1+ 19u/2gpa ° 
In Darwin’s treatment the planet is treated as a viscous liquid; in this case the 
elasticity factor is replaced by cos 2, and the maximum rate of change corresponds 
to 2e=}n. For actual planets, however, it is likely that elasticity is dominant. 
In the first place therefore I simply tabulate dé/dt without the elasticity and lag 


factors, which can be restored when required. For representative satellites the 
results are as follows. The unit is 1/1 sec. 


Moon 0°93 X 10716 Mimas 4°0 X10714 Ariel 23x10" 
Phobos —2:6 x 10714 Enceladus 1-6 x 107% Umbriel 1-4 x 107'* 
Deimos 1'o xX 107!" Tethys 44 x20" Titania 8-4 107" 
JV 0-7 xX1074 Dione I'o X 107" Oberon 1'0x 107" 
JI 49 X107% Rhea 270 xX 1075 Triton —6°5 * 10°" 
Jil 16 x<107'4 Titan 0°69 x 107}5 

jill as x10-% Iapetus 0-7 x 107° 

JIV o4 X10~%* 

JVI o8 x 107% 


Where the mass is not directly known I have inferred it from the diameters given 
by Russell, Dugan and Stewart; for Phobos and Deimos the density was taken 
as 3, for JV and JVI equal to that of JI, for the satellites of Uranus as 1. Rhea’s 
mass was taken as 3 x 10-* times Saturn’s. We see at once that the maximum 
possible effect for JVI and Iapetus, and therefore for other distant satellites, 
would be negligible in a time comparable with the age of the Earth, say 10!” seconds, 
and just appreciable for Deimos and Oberon. The largest values are for JI and 
Triton. 

On Darwin’s theory the maximum effect for the Moon would be 0-46 x 107!*/sec. 
Since integration brings in a factor 1/13, the time since the Moon was close to the 
Earth would be 5 x 10’ years, which is practically Darwin’s result. 

Correction for elasticity.—For the Earth the hydrostatic bodily tide, with the 
model used here, has an amplitude 2} times that of the equilibrium tide. The 
actual surface displacement makes the factor about 0-6; thus the elasticity factor 
is about 0°25. 

For Mars, if 4/p=(4°5 km/sec)*, the elasticity factor is 0-06. 

For the outer planets the variation of density with distance from the centre 
is great, and is mostly due to compression. Thus the bulk modulus and gpa 
must be of the same order of magnitude. If the planets are mainly solid p will 
be of the same order. Whether they are fluid or solid we should expect the 
elasticity factor to be larger than for the Earth, and it cannot be more than 1. 
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Phase lags for the Earth and Mars.—We had two possible solutions for the 

Earth on the hypothesis of elastic after-working, namely (Jeffreys 1952, p. 247) 

7’ =I, 7=18 days, 7’ =19°8 days, 
which we shall refer to as short scale damping; and 
: ee I 

Verwre= 30 S jeans 

which we shall call long scate damping. Denote the rate of rotation of the planet 

by w. For both the Earth and Mars the speeds are such that or must be treated 


as large, and 
I = a 
€= ——— |- - —}. 
2a(w—n)\r 7 


This gives the following results :— 





2e dé /dt 
Short scale Long scale Short scale Long scale 
Bot “nis , pl « 21 
Moon 2400 —— 0°97 X 10 0°gO x 10 
oe a _ —19 ai —20 
Phobos 5000 Goes 3x10 3x10 
Deimos ane = % 2x10"! =_ 12x 10-2 
480 5000 


Thus these three satellites are very unlikely to have had their distances 
greatly altered by bodily tidal friction. 
Observed secular accelerations.—We have 


Thus dé/dt is directly calculable from the observed secular acceleration. For 
the Moon (Jeffreys 1952, p. 237, eqn. (11)) 


er 
an? dt 4 
whence 
=) = 10718/sec = 3 x 107!" /year. 


Bodily tidal friction, at the most, accounts for about 1 per cent of the observed 
secular acceleration. 

For Phobos and Deimos B. P. Sharpless (1945) gives the quadratic terms in 
the longitudes 


Phobos (+0°-001882 + 0°-000171)(7'/year)? 
Deimos (-—0*°:000266 + 0°:000164)(7/year)?. 
For JV, A. J. J. van Woerkom (1950) gives 
+ (0°-00034 + 0°-00015)( 7/year)? 
and, for Mimas, Y. Kozai (1956) gives 
+ (0°-000278 + 0°-000119)( 7/year)?. 
All the uncertainties are probable errors. 


' 
; 
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These give for d&/dt 


Phobos —(1°1 x 107!®)(1 + 0°0g)/sec 
Deimos +(0°6 x 107!®)(1 +0°6) /sec 
JV — (0-3 x 1071®)(1 + 0°4) /sec 
Mimas -—-(o°5 x 107'*)(1 +0°4) /sec. 


The observed value for Phobos is about 10,000 times what is likely to be explicable 
by tidal friction. F. L. Whipple and F. J. Kerr, from the data in the second 
edition of The Earth, get a ratio of 800. The discrepancy for Deimos is by a 
factor of 10°. 

The observed values for JV and Mimas have the wrong sign to be accounted 
for by tidal friction. 

The general conclusion is that no observed secular acceleration of a satellite 
is due to bodily tidal friction. 

In any case the secular accelerations for Deimos, JV, and Mimas are 1 to 1°6 
times their standard errors, and therefore somewhat dubious, particularly since 
Whipple and Kerr quote from Clemence a comment that the uncertainties for 
Phobos and Deimos may reasonably be doubled. Thus the effect found for 
Deimos may conceivably have the wrong sign, and all the data would be explained 
qualitatively by a resisting medium or even by the Poynting—Robertson effect. 
Whipple and his collaborators find difficulties with an accretion hypothesis, but 
do not appear to discuss a gaseous resisting medium. 

Tidal lag in Jupiter.—'The great satellites of Jupiter have been observed much 
longer than any of the small satellites considered here, and no secular accelerations 
have been claimed for them. ‘This is noteworthy because the tidal effect on 
dé/dt for JI should be of the order of 70 times that for JV, and with a moderate 
tidal lag it should be easily observable. It is therefore interesting to try the 
consequences of supposing that JI has a secular term in longitude of the order 
of the probable error of those found for the above satellites, say (1°-5 x 10~*) 
(T/year)*. ‘This leads to 


S =(1°4 x 107*)/year =(0°5 x 107*)/sec. 


Comparison with the theoretical value gives 
sin 2€ = 10-4, 
(This includes the elasticity factor, but as this is of the order 1 it need not be 
considered in a rough comparison.) If we adopt the short time scale we get 
ns. 
t 7 400days" 

These results indicate that (1) it is possible (but not proved) that the orbit 
of JI has been considerably enlarged by tidal friction during the history of the 
system (2) the tidal lags must be very small, as if Jupiter behaves nearly as a 
perfectly elastic solid or as a perfect fluid. But if, as seems reasonable, sin 2« 
has similar values for tides raised by other satellites of the outer planets, the effect 


would be of secondary importance for Triton, Mimas and Tethys, and negligible 
for all others. 





160 Huntingdon Road, 
Cambridge : 
1957 August 25. 
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A STUDY OF THE MONOCHROMATIC POLARIZATION 
OF COMET AREND-ROLAND (1956h) 


D. E. Blackwell and R. V. Willstrop 


(Received 1957 August 27) 


Summary 


Measurements of the polarization of the Comet Arend—Roland have been 
made using combinations of glass and gelatin filters isolating respectively the 
continuum near 4530 A and the (0,0) CN emission band near 3850 A. 
Measurements in the continuum extend over a range of phase angle of 26 
and give polarizations between 5 and 20 per cent. There is reasonable 
agreement between the polarization curve of the Comet in the continuum and 
that for a metallic meteorite. Measurements on the CN emission band show 
that the plane of polarization is along the radius vector from Comet to Sun, and 
the amount of polarization is consistent with a fluorescence mechanism. 





1. Introduction.—Although there is now a considerable body of knowledge 
concerning the polarization of light from comets, the data have not so far found 
much application to cometary physics. The reason for this is that almost all the 
available data refer to wide spectral regions. In such regions several quite 
different processes contribute to the luminosity of the comet and the data obtained 
are virtually useless unless these processes can be separated. The stimulating 
work of Ohman (1), however, done with an objective prism, has demonstrated 
the possibilities of near-monochromatic measurements. 

The spectrum of the head of the comet consists of a continuum and emission 
lines of chiefly molecular origin (2). According to current theories the continuum 
is due to the scattering of sunlight by solid particles in the head of the comet*, 
and the line emission is from gas released from the solids when the comet is 
near to the Sun. A study of the continuum alone therefore gives information 
about the solid particles in the comet, and a study of the emission lines gives 
information about the excitation and emission processes. In this paper we present 
data relating to the continuum and the line emission separately. 

2. Previous observations.—Early observations of comets only establish the 
facts that the integrated light is strongly polarized and that the degree of 
polarization varies from one comet to another. The visual data of Wright (3, 4) 
gives polarization of up to 23 per cent for Comet 1881 III and up to 13°8 per cent 
for Comet 18811V. Danjon (5) made visual and photographic observations of 
Comet 1927 VII, obtaining values between 10 and 12 per cent. The first work 
of Ohman (6), that on Comet 1939 ITI, was done visually without filters and gave 
a value of about 6 per cent. Later, with Comets 1940¢ and 1941 c, Ohman used 
an objective prism giving a dispersion at 4200A of 308 A/mm, and attempted 
for the first time to distinguish between polarization effects in the continuum 
and in the line spectrum. His important data will be discussed later in this 
paper. 

3- Observing technique.—The 36-in. telescope of the Cambridge Observatories 
was used at its prime focus. This telescope is of relative aperture f/4-48, giving 


* The origin of the continuum is discussed in more detail in Section 5.2. 
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a plate scale of 50"-3/mm. A large calcite double image prism was placed a short 
distance before the photographic plate. The angle between the ordinary and 
extraordinary rays on emergence from the prism was about 7°, and an image 
separation of 5-5 mm, or 4’-6, could easily be obtained without vignetting. This 
was the image separation used during the observations; a larger linear separation 
is undesirable because of a rapid worsening of image quality, while a closer 
separation leads to an appreciable overlapping of the outer portions of the images. 
The prism was usually oriented so that the planes of polarization of the two images 
were parallel and perpendicular to the line joining the Comet and the Sun. 
The exposures, varying between 7 mins and 247 mins, were never sufficient to 
record more than the head and the brightest part of the tail, so that overlapping 
of the two images was quite insignificant. The investigation was necessarily 
carried out at short notice, and there was no time to construct any device which 
would permit direct guiding on the Comet. The plateholder was therefore 
moved at intervals of between 30 seconds and 20 minutes according to the Comet 
ephemeris, and the position of the Comet relative to the plate checked, usually 
every 15 min. ‘This technique produced rather elongated images of poor 
appearance, but this should not affect the accuracy of measurement of polarization. 





o.eF 
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Fic. 1. Spectral transmission of filter combination, Wratten No. 35 and Chance No. OY 18, isolating 
the continuous spectrum near wave-length 4530 A. 


The Comet was observed in the light of its continuous spectrum only, near 
wave-length 4530 A, by use of a combination of Wratten filter No. 35 and Chance 
filter glass No. OY18 placed between the prism and the plate. The spectral 
transmission of this filter combination is shown in Fig. 1. This filter was used 
with an Ilford Zenith plate, so that the presence of a second peak in the transmission 
curve in the red region is of no importance. The width at half-peak transmission 
is small for a glass—gelatin filter, being only 226A, but the peak transmission is 
only 5-1 per cent. As the Comet was unusually bright, this low transmission was 
adequate, but an interference filter would certainly be more suitable for fainter 
comets. The efficiency of this filter combination in cutting out emission bands 
can be seen by reference to the slit spectrograms of Comet Cunningham given 
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in the Liége Atlas of Comet Spectra (2). On the short wave-length side the 
C,(2,0) emission is completely absorbed, and on the long wave-length side the 
intense C,(1,0) emission with a head at 4737A is greatly attenuated. Details 
of photographs used in the work on the continuous spectrum are given in ‘Table [, 
which gives the exposure time, the distance of the Comet from the Sun, r, the 
distance of the Comet from the Earth, p, and the angle of scattering of sunlight, 6. 
This angle is defined as 7—«, where « is the phase angle. These data, 
together with the ephemeris magnitudes, which may be in error by at least 1", 
are taken from the ephemeris of Candy*. Observations using the double image 
prism were not made before May 10, partly because of cloud and also because 
of the inaccessibility of the Comet to this telescope. 


TABLE I 


Observations of the polarization of Comet Arend—Roland in the continuous spectrum at 4530 A. 








Distance Distance | Ephemeris| Exposure | Angle of Degree of 

Date | Earth-Comet |Sun-Comet} magnitude (min) scattering polarization 
May A.U. A.U. ; per cent 
10 1'O! 0'90 3°6 7 117°8 19°5 

II 1°04 0'92 4°2 7O 120°0 16°9 

13 I*10 096 4°5 120 123°2 10°5 
18/19 1°27 1°07 5°*3 60 130°0 79 

19 1°29 1°08 5°4 150 131°4 7's 

26 1°49 1°21 6-2 247 138°4 S* 





























Observations of polarization have also been made in the radiation from the 
(0, 0) band of CN near 3850A, using a Chance filter OV1 having a maximum 
at 3830A, a width at half-peak transmission of 370 A and a peak transmission 
of 56 per cent. The transmission curve of this filter is shown in Fig. 2. This 
filter was also used with an Ilford Zenith plate. Reference to the Liége Atlas 
shows that it isolates the CN emission from the C, emission near 4050A and 
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Fic. 2.— Transmission curve of Chance filter OV isolating the (0, 0) emission band of CN near 
3850 A. 
* U.A.LC. No. 1585 ; 1957 Feb. 20. 
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the (1,0) emission band of CN near 3590A. ‘There is also an unavoidable 
contribution from the background continuum, but its intensity is falling oft 
rapidly in this region and its total energy, integrated over the pass-band of the 
filter, is probably less than 10 per cent of the CN emission. As measurements 
have been made only on or near the head, there is no need to consider the 
contribution from the N,* molecule. Measurements were made in this region 
for two purposes, (a) to determine the plane of polarization and (5) to determine 
the degree of polarization. The first measurement was made by orientating the 
prism at 45 to the radius vector Comet—Sun, and determining the ratio of 
intensities of the two images. ‘The degree of polarization was then determined 
from a second photograph in the usual way by orientating the prism perpendicular 
to the radius vector Comet-Sun. The details of photographs used in this work 
are given in Table II. 








Tasce II 
Observations of the polarization of Comet Arend—Roland in the (0, 0) emission band of CN near 
3850 A. 
Ephemeris | Exposure | Angle of Purpose of 

Date magnitude (min) scattering photograph Result 

May to determine plane || to radius 

20/21 5°5 180 133°'0 plane of vector Sun—Comet to 
polarization within +5°. 

May to determine head 8-o per cent 

19/20 “9 135 131°°3 degree of tail 6-6 percent 
polarization mean 7°3 per cent 


























Calibrations were provided for all of these photographs by processing at 
the same time photographs of a step wedge which had been taken under identical 
conditions in the laboratory. 

4. Reduction of observations.—The plates were scanned with the Cambridge 
Observatories microphotometer using a scanning spot approximately 
o-I mm xo-1 mm (5”x5"). The degree of polarization is then defined as 


_h-I, 
at es % 


where /, and /, are the intensities of the two images. But two kinds of systematic 
error are involved in deriving the degree of polarization from the microphotometer 
traces. One is that the double image prism may, through a defect in its 
construction, divide an unpolarized beam into two unequal beams. The other 
is that the prism when used in an f/4:5 beam gives two images having slightly 
different foci and having slightly different structures. These errors were 
investigated by taking photographs of the planet Jupiter through the double 
image prism. Jupiter is a good test object for this purpose because it has 
approximately the angular diameter of the head of the Comet, and its integrated 
light has a polarization of less than 0-5 per cent (7). When these images were 
scanned with the microphotometer it was found that they showed a polarization 
of less than 0-5 per cent (corresponding to a value of A//J of less than 1 per cent), 
provided the scans were integrated over the whole image size. As, therefore, 
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in the Liége Atlas of Comet Spectra (2). On the short wave-length side the 
C,(2,0) emission is completely absorbed, and on the long wave-length side the 
intense C,(1,0) emission with a head at 4737 is greatly attenuated. Details 
of photographs used in the work on the continuous spectrum are given in Table I, 
which gives the exposure time, the distance of the Comet from the Sun, r, the 
distance of the Comet from the Earth, p, and the angle of scattering of sunlight, @. 
This angle is defined as s—«, where « is the phase angle. These data, 
together with the ephemeris magnitudes, which may be in error by at least 1™, 
are taken from the ephemeris of Candy*. Observations using the double image 
prism were not made before May 10, partly because of cloud and also because 
of the inaccessibility of the Comet to this telescope. 


TABLE I 


Observations of the polarization of Comet Arend—Roland in the continuous spectrum at 4530 A. 








Distance Distance | Ephemeris| Exposure | Angle of Degree of 

Date | Earth-Comet |Sun-Comet| magnitude (min) scattering polarization 
May A.U. A.U. 2 per cent 
10 ered 0-90 3°6 7 117°8 19°5 

II 1°04 0°92 4°2 70 120°0 16°9 

13 I"10 0°96 4°5 120 123°2 10°5 
18/19 1°27 1°07 3 60 130°0 7°9 

19 1°29 1°08 5°4 150 131°4 7°5 

26 1°49 1‘21 6-2 247 13874 . 





























Observations of polarization have also been made in the radiation from the 
(0, 0) band of CN near 3850A, using a Chance filter OV1 having a maximum 
at 3830 A, a width at half-peak transmission of 370A and a peak transmission 
of 56 per cent. ‘The transmission curve of this filter is shown in Fig. 2. This 
filter was also used with an Ilford Zenith plate. Reference to the Liége Atlas 
shows that it isolates the CN emission from the C, emission near 4050A and 
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Fic. 2.—Transmission curve of Chance filter OV isolating the (0, 0) emission band of CN near 
3850 A. 
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the (1,0) emission band of CN near 35g90A. ‘There is also an unavoidable 
contribution from the background continuum, but its intensity is falling off 
rapidly in this region and its total energy, integrated over the pass-band of the 
filter, is probably less than 10 per cent of the CN emission. As measurements 
have been made only on or near the head, there is no need to consider the 
contribution from the N,* molecule. Measurements were made in this region 
for two purposes, (a) to determine the plane of polarization and (b) to determine 
the degree of polarization. The first measurement was made by orientating the 
prism at 45° to the radius vector Comet-Sun, and determining the ratio of 
intensities of the two images. The degree of polarization was then determined 
from a second photograph in the usual way by orientating the prism perpendicular 
to the radius vector Comet-Sun. The details of photographs used in this work 
are given in Table II. 








TasBLe II 
Observations of the polarization of Comet Arend—Roland in the (0, 0) emission band of CN near 
3850 A. 
Ephemeris | Exposure | Angle of Purpose of 

Date magnitude (min) scattering photograph Result 

May to determine plane to radius 

20/21 5°5 180 133 '0 plane of vector Sun—Comet to 
polarization within +5°. 

May to determine head 8-o per cent 

19/20 5*4 135 131° °3 degree of tail 6-6 percent 
polarization mean 7°3 per cent 


























Calibrations were provided for all of these photographs by processing at 
the same time photographs of a step wedge which had been taken under identical 
conditions in the laboratory. 

4. Reduction of observations.—The plates were scanned with the Cambridge 
Observatories microphotometer using a scanning spot approximately 
o-I mm xo-1 mm (5”x 5"). The degree of polarization is then defined as 

I,-I, 

ae eo 
where /, and /, are the intensities of the two images. But two kinds of systematic 
error are involved in deriving the degree of polarization from the microphotometer 
traces. One is that the double image prism may, through a defect in its 
construction, divide an unpolarized beam into two unequal beams. The other 
is that the prism when used in an f/4:5 beam gives two images having slightly 
different foci and having slightly different structures. These errors were 
investigated by taking photographs of the planet Jupiter through the double 
image prism. Jupiter is a good test object for this purpose because it has 
approximately the angular diameter of the head of the Comet, and its integrated 
light has a polarization of less than 0-5 per cent (7). When these images were 
scanned with the microphotometer it was found that they showed a polarization 
of less than 0-5 per cent (corresponding to a value of AJ/J of less than 1 per cent), 
provided the scans were integrated over the whole image size. As, therefore, 
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the photograph of the Comet in the continuous spectrum showed virtually no 
tail, the degree of polarization was found by integrating the measured intensity 
over the whole of the image recorded on the plate (usually about 45” x 45”). 
The photographs in the light of the CN emission band on the other hand showed 
a considerable tail and with these a different procedure was necessary, Two 
separate regions have been examined for polarization, each of them large in 
comparison with the aberrations of the image, one in the head, of area 50” x 35”, 
the other in the tail, about 2’ distant from the nucleus, and of area 50” x 30”. 
Imperfect guiding on the Comet may give rise to another systematic error 
which would not be shown in the tests on Jupiter because these use short exposures. 
In addition to having the effect of reducing the effective exposure time of the plate, 
imperfect guiding also means that each part of the plate is exposed to a varying 
intensity. ‘This means at least that the exposure time of the calibration plate 
must be reduced by an uncertain amount. However, tests have shown that the 
shape of the characteristic curve of the plate does not change sufficiently with 
exposure time to affect the results appreciably, and we estimate that the effect of 
varying intensity will not appreciably affect the comparison of two closely similar 


images. 
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ANGLE OF SCATTERING 


Fic. 3.—Polarization in the continuous spectrum as a function of angle of scattering. 


Results for the continuous spectrum are shown in Table I and those for the 
CN emission band in Table II. In the determination of the plane of polarization 
of the emission CN lines of May 20/21, when the prism was orientated at 45° 
to the radius vector Sun—Comet, the intensities of the two images were found 
to be equal within the limits of experimental error. This means that the plane 
of polarization is within 5° of the perpendicular to this direction. All the values 
of polarization given in Tables I and II are unlikely to be in error by more than 
+ 1°5 per cent. 

5. The polarization in the continuous spectrum.—All the data obtained for the 
continuous spectrum near wave-length 4530A are plotted against angle of 
scattering, 9, in Fig. 3. In drawing a smooth curve through these observations 
it has been assumed that the polarization is zero for = 180° (back scattering), 
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5.1. Comparison with the observations of Ohman.—The only other polarization 
data of this kind are those of Ohman (1) obtained for the Comet 1941¢, using 
a low dispersion objective prism. In interpreting the results Ohman assumed 
that the spectrum of this Comet is continuous, and measurements were made at 
wave-length 4300A and 4700A on two consecutive nights. Reference to the 
McDonald spectrum given in the Liége Atlas (Plate XVIII 62) shows however 
that there are some emission lines in the spectrum of this Comet, and measurements 
at the two wave-lengths quoted will be influenced by the polarization in these 
emission lines, particularly in the C,(1, 0) band at 4700 A. The mean polarization 
found by Ohman for the two nights and the two wave-lengths is 22 per cent. 
But as a measure of the polarization of the continuum this value is fairly certainly 
too low, owing to the contribution of the C,(1, 0) band, which itself probably had 
a polarization of only about 8 percent. The correction needed is uncertain, but if 
it is supposed that the C,(1, 0) band supplies 18 per cent of the measured light, 
which is not unreasonable, the continuum polarization should be increased to 
25 per cent. This value has been plotted on the graph and it will be seen that 
it is not incompatible with the present observations; so our curve has been 
extended to include it. The complete curve of Fig. 3 is likely to be a faithful 
representation of the true polarization curve, except that the polarization near 
68=go° should perhaps be even higher than is shown. 

The large variation of polarization from one comet to another, found by the 
earlier observers who used large bandwidths, is now easily explicable in terms of 
a continuum polarization and an emission polarization, the relative properties 
of continuum and emission depending upon the comet. The observations of 
Wright on the rapid fluctuations of polarization in Comets 1881 III and 1881 1V 
are probably due to rapid changes in the relative proportions of continuum and 
emission. It is relevant that at the time of these observations the geomagnetic 
index of the Sun was high (8). 

5.2. Origin of the continuous spectrum.—Before the observations can be 
interpreted in terms of a comet model a clear decision is needed about the origin 
of the continuous spectrum, currently supposed to be due to the scattering of 
sunlight by solids in the head of the comet. This supposition is apparently based 
on the observation (for example, that of Mayall (g)) that the continuum spectrum 
is similar to that of the Sun; it shows Fraunhofer lines and has a distribution 
of intensity with wave-length which is similar to that of the solar spectrum*. 
But none of this evidence eliminates the possibility of scattering by a low 
temperature electron gas. Let us suppose, now, that 80 per cent of the light is 
due to scattering by dust, and the remainder to scattering by free electrons. 
Then, when 6=7/2, the electrons are scattering completely polarized light, 
and most of the observed polarization of 25 per cent would be due to the electrons, 
which would contribute to the whole a polarization of 20 per cent. Evidently 
an estimate of the importance of electron scattering in the formation of the 
continuous spectrum is desirable. In an attempt to decide this question we shall 
first calculate the order of magnitude of the electron density that would be required 
to account for the continuous spectrum of Comet Arend—Roland. 

Suppose that the Comet is a cube of side # and that the total flux from the Sun 
into unit solid angle is Bo. Then the flux on to the Comet from the Sun is 
Bo x h*/r*. 

* Measurements by Dr. Liller of Michigan show that the spectral distribution in the continuous 
spectrum of the tail differs from that of the solar continuum (private communication). 
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If the total number of electrons in the Comet is N,, the number per unit 
cross sectional area is N,/h?. Hence, the total flux of radiation scattered by the 


Comet in all directions 











h* WN, Bo 
=Bo x 7 x re xXxO,= ~ x N,o,, 
where a, is the Thomson scattering coefficient. 
, BoN.o, 
Hence flux on unit area of Earth from Comet= at F 
711" p 
; ‘ : Bo ‘ : 
The flux on unit area of Earth from the Sun= R2 where R, is the distance of 
the Sun. 
H integrated brightness of Comet in continuum __ &,? i BoN.o, 
_— integrated brightness of Sun Bo 4m1*p* 


. 

According to Wellman* the Comet had a very strong continuum on 1957 
April 23. On this date the ephemeris magnitude was o™-8 and we shall not be 
far wrong if we take the magnitude of the Comet, in the light of its continuous 
spectrum only, to be 2™. 

xy », 

Hence, qari? N,o,~ 6x 10°? 
and N,~ 4°3 x 10%, 
Most of the light of the Comet in the continuous spectrum was contained in an 
image of diameter 50”, or 2:2 x 10° cm. 

Hence the average electron density required to account for the whole of the 
continuum is 


| a heey ee} 
€ ~~ 10°6 x 107° 


= 4°06 x 10" cm-, 

According to Richter’s (10) model, the head of a typical comet contains between 
104 and 10° observable molecules cm-*. If, therefore, we assume that each 
molecule is singly ionized, and that the number of unobservable molecules does 
not too greatly exceed the number of observable molecules, the maximum electron 
density is only 10°cm-*, which is a factor 10° too smal! to account for the 
brightness of the continuous spectrum. We conclude therefore that scattering 
by free electrons does not make an appreciable contribution to the continuum, 
and that the observed polarization is due to scattering by solid particles. 

5.3. Interpretation in terms of a comet model.—The detailed interpretation 
of this polarization curve is not easy. It is likely, not only that we are dealing 
with a mixture of different kinds of particles, but also that the physical character- 
istics of the mixture change with distance from the Sun. We emphasize, too, 
that we have no knowledge whatsoever of the size or of the nature of the surface 
of the cometary particles. They may have a smooth glossy surface, or a deeply 
pitted surface, or a finely powdered surface. The polarization depends to a certain 
extent on the size and the kind of surface, and a comparison with laboratory data 
must be made cautiously. We suppose here that the particles, at least in the head 
of the comet, are large in comparison with the wave-length of light. 

* U.ALIL.C. No. 1596. 
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A long laboratory investigation of the polarization produced by reflection 
from various materials has been carried out by Lyot* (7). These materials were 
all of terrestrial origin and included inorganic crystals, sedimentary rocks and 
igneous rocks. In making a comparison between this kind of data and the cometary 
data it is unfortunate that comet polarization values are available only for 
scattering angles greater than 100°, but it will be clear from Fig. 3 that there is 
probably a maximum in the curve near 0=go°. This fact eliminates all the 
materials examined by Lyot because all that show a sufficient polarization near 
4 = 100° also have a steep gradient of polarization at this angle of scattering, and 
show a maximum polarization near #= 40°. 

















\ 
25} ~~ 
~ 
\ 
VN 
\ \ 
20F \ \ 
3 \ CONTINUOUS SPECTRUM 
. \ OF COMET 
a or } 
Nish METALLIC & | 
« METEORITE * | 
< 
} 
a 
w 
% 10 
< 
- 
z 
w 
& 
w o7 AYU 
a 
Sk ' ! 
SNOW ™ 
ow 
oO iL iL = 
°o 
80° 100° 120° 140° 160° 180 


ANGLE OF SCATTERING 


Fic. 4.—Polarization curves for (a) Comet 1956 h, (b) metallic meteorite material, (c) snow. The 
distance of the Comet from the Sun is also shown. 


Kloverstrom and Rense (11) have examined the polarization of light scattered 
by meteorite materials. Most of these show a maximum near 0=90°, the value 
of the maximum polarization being at least 25 per cent. The curve shown in 
Fig. 4 for a metallic meteorite (specimen “‘O” of Kloverstrom and Rense) is 
taken from their data: it reproduces the observations fairly closely. The stony 
meteorites show a polarization which is too great at all angles of scattering, and 
especially near 06=150°. We therefore draw the tentative conclusion that, 
at the time of observation, the continuum was chiefly due to scattering by 
meteorite material of a metallic nature. At the same time we do emphasize 
the inevitable uncertainty of this conclusion. 

It is evidently of great interest to compare these data with those appropriate 
to the Whipple model of a comet (12, 13). Essentially, this model supposes a 
comet to consist of a conglomerate of dust and a sublimate composed of condensed 
compounds such as CHy, CO,, NH3, C,N, but mainly H,O. Far from the Sun 
the dust particles are embedded in sublimate and the appearance of the comet 

* Lyot’s polarization values are expressed in terms of an “‘ angle de vision”’ which is equal to 
the phase angle «, or 7— 0, where @ is the angle of scattering. 
42 
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head is presumably due to reflection from the sublimate. Near to the Sun the 
sublimate melts, producing a cloud of dust and gas. The gas causes the familiar 
line spectrum and the continuum is now due to reflection by dust as well as by 
any remaining sublimate. According to Whipple the H,O ice of a comet is 
melted at about r=2A.U. from the Sun. A polarization curve for ice obtained 
from Lyot’s data* is also shown in Fig. 4, and clearly the polarization of Comet 
Arend-Roland is chiefly due to scattering from solids other than ice at least to 
a distance of 1:21 A.U. Measurements of polarizationwhen a comet is about 2A.U. 
from the Sun should provide a crucial test of the Whipple model especially if 
the comet is well placed to give the minimum scattering angle of 150°. Unfor- 
tunately it was not possible to observe Comet Arend—Roland at a phase angle 
greater than 138°. 

6. The polarization in the (0, 0) emission band of CN.—Ohman’s work on the 
polarization in the spectrum of Comet Cunningham (1940c) has already shown 
that the (0, 0) CN emission band is polarized to the extent of 10-1 per cent at 
an angle of scattering of go°. The theoretical polarization cannot be calculated 
with precision because it depends upon an accurate knowledge of the intensity 
distribution within the band. In particular, the lines of low J-number should 
have very high polarization. Ohman estimates a theoretical value of about 
8 per cent, but an examination of the intensity distribution in this band in typical 
cometary spectra in the Liége Atlas shows a large contribution from the lines 
of low J-number indicating that a rather higher value, perhaps even 12 per cent, 
may be expected. Using the expression given by Ohman for the variation of 
polarization with angle of scattering, 

Py sin? 6 
Po= T+ p,, cos®d 
90 
we find that our mean value of 6-6 per cent is equivalent to a value of 13-7 per cent 
at 8=go°. The agreement with Ohman’s value and the expectation of theory is 
satisfactory. 

The chief purpose of this work, however, was to measure the plane of 
polarization of the CN band which could reasonably be expected to be along the 
radius vector to the Sun (i.e. the electric vector is perpendicular to this direction). 
This is not necessarily the case, for we may recall that measurements by Lyot 
(14, 15) of the polarization of Ha and the D, line of Hei in prominences have 
demonstrated that the plane of polarization may deviate from this direction by 
as much as 30° for Ha and 40° for D3. These anomalies have been attributed 
by Thiessen (16) to the presence of a magnetic field. In Comet Arend—Roland 
our measurements show that the plane of polarization is parallel to the radius 
vector to the Sun to within +5°. We conclude therefore that the magnitude 
and direction of polarization of this band is consistent with the simple theory of 
fluorescence with no appreciable complicating factors. 

Acknowledgment.—The authors are much indebted to Professor R. O. Redman 
for discussion of the problems involved in this work. 


The Observatories, 
Madingley Road, Cambridge: 
1957 August 26. 


* Lyot gives four curves for dry snow of varying texture ; the average of these four curves, which 
are quite similar, is shown in Fig. 4. 
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MODEL SOLAR ATMOSPHERE DEDUCED FROM OBSERVATIONS 
OF LIMB DARKENING 


Antoni Przybylski 
(Received 1957 August 6) 


Summary 


A method of correcting model solar atmospheres by comparing the observed 
limb darkening and the intensities for the centre of the solar disk with the same 
quantities computed theoretically from the model solar atmosphere is des- 
cribed. The differences between the observed and computed quantities are 
represented as linear functions of the errors in the adopted values of the 
absorption coefficients and of the variational corrections of the temperature 
distribution at a number of chosen mean optical depths. The solution of the 
system of equations obtained in this way gives the corrected values of the mono- 
chromatic absorption coefficients and the corrected temperature distribution of 
the model solar atmosphere. 

The investigations show that : (i) the observations of the limb darkening 
may not be free of systematic errors; (ii) the determinations of the intensities 
for the centre of the solar disk are relatively poor, which makes it impossible to 
obtain reliable corrections for the absorption coefficient; (iii) the structure of 
the solar atmosphere is affected to a certain extent by convection; (iv) large 
errors in the values assumed for the absorption coefficient may be counter- 
balanced by relatively small changes in the temperature distribution; and 
(v) it seems to be impossible to explain all the results of the observation of solar 
phenomena without assuming temperature inhomogenities in the structure of 
the solar atmosphere. 





1. Introduction.—In a former paper by the present author (1), which will 
subsequently be referred to as Paper I, a model solar atmosphere has been con- 
structed by computing the effect of blanketing on the structure of another 
*‘initial”” model which had a practically constant net flux of energy, disregarding, 
however, the spectral lines. This ‘‘ corrected” model given in Table II of Paper I 


TaBLe I 
Differences between the observed and computed limb darkening for \=6980A. 
O-C O-C 

sin } initial model corrected model 
00836 —0°0034 —0°0035 
0°2226 + 0°0003 0*0000 
0°3693 +0°0015 +0°0006 
0°5165 +0°0075 t+-0*0060 
* 0°6531 +0°0094 + 0°0075 
0°7621 +0°0098 + 00085 
0°8369 +0°0096 + O°O102 
0°8835 +0°0084 +0o°o121 
09166 +0°0063 +O°O145 
09462 +0°0017 +0°0179 
09658 —0'0076 +0:'0200 
©°9775 —0°O173 +0°0233 


0°9857 —0°0344 +0°0242 





Wwe 


No. 6, 1957 Model solar atmosphere deduced from observations of limb darkening 601 


should come close to the real solar atmosphere if the adopted values of the absorp- 
tion coefficient are correct and if the effect of convection on the structure of the 
solar atmosphere is negligible. A comparison of the observed limb darkening 
of the Sun with the limb darkening computed from the model shows, however, 
considerable discrepancies, sometimes even larger than for the initial model 
without blanketing ; e.g. Table 1 shows the differences between the observed 
and computed values of the limb darkening for the initial and the corrected model 
for the wave-length A=6980A. ‘The large discrepancies between observation 
and computation show clearly that the two sources of error mentioned above 
are not negligible and their investigation is the subject of the present paper. 

The mean differences between the observations of the limb darkening made 
for a certain wave-length by different observers are of the order of 0-008. The 
series of observations used in the present investigations seems to be, however, 
more precise ; this can be expected from the fact that the differences between 
the observed and computed values of the limb darkening, shown in Table I 
for one wave-length, are not scattered at random but show a clear systematic 
run. ‘There can therefore be no doubt that these differences exceed the experi- 
mental errors of the observation and that appropriate corrections to the temper- 
ature distribution of the model solar atmosphere can reduce the discrepancy 
between the observed and the calculated limb darkening. Indeed, the final 
results of the present investigations show that the mean error of one observation 
is of the order of 0-003, i.e. considerably less than most of the differences shown 
in Table I for A= 6980A. 

2. General remarks on the limb darkening.—The intensity I,(o, #) of the solar 
radiation for the frequency v emerging at the angle # to the normal to the solar 
surface is given by the integral equation 


1,(0,8)= | ” B,(r,)e~ 7°? sec Pdr, (1) 


where B,(7,) is the source function and 7, the monochromatic optical depth. 
The introduction of 4 =cos # into equation (1) gives 


1,(0,n)= | By(r,Je-"*" dr, | (2) 
which divided by /,(0,0) leads to the equation of limb darkening 
ee 
1,{0u4)/1,(0,0)=4.(u)= [Ao = dry (3) 


The limb darkening ¢,() and the intensity /,(0, 0) in the centre of the solar 
disk are known from observations, and consequently equation (3) gives the 
possibility of investigating the temperature distribution in the solar atmosphere, 
which cannot be done directly. The analysis of the structure of the solar atmos- 
phere by means of equation (3) is, however, a difficult task since the results are 
seriously affected by observational errors and since the contribution of the deeper 
layers of the atmosphere to the integral on the right-hand side of equation (3) 
decreases rapidly on account of the factor e 7””. 

3. Solution of the equation of limb darkening.—Many attempts have been 
made by various authors to solve equation (3) ; all the solutions can, however, be 
grouped in three classes described below. 
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3.1. Inversion of the Laplace transform.—As equation (3) shows, the function 
pu. $,(u) is the Laplace transform of the function B,(7,)/1,(0,0). The latter can, 
therefore, be found from the inversion of its Laplace transform : 


ra+t 


B,(z,)/1{0,0)= [ewan de. (4) 


The integration in equation (4) is to be performed along an axis parallel to 
the imaginary axis or along this axis if a=o. This procedure assumes that the 
function ¢,(u) of the limb darkening exists in the whole plane, which is not true, 
since for physical reasons it is defined only in the interval from o to 1 along the 
real axis. A theoretical error is made, therefore, if one represents the observed 
values of the limb darkening by a function and extends its validity to the whole 
plane. The source function found in this way depends on the analytical form 
of the function used in order to represent the observed limb darkening. 

The above procedure has been used by Busbridge (2) who represented the 
limb darkening curve by the function 


$(u) = {m/p + (1 —m)5* 
with two disposable parameters m and a>o. ‘The source function can then 
be expressed by ['— functions. 

The method described by Sykes (3) will be seen to belong essentially to the 
same group if the limb darkening curve be represented by a linear combination 
of the powers of log p. 

3-2. Methods using analytical expressions for B,(r,)/I,(0,0).—Equation (3) 
can be solved, assuming that B,(7,)/J,(0,0) can be represented by a function 
with a number of disposable parameters if the Laplace transform of this function 
is known. The solution consists then in the computation of these parameters 
from the system of equations (3) written for a sufficient number of the values 


of ». Milne (4) and Lundblad (5), for instance, expressed the source function 
in terms of the form 


a-— 


a,7,* 
which lead to terms of the form 
a,.k! pk 
in J,(0, #). This method has been used by Lundblad (5) and later by Chalonge 
and Kourganoff (6) and by Barbier (7) for the investigation of the temperature 
distribution of the solar atmosphere. 

It is also possible to expand the source function in a power series or to represent 
it by means of the exponential integral functions or other functions. Examples 
of such procedure were published by Pierce and Aller (8) and by Mineur and 
Peyturaux (9). 

The disadvantage of these methods is that they determine very poorly the 
source function for large values of r. This is due to the fact that the deeper 
layers of the solar atmosphere contribute only little to the emerging radiation on 
account of the factor e~”” under the sign of the integral in equations (2) and (3). 
Attempts to express the source function in an analytical form hide this fact 
completely. 

3-3- Numerical methods.—In analytical methods the source function is 
represented by various functions, each of which gives values for every level of 
ty In numerical methods the source function is defined only for certain 
discrete optical depths or is treated as constant over certain intervals. 
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Thus, for instance, Plaskett (10, 11) divided the photosphere into layers 
and assumed the source function to be constant within each layer. A layer 
(7,,i Ty, 441) in which the constant value of the source function is C, ; gives a 
contribution 

Cy,de~ 7#—e- Tom) 
to the emerging intensity J,(o,4). This means that the problem of solving 
equation (2) is reduced to the solution of a system of linear equations with the 
unknowns C, ;. This method is fairly simple but has a big disadvantage in that 
the system of linear equations is often indeterminate ; this is due to the fact that 
little error is made if a constant is added to a C, ; in one interval and the same 
constant is subtracted in the next interval to C, ;,, (or C, ;,)—provided that 
both intervals are equal. In order to obtain a reasonably well-determined 
solution it is necessary to restrict the number of unknowns C,, but then the 
determination of the source function is very rough. 

Other methods are based on the use of some quadrature formulae for the 
evaluation of the integral on the right hand side of equation (3) by means of a 
linear combination of the values of the integrand at some, preferably equidistant, 
points. Every observation is thus expressed in terms of some discrete values of 
B,(7,,;) which are to be found by solving the system of equations (3). The 
most advantageous seems to be the use of quadrature formulae of the Newton- 
Cotes type applied to a finite interval, for instance from 7, =0 to 7, =3°0 and the 
approximation of B,(7,) by a linear function for 7,>3-0. ‘This method has been 
used by Sykes (3). 

An attempt was also made by Sykes (3) and Mineur and Peyturaux (9) to 
evaluate the integral on the right hand side of the equation of limb darkening by 
means of the Laguerre polynomials, but it failed completely since it led to a source 
function which certainly does not correspond to the actual temperature distri- 
bution in the solar atmosphere. 

4. Method used in the present investigations. 

4.1. Introductory remarks.—The methods described in the preceding section 
give results necessarily affected by observational errors in the determination 
of the monochromatic intensities /,(0,0) for the centre of the solar disk as well 
as by the errors in the determination of the limb darkening. They must be 
applied to every wave-length separately and this means that the results are 
seriously affected by systematic errors. For this reason the present investigations 
make use of a new numerical method which allows the possibility of combining 
observations in different wave-lengths and reducing the effect of errors in the 
determination of central intensities. 

In this method the differences ‘‘O-C”’ between the observed and computed 
limb darkening are represented as a linear function of the errors in the adopted 
values of the monochromatic absorption coefficients and of the corrections to the 
temperature distribution at certain optical depths. The solution of these 
equations then gives the corrected values of the absorption coefficients and the 
corrected temperature distribution. 

4.2. Computation of the limb darkening.—For the purpose of computing the 
limb darkening the monochromatic depths have been computed for the following 
wave-lengths in A: 3646+, 4235, 5050, 6252, 8206—, 8206+, 9214, 10503, 
12213, 14588—, 14588+, 16031, 17790, 19984, 22794—, 30392, 45589 and 
91185. The monochromatic absorption coefficients used in these computations 
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were taken from the author’s own tables, which differ only slightly from the 
values represented in diagrams by Vitense (12). Subsequently the limb 
darkening has been computed in these wave-lengths for the following values of 
sin # : 0°0, 0:2, 0°3, 0°4, 0°5, 0°55, 0°60, 0°65, 0°70, 0°75, 0°80, 0-825, 0°85, 0°875, 
0°90, 0°92, 0°95, 0°96, 0°97, 0°98 and 0-99. ‘The table of the values of the limb 
darkening obtained in this way was sufficient to allow the interpolation of the 
value of limb darkening for all intermediate values of A and #. 

The values of J,(0, 0) and J,(0, 1) have been computed by numerical evaluation 
of the integrals in equation (2) using Cotes’ formulae. ‘The computations were 
made to five figure accuracy for reasons explained in Section 4.4., although the 
last digit at least has no real significance. 

4.3. Comparison with observations.—'The computed values interpolated to the 
proper wave-length and sin? were compared with the observations of limb 
darkening made by Pierce, McMath, Goldberg and Mohler (13). ‘This series 
of observations was chosen for comparison because it was hoped to use its deter- 
minations of limb darkening in the infra-red for which no other determinations 
were available. In addition, these observations have been used already by 
Pierce and Aller (8) in their investigations of the solar atmosphere, which makes 
it possible to compare their results with those obtained in the present 
paper. 

The results of the comparison of these observations with theory are shown 
in Table I for one wave-length (A= 6980A) in order to give an idea as to the order 
of magnitude of the discrepancies between theory and observation. 

4.4. The effect of errors in the absorption coefficient.—In order to find the 
effect of a systematic error in the adopted values of the absorption coefficient 
the theoretical limb darkening has been computed once more with the mono- 
chromatic absorption coefficients increased by 10 per cent. The differences 
between the values of the limb darkening computed in this way and those obtained 
in Section 4.2. will be designated by «, in further investigations ; they represent 
the change in the limb darkening due to a 10 per cent increase in the mono- 
chromatic absorption coefficient. In further investigations it will be assumed 
that the actual error in the monochromatic absorption coefficient is k, times 
10 per cent, i.e. its total effect on the error “‘O—C” of the computed limb darkening 
is x,.k,. ‘This means that it is assumed that the monochromatic absorption 
coefficient for a certain wave-length A is wrong by a constant factor (1 +k,/10) 
throughout the whole photosphere. ‘This assumption can hardly be considered 
as an ideal one, but any other more complicated assumption would make the 
computation too difficult or even impossible. For this reason the corrections 
of the absorption coefficients obtained in the present investigations must be 
considered as fairly rough. 

The values of the coefficients «, are small and therefore five figure accuracy 
was necessary in the computation of limb darkening. Examples given in Section 
4-6. give an idea of the order of magnitude of these coefficients. 

4.5. Corrections to the temperature distribution.—In order to make allowance 
for the change in the limb darkening due to the change in the temperature distri- 
bution the model solar atmosphere has been divided into 7 layers by planes 
corresponding to the mean optical depths 7=0, 0:0593, 0174, 0°4*7, 0-984, 
2°145, 3°508 andoo. (The values of the mean absorption coefficient wsed in the 
computation of the mean optical depths are shown in Table I of Paper I.) It 
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has also been assumed that in every interval 7; . . . 7;,, the reciprocal temperature 


5040 
— 
has to be changed by a certain number n; of units of the second decimal place 
in order to obtain the true temperature. ‘These numbers n; are the unknowns to 
be found by the comparison of the computed limb darkening with the observations. 

Suppose that the change in the source function over the interval 7; . . . 7;,, 
corresponding to a change of one unit in the second decimal place in the value of 
is AB, ; ; this change, AB, ,, is of course variable over the interval but a good 
mean value can be obtained if one chooses the value corresponding to the point 
situated for instance } down through the layer. This choice of representative 
points gives more weight to the outer layers of each interval, which, of course, 
contribute more to the emerging intensity. Actually the representative points 
were chosen so that the reciprocal temperatures © corresponding to them were 
numbers which could be expressed exactly with only two figures following the 
decimal point. In the last (infinite) interval the representative point was 
chosen at the mean optical depth t=4-665. The mean optical depths of the 
representative points are : T=0-015, 0°0852, 0°2281, 0°581, 1°441, 27491, 4°665 
for which the reciprocal temperatures are © = 1-20, 1-01, 0°93, 0°86, 0°76, 0°69 
and 0-61. 

Assuming now that for the frequency v the interval r,,; . . . 7,,;,, corresponds 
to the interval 7; . . . 7;,, of the mean optical depth, we find that the contribution 
to the emerging intensity in the direction » due to the change of the reciprocal 
temperature by n; units of the second decimal place is 


Ti+ I _, ar 
n,AB, i| e7 '# aa a n,AB, fe —Ty,i/u —e~ Tr, i+] “). 
mn , 


O= 


The corrected emerging intensity /,*() is then 

T,*(u)=1(u) +n AB, (e~ 7 ie —e7 74+ Un) (5) 
where /,() is the emerging intensity for the model to be corrected. Dividing 
this by 

"(w= 1)=1,(1)+n,AB, (e~™'-e ~~ +1), 
and disregarding the terms of higher order, we obtain the following equation 
for the corrected limb darkening 











= I,*() =i T,(#) AB, ; a— Tp tly —Ty §4-1/ I,(u) 
Oy (H)= Ta) = T(1) +l Te ee 
(e~ 7», t/u —e~ Ti+ Mn) >, 


Using the designation A,, for the expression between the brackets and 
$,(u)=1,()/1,(1) for the limb darkening of the uncorrected model we obtain 


$,*(u)=d(u)+A, ;. 0). 


‘Taking into account the contribution of the corrections applied to all the intervals 
in question we obtain the generalized equation 


bw) =$.u) + SA, om (6) 


The coefficients A, ; were computed for the same values of A (or v) and © for 
which the limb darkening was computed (Section 4.2). The values of A, ; 
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corresponding to the values of A and © for which observations were available 
were obtained by interpolation. 

For simplicity, the fact that the intensities J, were computed for 7=0 has not 
been indicated explicitly in the formulae of this section. 

4.6. Equations of condition.—Incorporating the effect of a change in the 
absorption coefficient on the limb darkening into equation (6) we obtain 


$,*(u) -_ $() +k,k, + 2A vj 
or 
Kk, + Ay im +A, ong+ A, gng+A, mgt A,, 5M5 + A, gtgt+ A, qn, 


= $,"(u)— $,()- (7) 


If the right hand side of this equation is equal to the difference between the 
observed and computed limb darkening this equation can be used for the deter- 
mination of the corrections k, and n; Equation (7) can be written for every 
observation and since the number of equations is greater than the number of 
unknowns the system of equations (7) must be solved using the method of least 
squares. 

In order to convey to the reader an idea as to the order of magnitude of the 


coefficients of equation (7) this equation is shown below for three observations 
made for the wave-length 


A= 6980A and sin # = 0-2226, 0°7621 and 0-9658 : 
+ 22Reos9— 2%,—- Jmg— 17M,— 16m,+ Qnz;+ 21mM,+ 17N,;= 0°0000 ; 
+321Rg999— 44M, — 141N,— 281n,— 167, + 271N, + 270N, + 139N,= +0°0085 ; 

+ §26Reog9 — 186, — 472, — 528, + 164n, + 723; + 291M, + IL IN, = +0°0200. 

The coefficients on the left hand side of the above equations are given in 
units of the fifth decimal place. The wave-length A is used instead of the 
frequency v as a subscript of the unknown k’s. 

The method described above bears a certain similarity to the numerical 
method of Plaskett (Section 3.3), which assumes that the source function is 
constant in each interval. Here, however, the source function is not constant, 
and only the correction of the source function superimposed on an approximate 
initial solution is considered to be constant in some intervals. The results of the 
solution of the system of equations (7) can, however, be represented finally in the 
form of a smooth continuous mean curve. In comparison with Plaskett’s method 
the present method seems to be more advantageous since it allows us to combine 
the observations in different wave-lengths for the purpose of the determination 
of the temperature distribution. The determination of the monochromatic 
absorption coefficients is, however, more difficult than the determination of the 
temperature distribution, since the observations in different wave-lengths cannot 
be combined for this purpose. For this reason it is easy for systematic errors 


to affect seriously the determination of the k,’s. 
The method described above is likely to give good results if the expected 


corrections are reasonably small, i.e. if the investigated model solar atmosphere 
is a reasonably good one. 


5. Numerical results. 


5-1. First solution The observations of limb darkening by Pierce and his 
collaborators (13) were made in 13 wave-lengths. In all, 168 observations were 
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made, which gave 168 equations with 20 unknowns (thirteen k,’s and seven 
n;’s). ‘The system of these equations proved to be very poorly determined and 
therefore the results were seriously affected by observational errors. The 
values found for the n;’s did not show a systematic “‘ run”’, but were scattered at 
random. This means that the method employed here suffers from the disadvan- 
tage mentioned in Section 3.3. in connection with Plaskett’s method. 

It seemed therefore necessary to introduce further equations which could be 
obtained from the comparison of the observed and computed intensities for the 
centre of the solar disk. 

Taking into account the effect of a change in the temperature distribution 
in each interval, we can write equation (5) in the form 


1,*(1)=1,(1) + Un AB, (e~ 7. i—e7 7, i+), 


It remains to include the effect of a change of the monochromatic absorption 
coefficient in this equation. Using the designation 6/,(1) for the change 
of the central intensity due to a 10 per cent increase of the absorption coefficient 
(which is known from previous investigations) we obtain 

1,*(1)=1,(1) +R,61,(1) + Un AB, (e7 7 i—e7 ™i+1), (8) 
Dividing this equation by /,(1) for reasons explained below we obtain 


51, yy oo reste 1M(1)—11) 





7,1) 7,() T,(1) 
Using the designation D, for 5/,(1)/I,(1) and E, ; for the content of the braces 
we obtain 

Dk, + XE, n= (1,*(1) — 1,(1)}/1,(2)- (9) 


The coefficients D, were computed from the results of previous calculations. 
The coefficients E,, ; were computed additionally in the same way as the coefficients 
A, ; have been computed in Section 4.5. 

The right-hand side of equation (7) is equal to J,*(u)/J,*(1)—J,(w)/7,(1), 
whilst the right-hand side of equation (9) is {/,*(1)—J,(1)}/J,(1). It can easily 
be seen that the right-hand sides of both equations represent values of similar 
order of magnitude and this is the reason why equations of the type (9), rather 
than of the type (8), are used with the system of equations (7). The use of equa- 
tions (8) would give too much weight to observations of the intensities at the 
centre of the solar disk. 

The observations of Pierce (14) and Peyturaux (15), interpolated for the 
appropriate wave-lengths, were used in the present investigations. An additional 
constant —1°264 was applied to the relative measurements of Pierce in order 
to adjust them to the absolute values of Peyturaux. Since no observations 
were available for the far infra-red the last three wave-lengths (A= 35 000A, 
83 000A, 102 o00A) were not taken into account for these and the following 
investigations. For A=5485A the value of the central intensity given by Mulders 
(16) was used since Peyturaux’s observations do not extend as far into the visible 
region of the spectrum. 

For purposes of illustration the following equation for the correction of the 
central intensity at A= 6480A is shown : 


— 4214999 — 642, — 232m. — 6383 — 1003n, — 1356n. — 504M, — 191N, = +0°0882. 
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The coefficients on the left-hand side of this equation are again given in units 
of the fifth decimal place. 

Ten new equations were thus added to the 130 equations of limb darkening ; 
the number of unknowns is now 17 (ten k,’s and seven m,’s). 

An attempt to solve this system of equations again resulted in failure since 
the corrections m, for the deeper layers were scattered at random. It was therefore 
found necessary to introduce the further quite reasonable assumption that the 
unknowns n; could be represented by a polynomial 

n;=a+ Br,+y7;, 2<1<7 (10) 
where 7, are the mean optical depths of the representative points mentioned in 
Section 4.5. In this way it was hoped that the values obtained for m, would not 
be scattered at random. Since, however, the previous solutions showed that 
the correction n, for the first interval is large in comparison with the other values 
n, this correction was left unaltered in the equations. If now the new values of 
n;, expressed in terms of the unknowns «, f and y, are introduced in the equations 
of condition, the number of unknowns is reduced to 14 (ten k,’s, «, B, y and m,). 
The three equations quoted in 4.6. and the one equation given in this section 
are shown below in the new form : 


22 Rego — 2%,+ ‘Jat 1308+ 513y= 0°0000 ; 
321... — 44m, + giat+15388+ 51g2y= +0°0085 ; 
526 ... —186n,+ 28gx+22198+ 5749y = +0°0200 ; 
— 4214 — 64n,— 3924 — 48488 — 10476y = + 0°0882. 


The least squares solution of these equations gave the results shown in the 
second column of Table I]. In addition to the unknowns found in the solution 


Tasie II 
Results of the solutions 
Unknown Ist solution 2nd solution 3rd solution 4th solution 
Rsaso +442 +1°84 —3'84 +1°03 —rag +1°48 —5°38 
Reiss +1°33 +1°86 —4°06 +1°04 —1'89 +1°54 — 7°05 
Revso +058 +1°89 —495 +1°05 —2°25 +1°60 —7°88 
Roses +097) +1°84 —4°49 +103 —0°97, £1°55 —4°99 
Roo20 +2°43 +1°89 —4°01 +1°06 —o'70 +1°56 —5*69 
Riss = + 3°64 $199 —2°96 +1°05 +o'16 +1°58 —4°64 
Risser = + 5°69) 42°39 +289 +1°29 +415 +1°91 +or72 
Rinses = + 5°13, 2°53 +4°74 £1°37 +5°07 42°04 + 2°40 
Roees +6°87  +2°49 +3°46 +1°28 +5°25 +1°94 ° +1°28 
Reser «= + 6°81 = + 2°39 +037 +1°29 +3°35 +1°93 — 1°60 
a —2°771 +1°220 +1°519 +0°719 —0'93I +1°012 — 1688 
B —o'8809+1°1548 =+1°3457+0°5964 +0°3529+0°9119  +2°8375 
Y +0°4216+0°2277. —0-0105+0°1175 +0'1846 +0°1799 —0'2747 
ny —10°682 +3042 -—7-658 +1'589 —8-822 +2:352  —8-822 
Ne —2°843 +0°889 +1634 +0°527 —o'g00 +0°736 —1°428 
Ns —2'950 +0°907 +1°825 +0°536 —o'841 +0°751 —1°O55 
Ng —3°141 +0°958 +2'297 +0°560 —0°664 +0°791 —0°132 
Ns —3°165 +1-092 +3°437 +0°628 —0°039 +0900 +1°831 
Ng —2°349 +1°245 +4'806 +0°713 +1°094 +1°022 +-3°676 
n> +2296 +1°712 = +7°569 +0°985 +4°734 +1°380 + §°§70 


Each unit in the solutions for the j’s indicates an increase of the absorption coefficient of 10. 


~~ = while each unit in the solutions for nj indicates a change of the reciprocal temperature 
y +o-0r, : 
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the values of m,, computed with the aid of equation (10), are shown there together 
with their mean errors. The other columns of this table will be explained in the 
next section. 

5.2. Second solution.—The investigations of the preceding section show 
clearly that the results are seriously affected by any observational errors. There 
is no doubt that our knowledge of the intensities for the centre of the solar disk 
is poor and therefore further investigations were undertaken in order to get an 
idea as to the reliability of the results. 

In a solution, in which six of the equations for the central intensities were 
neglected and only the first four were taken into account, results were obtained 
which are shown in the third column of Table Il. They differ a great deal from 
those obtained in the first solution, especially in the determination of the un- 
knowns k,, where the differences between the results are often larger than the 
probable errors. ‘This means, of course, that the equations for central intensities 
must contain some systematic observational errors of a fairly large order of 
magnitude, a fact which after all is known from other investigations (3, 8). There 
is, of course, no reason to suppose that the four equations retained in the second 
solution are better than the other six disregarded in the computations, and therefore 
the second solution cannot be regarded as better than the first one. 

The fact that the omission of only six equations out of the hundred and forty 
causes a serious change in the results of the solution can easily be explained if the 
equations for the central intensities are seriously affected by observational errors. 
As the examples quoted in the preceding section show, the coefficients of the 
equations for the central intensities are larger than those of the equations of 
limb darkening. In fact, by the addition of ten equations for central intensities 
to the system of equations of condition some coefficients of the normal equations 
are increased by more than twenty times and this explains why errors in the 
observations of central intensities affect the solution so seriously. 

5.3. Third solution.—It is not possible to obtain a solution without taking into 
account the equations for central intensities, which on the other hand may 
seriously affect the solution if they contain systematic observational errors. In 
order to minimise the effect of these errors on the results another solution was 
made in which, however, the equations for central intensities were taken into 
account with the weight +, only. In this way the contributions of the equations 
for central intensities and for the limb darkening to the coefficients of normal 
equations are held in better balance. ‘The results of this solution are shown in 
the fourth column of Table II, together with their mean errors. 

5.4. Discussion of results —The results of the third solution seem to be more 
reliable than those of the first two solutions. In order, however, to check how 
trustworthy they are, another solution was made in which only the equations 
for the limb darkening were taken into account and the value n, = 8-22, obtained 
in the third solution, was accepted for m, in advance. The results of this solution 
are shown in the last column of ‘Table II. 

The corrections of the temperature distribution obtained in this solution do not 
differ greatly from those obtained in the third solution. The same, however, 
cannot be said about the &, corrections and this means that observations of limb 
darkening also contain systematic errors, which consequently affect the solution. 

The systematic errors of the observations are more likely to affect the results 
for the k, than those for the corrections of the temperature distribution. This is 
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due to the fact that the corrections for the temperature distribution are present in 
all the equations of condition whilst those for the absorption coefficients are 
present only in 12 or 13 equations, i.e. in a number of equations equal to the 
number of observations at a given wave-length. The results of the solutions 
of the equations of condition confirm that the corrections of the temperature 
distribution are better determined than the k, corrections. 

The results of the third solution seem to be more trustworthy than those of 
other solutions. The corrections of the temperature distribution obtained in 
this solution will be taken as the basis for the computation of the final model. 
Little can be said about the reliability of this solution as far as the corrections to 
the absorption coefficients are concerned except that their determination is less 
trustworthy than the mean errors would suggest. The results obtained in the 
third solution would decrease the mean absorption coefficient by 2-6 per cent 
for © = 1-0 and by 47 per cent for@=0-8. The results obtained for the corrected 
temperature distribution must also be accepted with some reservation. The 
formulae used for the computation of the coefficients A,,; are valid provided that 
the assumed values of the monochromatic absorption coefficients are not seriously 
wrong, which, however, is not certain. 

The mean errors of one single observation obtained in the three solutions 
are 0:00685, 0:00328 and 0-00494 respectively. After the smoothing of all 
observations by means of the least squares method (third solution) the remaining 
errors for the central intensities in order of the increasing wave-lengths are : 
— 13%» — 41%, — 12%, +2°B%, — 21%, —2°8%, +4°5%y» +8-3%, +4°0%, 
—0o°5%. These remaining errors are much higher than the mean error of one 
observation obtained in the solution (0°-00494) which shows clearly that the 
results of the investigations could be improved considerably if our knowledge of 
the distribution of the intensity in the solar radiation were better. 

6. Comparison with the results of Pierce and Aller.—The observations used 
in the present investigations have already been used by Pierce and Aller (8) in 
their investigations on the solar atmosphere in which however another method 
was used. Their paper was mainly concerned with the determination of the 
absorption coefficient from the observational data and their results are very much 
the same as those of the present investigations. The results are seriously affected 
by errors in the determination of the intensities for the centre of the solar disk. 
The values obtained from observations differ considerably from those found 
theoretically, sometimes to such an extent that this disagreement must be con- 
sidered asincredible. The present investigations show that the effect of systematic 
errors in the observations is even more serious than was assumed hitherto, 
because large errors in the determination of the absorption coefficients may be 
compensated by relatively small errors in the temperature distribution. 

7. Adopted final model solar atmosphere.—The present investigations show 
that the determination of monochromatic absorption coefficients from observa- 
tional data gives very uncertain results and consequently no reliable mean 
absorption coefficient can be deduced. Hence, for the purpose of computing 
the final model solar atmosphere, it is assumed that the mean absorption coefficient 
remains unchanged. As already stated in Section 5.4, this assumption is fulfilled 
within a few per cent by the results of the third solution: The corrected temper- 
ature distribution obtained in this solution is also accepted as the basis for the 
computation of the final model. For the mean optical depth +>0°0593 the 
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corrections of the temperature distribution are given by formula (10). In the 
interval from t=0 to r=0-0593 the corrected temperature was found by extra- 
polating graphically the temperature distribution obtained for 7> 0-0593 so that 
the integrated temperature correction corresponds to a uniform change of the 
reciprocal temperature © by —0-0882 found in the third solution for the outer- 
most interval. 

The procedure of computing model stellar atmospheres has been described 
in previous papers of the author (1, 17) and the same method was used here. 
The results of the computations are shown in Table III. The mean optical 
depth 7 in terms of the mean continuous absorption coefficient (taken from Table I 
of Paper I) is shown in the first column of this table ; a uniform increase of the 


Taste III 
Adopted final model solar atmosphere 


(p indicates total pressure due to hydrogen alone) 


T a To Y is log p log p, 

° ° ° 4598° 

o'ol o’oll 0'0092 4649 3°87 1°97 
0°02 0°022 00185 4701 4°03 o'12 
0°03 0033 0°0276 4751 4°13 o'21 
0°04 0°044 0°0368 4799 4°20 0°29 
0°05 0°055 0'0460 4852 4°25 0°35 
0°06 0-066 O°0551 4907 4°30 0°40 
0°07 0°077 0°0642 4959 4°33 0°45 
0°08 0-088 0°0733 5010 4°37 0°49 
0°09 0099 0°0825 5063 4°39 0°53 
o'10 O'110 00916 5114 4°42 0°57 
o'12 0°132 0°1099 5204 4°46 0°63 
O'14 O°154 0°1284 5272 4°50 0°69 
0°16 0176 0°1468 5329 4°53 0°73 
o'18 0198 o°1651 5376 4°56 0°77 
0°20 0°220 0°1833 5418 4°59 o'81 
0°25 0°275 0'2289 5505 4°64 0°89 
0°30 0°330 0°2747 5579 4°68 0°96 
0°35 0°385 0°3202 5644 4°72 1°02 
0°40 0440 0°3657 5705 4°75 1°08 
0°45 0°405 0°4113 5762 4°77 1°13 
0°5 0°55 0°4566 5818 4°79 1°18 
0°6 0°66 0°5471 5926 4°83 1°28 
o'7 0°77 0°6374 6028 4°86 1°37 
o'8 0°88 0°7275 6125 4°88 1°46 
°o'9 0°99 08182 6216 4°90 1°54 
1‘o I'Io 0*9090 6302 4°92 1°61 
1°2 1°32 I‘OgIo 6469 4°95 1°75 
1°4 1°54 1°2724 6607 4°97 1°88 
1°6 1°76 1°4528 6739 4°98 1°99 
1°8 1°98 1°6317 6859 5°00 2°09 
2°0 2°20 1°8088 6966 5°O1 2°17 
2°5 2°75 2°2451 7194 ; 5°03 2°35 
3°0 3°30 2°6732 7368 5°04 2°48 
3°5 3°85 3°0967 7500 5°05 2°57 
4°0 4°40 3°5163 7589 5°06 2°64 


4°5 4°95 3°9347 7654 5°07 2°68 
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absorption coefficient by 10 per cent shown in the next column allows approxi- 
mately for the change of the mean absorption coefficient due to the absorption 
lines. ‘The third column of the table gives the monochromatic optical depth 
for A= 5000A. 

The final model solar atmosphere is shown graphically in Fig. 1 together with 
the initial model, which has a constant net flux of energy if spectral lines are 
disregarded, and the corrected model, which takes into account the effect of 
blanketing on the structure of the solar atmosphere. The graphical representation 
of all three models shows clearly that convection has an effect on the structure 
of the solar atmosphere in that it lowers the temperatures in the deeper layers of 
the photosphere. This effect seems to be a fairly well established fact since it is 
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Fic. 1.—Comparison of the initial, corrected, and final models: 
———————_ Initial model 
— ——— — Model corrected for blanketing 
— tet Final model 


shown by all solutions of the system of the equations of condition. The 
lowering of the temperature in the deeper layers of the photosphere is compensated 
partly by its rise in the outer layers and partly by the decrease of the absorption 
coefficient in the visible region of the spectrum, where the intensity of the solar 
radiation is strongest. Generally speaking one can say that convection and blank- 
eting effect modify the structure of the solar atmosphere in opposite directions. 

The present investigations show that the sharp drop of temperature in the 
outermost layer of the atmosphere deduced from the theory of the blanketing 
effect and from the observations of the limb darkening of certain spectral lines 
(1, 18, 19) is not substantiated by the limb darkening in the continuum. This 
indicates that the structure of the solar atmosphere is more complex than is 
usually assumed in the computation of model solar atmospheres. Indeed, there 
can be no doubt that there are temperature inhomogenities in the solar photo- 
sphere and they may be responsible for some contradictions in the results 
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derived from observations. It is intended in a subsequent paper to investigate 
the possibility of allowing for these inhomogenities in the present investigations. 

Fig. 2 shows a comparison of the finally adopted model solar atmosphere, 
with Barbier’s model (20) derived from the limb darkening in the continuum, 
de Jager’s model VII (18) derived from the limb darkening of the hydrogen lines 
and Vitense’s model I (19) derived from the investigation of metal lines. 
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Fic. 2.—Comparison with empirical models: 

— — — — — Véitense’s model I 

—--— se de Fager’s model VII 

—:-—: -—: -— Barbier’s model 
Model under discussion 
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Summary 


Results of self-consistent field calculations with exchange for the ground 
state (3s)? (3p) and the first excited state (3s) (3p)* of Fet+!* are presented. 
The differences between energy levels of the terms are given, including the 
separation due to spin-orbit interaction. Integrals are given for use in the 
calculation of transition probabilities for lines of importance in the solar 
corona. 





Introduction.—Some of the emission lines in the spectrum of the solar corona 
have been identified by Edlén (3) as due to transitions from metastable states of 
highly ionized Fe. The structure of these ions is important in the development of 
a theory to explain the physical conditions in the solar corona such as relative 
abundancies and the formation of metastable states through collision. Calcula- 
tions without exchange have been done previously (Gold, (§)) but calculations 
with exchange, besides being more accurate, should allow a more detailed analysis 
of transition probabilities. ‘The results would also be useful in the calculation of 
collision cross-sections and in the interpolation of wave functions for atoms with 
the same electron configurations. 

Method of solution—The method used for determining the self-consistent 
field solutions of the Hartree-Fock equations was similar to that described in a 
previous paper (4). The only essentially different feature is that the (3s) (3p)? 
configuration gives rise to several terms, namely *P, *D, ?P, 2S, each with its 
own set of differential equations. ‘The general form of the Hartree-Fock 
equation assuming €,, ,,;=0 is 


i W(d+1 2 
[ StF EVV OHA} ens SE | Pease) + 2X (sr) 0, (1) 


Minimizing the energy of each term (1, p. 192) determines Y (nl; r) and X (nl; r) so 
that, in the usual notation (6), 
Y (38;r)=2Y (15, 18) +2 Yo (25, 25) +6 Yo (2p, 2p) + 2 Yo (3p, 3p) 
X (3857) = Yo (ts, 38) P(18) + Yo (28, 38) P(2s) + Y; (2p, 38) P(2p) 
+a Y, (38, 3p) P(3p) 
Y (3p;7r)=2Y (15, 18) +2 Yo (2s, 25) +6 Vo (2p, 2p) + Yo (35, 35) 
+ Yo(3p, 3P)— 4 ¥2(3p, 3p) 
X (3p 517) =(1/3) Yi (1s, 3p) P(1s) + (1/3) Yi (2s, 3p) P (2s) + Yo (2p, 3p) P(2p) 
+ (2/5) Y2(2p, 3p) P(2p) + (a/2) Yi (3s, 3p) P (35), 
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where the coefticients a and 6 depend on the term and are given in ‘Table I. ‘The 
ground state configuration (3s)? (3p) gives rise to only one term, the *P, and 
coefficients in the differential equations are the same as for a complete group. 


TABLE I 
Coefficients in Hartree-Fock equation 
Term i 2D ? aS 
a 2/3 1/3 1/3 1/3 
b 1/5 — 1/25 1/5 2/5 


In all the calculations it was assumed that the core wave functions remained 
the same as those for Fe*+!*(4) and so only the (3s) and (3p) wave functions had to 
be made self-consistent with each other. ‘The wave functions self-consistent to 
0-0002 are given in Tables V and VI. 

Term energies.—The Hartree-Fock equations are derived by applying the 
principle that the energy is stationary with respect to small perturbations of the 
wave functions. Since the variation between terms of the wave functions for the 
outer electrons is not large, the variation of the core wave functions will be small 
and so we may assume that the total energy of the core is the same for all terms, 
say Ey. By substituting for /(nl) in terms of €,, ,,,, F;,(nl,n'l’) and G,, (nl, n'l’) in 
the expression for the total energy of the system, we obtain the following equations 
for the term energies * : 


(3s)? (3p) 7 : E- Ey oe €3s, 38 ‘ (1/2) €3p, 3p ~ Fy (35, 35) 
— 2F 4 (35, 3p) + (1/3) Gi (35, 3p) 


(38) (3p)? *P 
2D | . E- E,= — (1/2) €35, 3s — €3p, 3p — 2 Fy (3s, 3p) 
71 — Fy (3p, 3p) + 4 Fy (3p, 3p) + a G, (35, 3p) 
3c CI 
= 


From Table II it is seen that the ?P level is slightly higher than the 2S. _Calcu- 
lations without exchange where €,, 3, and €3,, 5, are the same for all terms would 
predict that the *S level was the higher. Comparing these energy differences 
with the observed ones for Ca*’, it is seen that the 2P level should be the higher one 
and that the separation of the levels is not great enough. In these calculations, 
superposition of configurations has been neglected and could quite possibly 
account for this discrepancy. 








TaBLe II 
Contributions to term energies (atomic units) 

(35)? (3p) r (3s) (3p)? ‘ 

Term sy 4p 2p 2p 8S 
Eas, 39 30°770 31608 317046 29°910 31°05! 
Esp, 3p 28552 28-912 28°176 28:076 27°494 

F (38, 35) 1‘280 
Fo(3p, 3p) 1297 1'292 1291 1:288 
F3P, 3p) 0°665 0662 0:660 0659 
F4(35, 3p) 1°287 1°291 1°286 1°280 1°284 
G (35, 3?) ©°535 0°536 0°534 0°532 0°533 
- ae —48°722 — 48-105 -47°426 — 46-928 -46°961 
° 0617 1'296 1°7 1°761 
_ E—E(P) one : 

obs. in Cat? ° o'581 0°769 1°47 0°074 





* The notation is slightly different from that of (x). Here the D factor (1, p. 177) has been 
incorporated into the coefficients instead of the F,, and G,. functions, 


a 
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Spin-orbit interaction.—The electrostatic interactions break a configuration 
up into terms each with a given L and S. When spin-orbit interaction is included 
levels are characterized by their J values and the eigenvalues S and L are no 
longer accurate quantum numbers (1, p. 193). In the case of Russell-Saunders 
coupling where spin-orbit interaction is small compared to electrostatic interaction, 
a given term characterized by S and L will split into a small group of levels each 
with a different J value. If I (ySL./) specifies a level then, by the Landé interval 
rule, 

VP (ySLJ)—V(ySLJ —1)=J C(ySL). 
For the configurations considered here, {(ySL) is directly related to ¢;,, (1, pp. 122, 
199). If 
P' (3p 57) + (2U (7) —€— (1+ 1)/r*) P(3p)=0 (2) 
then 


> eat ae d 
Ss) = — Ra? | 7; (3p) FG U(r) dr. (3) 


For wave functions which satisfy the Hartree-Fock equation (1), the potential 

term, U (r), consists of two parts (7), which for the (3p) wave function are 

(1) the coulomb potential [N — Y (r)]/r 

(ii) the exchange potential X (3p; 1r)/rP(3p;r). 
Physically, the exchange potential can be interpreted as a correction for exchange 
between an electron at radius r in a (3p) wave function and the other electrons in the 
same spin state. If the probability of finding the (3p) electron at r is zero, that is 
P(3p;r)=0, there would be no exchange, so the nodes of X(3p;r) would be 
expected to coincide with those of P(3p;r). On account of the physical approxi- 
mations of the theory on which the numerical work is based, they may not coincide 
exactly. However, the factor P?(3p;7r) in the integrand in (3) ensures that the 
whole integrand remains non-singular at a node of P(3p;7); indeed this expression 
for C;,, can be written 


“a /( 


2 2 > , , 

bay = Ra? | L {5 W- Y)+ $ y’- (5) + = far, (4) 
Values of ¢,, as computed from the wave functions are given in Table III where 
Ra? was taken to be 5851 cm~! according to the values of R and « given by DuMond 
and Cohen (2). ‘The energy difference between the ?P,,. and ?P,,, levels of the 
ground state is (3/2)¢;,. From Table III this calculated energy difference is 
18,750cm~! which compares well with the observed difference of 18,852°5 cm='; 
the value calculated without exchange is 19,390 cm (5). 


Tasie III 
Values of 3, 


Term 
(35)* (3p) ? 12,500 cm! 
(3s) (3p)" ; 12,500 
2D 12,460 
sp 12,500 
2S 12,400 


Transition probabilities.—'Vransitions between states are accompanied by the 
emission or absorption of radiation. All lines due to electric dipole radiation arise 
from transitions between states of opposite parity and may be thought of as resulting 
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from single electron jumps such that /= +1. Three transitions (35)(3p)* > (35)*(3p) 
are allowed and their transition probabilities depend on the radial integral, 


| r P(nl;r) P(n'l—1;1)dr. 
The value of this integral is 0-696 for ?P, 2S >?P ; 0-695 for ?D — *P; ando-7o1 for 
wave functions without exchange (5). 

Transitions which are forbidden for dipole radiation may occur due to a 
quadrupole moment though the lines will be much weaker. ‘The transition 
probability depends on the radial integral, 


“oO 
| r? P(nl;r) P(n'l ;r) dr, 
0 
the values of which are given in Table IV. ‘These are to be compared with 0°565 
for wave functions without exchange. 

Transitions due to a magnetic dipole moment may also occur but their tran- 
sition probabilities depend only on the angular momentum of the particles 
concerned and not on their radial distribution. 














Taste IV 
Values of : r= P (nl; r) P(n'l’; rv) dr. 
0 
(35)*(3p) (3s)(3p)* 

Term 1p “p "=D “pp ag 

sy 0°552 

4p 0°550 

*D 0°555 0°556 ©°557 

P 0°556 0°558 

25 °°559 

TABLE V 
Normalized 3s wave functions for Fe*™ 
P(3s) 
r (35)*(3p) (35) (3p) 
2p ‘p 2D 2p 2S 

0°0025 00816 0:0826 00819 0°0807 00819 
0°00375 0°1186 ©°1199 o°1190 O°1172 01190 
©0050 0°1529 0°1546 0°1535 O'1S11 ©°1535 
0°0075 0°2144 0°2167 oO’2151 0°2118 oO'2151 
0’O10 0°2669 0°2696 0°2676 0'2635 02676 
O°O15 0°3473 0°3510 0°3484 0°3429 0°3434 
0°020 0°3992 0°4034 0°4004 0°3941 0°4004 
0°025 0°4269 0°4315 0°4282 0°4216 0°4283 
0°030 0°4344 0*4390 0°4358 0°4289 0°4358 
0°035 0°4250 0°4295 0°4263 0°4196 0°4263 
0°040 0°4016 0*4060 0°4028 0°3964 0°4030 
0°045 0°3668 0°3709 0°3680 0°3621 0°3680 


0050 0°3230 0°3267 0°3241 03188 0'3242 
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TAaBLe V (continued) 
Normalized 3s wave functions for Fe+™ 





P(3s) 
r (35)°(3p) (3s) (3p)? 

2p sy 2D 2p 2S 
0°06 0°2161 0°2187 0'2168 0°2130 0°2168 
0°07 0°0934 0°0949 0°0938 0°0gI9 0°0938 
0°08 —0°0349 —0°0348 —0°0348 —0°O351 —0'0348 
0'09 ~0'1619 01630 ~0'1622 —0'1605 —0o'1622 
o'10 —o0°2821 —0°2843 —0°2827 —0'2792 — 02827 
orl! —0°3920 0°3954 —0°3929 -0°3879 —0°3929 
o'l2 04890 9°4935 0°4903 0°4837 —~0°4904 
0°13 °°5719 —0°5771 0°5734 0°5656 ~0°5734 
O°14 0°6398 -0°6458 0°6415 0°6326 —0°6416 
O'ls 0°6923 -0°6990 —0°6943 —0°6845 —0°6944 
O°175 0°7596 0°7674 0'7619 0°7507 0'7620 
0*200 0°7445 0°7526 0°7471 0°7356 —0°7471 
0°225 —0°6625 —0°6705 0°6651 —0°6541 -0°6653 
0°250 0°5304 0°5379 0°5330 ~~ 0°§232 —0°5330 
0'275 0°3639 0°3703 0°3662 0°3580 0°3662 
0°300 0°1765 o'1819 0°1787 0°1723 -0°1785 
0°325 00201 00160 0°0183 0°0226 0°0183 
0°350 0'2168 O'2141 0°2155 0'2176 O°2155 
0°375 0°4066 0°4053 0°4056 0°4057 0°4057 
©°400 0°5843 0°5843 0°5836 05818 0°5838 
0°45 08899 0°8923 0'8g00 08849 o'8go02 
0*50 I‘rigi 1°1234 1*1200 11125 I‘1201 
o's5 1‘2711 1°2764 1'2726 1°2639 1'2726 
0°60 1°3531 1*3586 1°3548 1°3463 1*3548 
0°65 1°3765 1°3814 1°3780 1°3708 1°3780 
0°70 1°3535 1°3570 1°3547 1°3496 1°3547 
0'75 12961 1°2979 1°2968 1°2944 1°2967 
080 1°2149 1°2148 I‘2151 1'215§7 1°2151 
°o'9 1°0156 I‘O117 1°O147 1'0210 1°0146 
I‘o 08072 08002 08052 0'8157 08051 
I'l 0°6174 06088 06148 06276 0°6148 
I'2 0°4584 0°4492 0°4555 0°4691 0°4555 
¥*3 0°3322 0°3235 0°3294 0°3423 0°3204 
I*4 0°2359 0°2283 0'2335 0°2449 0°2336 
1°5 0'1648 0°1584 0°1628 0°1723 0°1628 
16 O°1135 0°1083 o-rr118 o'1196 o-r118 
1°7 0°0772 0°0732 00760 00819 0'0760 
18 0°0520 00490 00510 0°0555 00510 
19 0°0346 0°0323 0°0339 0°0374 0°0339 
2°0 0°0229 0°0212 0°0224 0°0249 0°0224 
2°25 0°0079 0°0072 0°0076 00089 0°0076 
2°50 00026 0°0023 00026 0°0031 00026 
2°75 00009 00008 00008 o'oor! 00008 
3°00 00003 00002 070003 00004 0°0003 
3°25 ©°0001 00000 00000 0°0001 00000 
3°50 00000 00000 ©0000 00000 00000 
E35. 38 30°770 31608 31°046 29°910 31051 
ao(35) 34°910 35°280 35°018 34°476 35°023 
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0'0025 
0°00375 
0°0050 


0°0075 
o’olo 


O'O15 
0'020 
0°025 
0°030 
0°035 
0°040 
0°045 
0°050 


0°06 
0°07 
0°08 
009 


II 


13 
14 
15 


ooooo 


"175 
0°200 
0°225 
0°250 
0°275 
0°300 
0°325 
0°350 
0°375 
0*400 


0°45 
0°50 
0°55 


0°65 
0°70 
0°75 
o'8o0 
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Tas_e VI 


Normalized 3p wave functions for Fe+™ 


(35)*(3p) 


2p 


0'0014 
00032 
00056 


00124 
0°0214 


00450 
0'0751 
0*1 100 
01485 
0°1895 
0°2321 
0°2753 
0°3188 


0°4033 
0°4822 
0°5528 
0°6133 
0°6628 
0*7008 
0°7274 
0°7425 
0°7466 
0°7405 


0°6850 
0°5820 
0°4437 
0°2820 
0'1068 
—0°0735 
—0'2521 
— 0°4237 
—0'5846 
—0°7317 


—0°9784 
—1'1581 
—1'2733 
— 1°3309 
— 1°3405 
—1°3120 
—1°2552 
—1°1784 
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‘p 


00014 
0°0032 
0°0056 


0°0124 
0°0214 


0°0450 
0'0751 
0°1099 
01485 
0°1895 
0'2321 
0°2753 
0°3188 


0°4034 
0°4823 
0°5530 
0°6136 
0°6631 
0°7013 
0°7278 
0°7433 
0°7475 
0°7417 


0°6866 
0°5840 
0°4462 
0°2847 
0°1099 
—0'0703 
—o'2489 
—0°4207 
—0°5817 
—o'7291 


—0°9769 
—1°1§78 
—1'2740 
— 1°3324 
— 1°3425 
—1°3142 
—1°2§72 
—1°'1801 


(38) (3)’ 
2D i 
00014 00014 
0°0032 00032 
00056 070056 
00124 0°0124 
00214 0°0214 
0°0450 0°0450 
0°0749 0°0751 
©1097 O'1100 
0°1482 0°1485 
0°1892 0°1895 
0°2316 0°2321 
0°2750 0°2754 
0°3183 0°3188 
0°4027 0°4033 
0°4814 * 0°4822 
0°5519 0°5527 
06122 0°6131 
06616 06625 
06996 0°7004 
07260 0°7268 
0°7410 0°7417 
0°7451 0°7457 
0°7390 0°7396 
0°6834 0°6837 
0°5804 0'5804 
0°4424 0°4422 
02809 02806 
0*1060 0°1057 
—0°0739 —0°0740 
—O°2519 —0'2519 
—0°4231 —0°4228 
—0°5833 —o'5829 
—0°7300 —0°7292 
—0°9758 —0'9746 
— 1°1§49 — 11533 
— 1°2697 — 12679 
— 1°3274 —1°3255 
—1°3374 —1°3358 
— 1°3097 — 1°3084 
—1°'2538 —1°2531 
— 1°1783 —1°1780 


0°0014 
00032 
0°0056 


0°0124 
00212 


0°0449 
0°0748 
0°1096 
0°1479 
01888 
0°2312 
0°2744 
0°3175 


04018 
0°4802 
0°5504 
06105 
0°6598 
0°6975 
0°7236 
0°7385 
0°7425 
0°7361 


0°6802 
0°5769 
0°4388 
0°2773 
O°1025 
—0'0772 
—0°2550 
—0°4257 
—0°5853 
—0°7312 


—0°'9755 
— 1°1§30 
—1°2665 
— 1°3232 
— 1°3330 
— 1°3057 
—1°2506 
—TI 1763 
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TABLE V1 (continued) 


Normalized 3p wave functions for Fe+' 





P(3p) 
r (35)*(3) (38) (3p)? 

t sainitie ’ 

2p 4p 2D 2p aS 
o'9 —0°9937 —0°9940 0°9955 09964 09964 
1°0 0°7999 —0°7988 —0°8035 08051 —0o'8072 
I'l —0°6215 —0°6194 06261 —0'6282 ~0°6317 
1°2 —0°4695 —0°4668 —0°4745 0°4767 —o°4811 
1°3 —0°3467 -0°3436 —0°3516 ~0°3537 —0°3584 
1°4 —0'2512 —o0°2481 —o'2556 0°2574 —o'2621 
1°5 —0°*1790 —o'1764 —o'1828 —0°1845 —0o'1886 
16 —0'1259 —0°1236 —~O'1291 0°1303 —0°1340 
he f —0°0874 —o'0856 —0°0g00 —O'ogiI —0'0940 
1°8 —o'o601 —0o0'0586 —o'0621 0'0629 —0°'0653 
1'9 00409 00398 0°0424 00430 00450 
2'0 —0'0276 —0°0269 —0'0288 -0°0293 0°0307 
2°25 O’oIo! 00098 O’O107 o-o0108 O’Oll4 
2°50 —0°0035 —0°0034 0°0038 —0°0038 0°0041 
2°75 —0°Oo12 —0°OO12 —0'0014 00014 00015 
3°00 —0©°0005 ~-0°0004 —©°0005 00005 00005 
3°25 —0*0002 —0°0002 0'0002 00002 —0°0002 
3°50 — 00000 —0*0000 —0*0000 —0©"0000 —0*0000 
Esp, 3p 28°552 28-912 28°176 28:076 27°494 

a(3p) 243°09 243°O7 242°74 243°16 242°22 


The writer is grateful to Professor D. R. Hartree for helpful discussions and to 
Dr M. V. Wilkes for permission to use the facilities available at the University 
Mathematical Laboratory, Cambridge. 
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ABSOLUTE OSCILLATOR STRENGTH MEASUREMENTS IN 
Mg, Ca AND OTHER ATOMS 


C. W. Allen 
(Received 1957 July 27) 


Summary 

The diluted copper alloy method of determining absolute oscillator 
strengths has been extended to Li, Na, K, Mg and Ca. Results could be 
obtained for Mg and Ca only. The original absolute intensity scale was 
confirmed. The results show fair agreement between observations and calcu- 
lations but some disagreements have been noted. A re-examination of earlier 
measurements show that the oscillator strengths of the heavy atoms Ga, 
Ag, Sn, Pb, and Bi are systematically greater than calculated values. 





The purpose of this work was to extend the earlier method of determining 
oscillator strengths (1) to some light atoms whose absolute values may be 
considered known already. This would check the generality of the method, 
re-evaluate the absolute scale and perhaps provide useful measures for the new 
atoms. With this in mind some copper rods containing known amounts of 
Li, Na, K, Mg, Ca and Fe were obtained from Johnson, Matthey and Co. Most 
of these atoms have lines with well known absolute oscillator strengths. Fe was 
included in order to check back with the original work. 

1. Method.—The experimental procedure was precisely the same as in the 
earlier work (1). Spectra of the diluted copper alloys were photographed and 
examined with a recording microphotometer. Photometric calibration was 
provided by a step filter, and colour sensitivity corrections were derived-using the 
positive pole of a carbon arc and also a tungsten strip lamp. The excitation 
temperature (4300 °K) and the absolute scale were taken from the earlier work, 
and possible corrections to them are discussed. 

The new alloys were supplied in two grades of which the “ strong” grade 
contained roughly 0-o1 per cent of each element and the “‘ weak ” grade nearer 
o-oo1 percent. Since the weak grade was obtained by dilution of the strong with 
pure copper one might have expected the concentration ratios (strong/weak) 
to be constant from one element to another. However the ratios, as supplied by 
the makers from chemical analysis, were Na 5, K 3, Li 11, Ca 135, Mg 33, and 
Fe 4. It is quite clear from the intensity measurements also that the ratios were 
not constant and that it was better to use the makers’ values than to adopt one 
ratio for all atoms. These figures show that the metalurgical problem of making 
the new rods was more difficult than for the previous ones (mostly elements of the 
iron group) and some errors resulting from the chemical analysis figures may be 
expected. It is thought that this is a significant source of error for Mg and Ca. 

Since the ionization potential was low for some of the atoms it was necessary 
to find whether ionization was great enough to affect the population of neutral 
atoms in the arc plasma. It is fortunately possible to estimate the degree of 
ionization by comparing the intensities of the neutral and ionized Calines. From 
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a number of spectra the relative intensities averaged (Ca I, 4226)/(Ca 11, 3933)=1°6, 
leading to 0-32 as the ratio of neutral to total calcium, With a temperature of 
4300 “K the Saha equation gave an electron pressure of 54 dynes cm~*. Using 
these values in the Saha equation again it has been possible to estimate the degree 
of ionization for the other elements. Li, Na, and K were mainly ionized and 
Mg, Fe, and Cu mainly neutral. Corrections for ionization were small for all the 
elements in (1). 

In (1) we assumed that the atomic composition of the plasma was the same as 
that of the pole and we checked this by studying spectra in which the arc was 
focussed onto the slit. The similarity of intensity distribution along the slit for 
all the elements gave us confidence in the assumption. When focussed spectra 
were obtained with the new poles it was found that the emission for Li and Na was 
strongly concentrated near the negative pole. The behaviour in the arc of Li 
and Na was evidently quite different from the other atoms and it would not be 
safe to assume they had the normal excitation temperature or diffusion conditions. 
In addition their ionization correction is large and therefore doubtful. It has 
been found in fact that the oscillator strengths derived from measurements in Li 
and Na were very low and inconsistent and these results have been discarded. 
No K lines were seen. The lines of Mg and Ca appear to be normal and should 
give usable oscillator strengths. 

The intention of the observations has been to study absolute values, and re- 
lative measurements have been somewhat sacrificed to this end. Because of the 
low concentration of alloy atoms the lines measured are often so faint that 
blending ambiguities become common and appreciable. It often happens also 
that the intensities of some lines are derived from one grade of alloy concentration 
and other lines from the other grade. ‘This has a rather adverse effect on relative 
accuracy. 

2. Results for Mg and Ca.—It is intended that the results shall (i) check the 
original absolute scale adopted in this work, (ii) compare present results with 
both observations and calculations from other sources, and (iii) derive “‘ best 
values” for the absolute oscillator strengths so that these may be available for 
astrophysical purposes. We follow recent practice in using the product of the 
weight w and the oscillator strength f to represent the resulting value for each 
spectrum line. 

In order to reduce the amount of tabulation the spectrum lines are not 
quoted separately; instead we tabulate log Lwf where Laf is the sum of wf for 
all lines of the multiplet whether measured or not. Since LS intensity rules are 
followed within nearly all the multiplets concerned it is easy to convert measure- 
ments to Lf and to derive individual line values from the tables. 

wf measurements in Mg and Ca were first derived on the original absolute 
scale previously used for the results in Fe, Cr, Mn, Co, Ni and other atoms (r). 
By comparing these with calculated values it was found that the new Mg measure- 
ments were about 0-40 dex high and the new Ca about 0:23 dex low*. These 
differences are no greater than might have been expected for absolute values by 
this method. As Mg is high and Ca low they give no indication of any significant 
error in the original scale; indeed they confirm the original scale to an accuracy 


about +0:2 dex. 


* Since results are expressed logarithmically throughout it is convenient to use the term “ dex ”’ 
(2, 3) instead of ‘‘ in the ten based logarithm ’’. 
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Although we are able to establish order-of-magnitude agreement between the 
observed and calculated values in Mg and Ca it is thought that the calculations 
give a better indication of absolute values than observations, and therefore 
in Tables I and II the observations have been adjusted to the calculated absolute 
scale. ‘The adjustments used are shown at the top of the columns. 

The columns headed “‘ Sun ”’ are derived from Fraunhofer line intensities. 
They are taken directly from the N.f tabulations in (4) and are not corrected 
for modern views on absorption coefficient and excitation temperature. In 
(4) the absorption coefficient adopted was proportional to A*, and the excitation 
temperature was 5600 °K. Absolute values are adjusted to agree with calcula- 
tions. 


Taste I. Values of log Xmf for Mg multiplets 











Observations Calculations 
Terms A CWA _ Ker Sun Trefftz B.D. B 

est 

(6) (4) (8) (9) (5) poor 

adjusted 
—0'°40 —O'31 —7°95 

7a ~Fe" wn —~sot —$°53 5°7 
—3'P® 2852 >—O'7 +O'21 +-0°22 +0°22 +o°22 

3°P° —3°D 3838 +075 +0°72 +o'79 +-0°77 +0°64 -or7I 
4°D = 3096 +O'25 +0°03 o'll 

—§°D 2851 —0°42 -0°36 -0°42 

6D 2736 0°79 o°5? o'8 
4S 5183 0°07 +0°17 +-0°32 0706 +008 —o0'18 0°02 

—5§°S 3336 —o'58 —o'78 112 —0°83 

—6°S 2942 —1°'20 —1°60 —1'4 
—p*?*P 2779 +078 +083 +0°78 -+-0'60 -O°72 
3'P® —3'D 8806 +005 —O'02 +0°34 0°04 
—4'D 5528 —0o'58 —0°54 —o'10 —1°OQ —1°32 —2°3? —o'75 
—5'D 4703 —0o'75 —0°43 +0°16 —1'0§ —1'20 —1°6? -0°73 
—6'D 4351 —1°Ol —0O'50 —1°'09 —1°8? -1'O 

~4'S 11828 —0°32 —0°32 -0°3 

~s'S 5711 <—ky. ~295 — 1°54 “$s 
—6'S 4730 <—-18 —2°71 2°08 —2°1 




















Calculated data are included where available, and Bates~Damgaard tabulations 
(5) have been used to obtain further values (labelled B.D.). When the application 
of the tables is doubtful a “?” is added. The B.D. calculations in Mg and Ca 
give fairly consistent results and a certain amount of extrapolation can be used. 

In order to obtain the Tables I and II “ Best value ” columns both laboratory 
observations and calculations have been considered, but not the Sun values. 
Before taking a weighted mean some minor adjustments have been made to the 
earlier columns. The Kersten values of Mg were further reduced (0-09 dex) 
to take B.D. calculations into account and the excitation temperature for the 
CWA measurements was increased from 4300 °K to 4700 °K without changing 
the mean Zwf. Reasons for this are discussed in Section 3. 
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Magnesium.—Kesults for Mg are in the CWA column of ‘lable I. The 
measurements were made almost entirely in the “ strong” grade of alloy since 
in the ‘‘ weak ” grade the lines were barely visible and badly blended with ghosts 
and band lines. The one line (A 5172) that was free in the weak grade was much 
weaker (about 0-8 dex) than expected on the basis of the strong grade spectrum. 
This left an ambiguity in the absolute evaluation. In an endeavour to give values 
on the original scale the strong grade values were reduced by 0-1 dex. We have 
already seen that they require a further reduction of 0-40 dex to agree with 
calculated absolute values. 


Tasie II. Values of log =awf for Ca multiplets 











Observations Calculations 
Terms A CWA Sch Sun B.D. Other Best 
value 
(10) (4) (5) (12, 5, 
adjusted 13, 14) 
70°23 FOr75 “Saf 
455 —4*P° 6572 ~ 3°65 - 4°07 4°16 | —471 
+o'16 
—4'P° 4226 +-O°05 —0°07 0°26 { + 0°37 +0°'22 
3 3g ' : L 0°06 ; 
4°P° 5°S 6162 -0'22 ~0'06 0°40 O15 O'12 
69S 3973 0°72 0'99 1°06 0°98 
7S 3487 1°52 1°54 —~I33 
-4D 4454 0°55 0°31 0°56 0°48 0°46 
-53D 3644 -0°33 o'21 oll o'21 
—6°D 3361 0°44 -0°44 ~0'46 
—p?*P 4302 0-70 0°46 +058 0°65 
37D 4°F° 4585 O'17 0°31 -0'5? o'2 
5°F° 40098 0°59 0°77 waite o'7 
dp*F® 6439 0°44 + 0°65 0°76 +0°64 
dp*D® 5588 0°66 +0"70 + 0°67 +0°68 
—dp*P® 5270 0°52 + 0°26 +0°46 Lo50 
dp'D® 6450 0°72 0°40 —o'6 
3'D 5'P° 5041 -0°76 0°70 ~2°6? —o'8 
6'P® 4526 0°89 0°78 — —o'9 
4'F° 4878 0°34 0°38 — —0°4 
dp'F® 5349 0°38 0°52 0°36 03 
4'P? —s5'D- 5188 0°31 O'24 0°46 —0'4 
p?'D 5857 O'14 0°22 0°30 +-o'22 
-p?'S 5867 1°80 —0°39 
eS  §sr2 0°59 1°39 




















For Mg there are several independent estimates of absolute wf. Kersten’s 
measurements (6, 7) were originally evaluated by fitting them to known value of 
the intercombination line A 4571. ‘Trefftz has made two calculations, (8) without, 
and (g) with, allowance for the influence of polarization. Of these the calcula- 
tions without polarization are in slightly better accord with experiment but the 
difference is not significant. Trefftz showed that Kersten’s measurements agreed 
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better with calculations if reduced by 0-31 dex and this adjustment has been made 
in the “ Ker” column of Table I. A third set of absolute values comes from cal- 
culations based on the Bates-Damgaard (B.D.) tables (5). For the multiplets 
that have been measured these are systematically 0-17 dex less than the Trefftz 
calculations but this is evidently an effect of selection because other multiplets 
(some not tabulated) restore the balance. 

The differences between the ‘‘ Sun ”’ values and the “‘ Best values ”’ are rather 
large but this is doubtless associated with damping constant variations. 

Calcium.—Results for Ca are in the CWA column of Table II. The reported 
ratio of the strong weak grade concentrations was 2:13 dex compared with a 
measured difference of 1-83 dex. For deriving a mean on the original scale 
double weight was given to the strong grade (i.e. corrections were + 0-10 dex for 
the strong and —o-20 dex for the weak). A further adjustment of +0-23 dex 
has been applied to bring into mean agreement with calculations. 

The measurements of Schuttevaer (10, 11) had been put on an absolute wale 
by comparison with the intercombination line A 6572 for which the transition 
probability was approximately known. However the comparison of these 
observations with B.D. calculations showed—with very good consistency—that 
a further correction of + 0-75 dex was required. This adjustment has been made 
in the “ Sch ” column of Table II. 

The “ Sun ” values show a good agreement with both observations and cal- 
culations. ‘The last two lines of the table have not been measured in the 
laboratory, but are added for comparison with the calculations. 

Except for the resonance and intercombination lines (5, 13, 14, 12) the only 
calculations available are by the B.D. tables. 

3. Discussion of Mg and Ca results —The Mg and Ca results give no indica- 
tion that there should be any substantial change in the original absolute scale. 

Fair agreement has been found between the independent measurements of 
Kersten (6, 7), Schuttevaer (10, 11) and myself. Comparison with calculation 
suggests that the observations of Kersten and Schuttevaer are somewhat more 
consistent than my own. 

There is reasonable agreement between measurements and calculations in 
both Mg and Ca. The irregularity of the intensity decline in the Mg series 
3'P°—n'D shows up in both the observations and the calculations, but the Trefftz 
calculated values for AA 5528 and 4703 appear to be too low. The B.D. calcu- 
lations for this series are unsuitable. In Ca the lines AA 5349, 5867, and 5512 
are probably examples of a significant difference between true and B.D.-calcu- 
lated f-values. For AA 5349 and 5867 the discrepancy appears to be associated 
with failure of the array multiplet intensity rules rather than the B.D. Coulomb 
approximation. 

In Fig. 1 the logarithmic differences between the calculated and measured 
values are plotted against E,, the upper excitation potential. Doubtful B.D. 
results are omitted from this diagram. ‘There is not much evidence for a large 
excitation potential effect such as found for the atoms of the iron group (1 5) and 
the slope in the CWA results is thought to indicate some error in the adopted 
arc temperature 7. There is no appreciable slope when Kersten or Schuttevaer 
observations are used. 

The sloping lines drawn on Fig. 1 are considered to represent the effect of 
error in T, and the slope may be corrected by increasing T from 4300 °K to 





SER hi PES MERA TAN SID i Nae “itis Be aly Po 


2 se atihg hens io 


Lact slr Te lst calatair ene ines Bt eat 











less eS Phatlste aaatrentaies ine 





hk S086 A RA ct St 





PW cl mele ete Babli oe 








No. 6, 1957 strength measurements in Mg, Ca and other atoms 627 


4700 °K. As explained in (1) the original temperature was obtained almost 
entirely by comparison with furnace observations by King and his collaborators. 
The higher temperature now suggested would reduce the excitation potential 
effect that still persists in Fig. 4 of (x). It would also increase the derived solar 
reversing layer temperature from 4850 °K to 5370 °K which is nearer the probable 
solar temperature at the depth (r=0-3) where Fraunhofer lines are produced 
(data from 3, pp. 133-4). It is thought that these points justify the higher tem- 
perature and this has therefore been used to amend the CWA results for deter- 
mining the ‘‘ Best value” in Tables I and II. We will find that the change intro- 
duces some improvement in Section 4. 
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Fic. 1.—Relation between calculated and measured oscillator strengths plotted against excitation potential. 


4. Applications to other atoms.—The results given in (1) lead to the conclusion 
that the probable error of the present method when used for transferring absolute 
oscillator strengths among elements of the iron group is about +0-2 dex. 
Applying the method to Mg and Ca gives ambiguities suggesting a greater 
probable error, perhaps +0-3 dex, while for Li, Na, and K the method is not 
applicable. ‘The question now arises whether the published values (1) for Al, Si, 
Ga, Ag, Sn, Pb and Bi can be used with any confidence. 

The results for these atoms are collected in Table III where they are compared 
with calculations. ‘The CWA column is on the original scale of (x) but the spectra 
have been remeasured and some considerable alterations have been found 
necessary for Ga, Ag and Bi. B.D. calculations have been attempted but these 
are not very satisfactory for the heavy atoms because term energy levels are widely 
separated and LS coupling rules not obeyed. However these inaccuracies should 
not have a systematic effect on the calculations. Wave-function calculations 
have been made for Al only (16) and these do not agree well with B.D. 

Heavy atom oscillator strength measurements are consistently higher than the 
calculations and we must seek some explanation for this difference. There is no 
reason for thinking that the heavy atoms behave abnormally in the arc because, 
unlike Li and Na, the distribution of heavy atom line intensities throughout 
the arc plasma is much the same as for atoms of the iron group. One possibility 
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to consider is that the heavy atoms remain longer in the arc plasma because the 
diffusion velocity is smaller. A reasonable correction for diffusion speed might 
be obtained by introducing a concentration factor inversely proportional to the 
thermal velocity, i.e. proportional to the square root of the atomic weight. Such 
a correction improves the agreement. Some further improvement may be ob- 
tained by adopting T=4700°K as suggested in Section 3. Both of these 
corrections are introduced into the column labelled ‘‘ Corr” of Table III. The 
observations for heavy atoms are still systematically high but the improved 
agreement is an indiction that the corrections were necessary. Some inde- 
pendent experimental confirmation is desirable and for this purpose absolute 
measurements in Pb are contemplated. 


Tasie III. Values of log Xwf for other atoms 











log Loaf 
Atom | Configurations Terms A Observations Calculations 
CWA Corr. B.D. Wave 
Al 3p-4s *pe-2S 3961 “I'l —o'9 —0o'52 —0'09 
Si 3p*-3p45 sda il 3905 | —o'91 —o'94 | —0°87 
Ga 4p-5s 2pe-25 4172 +o0:08 +005 | —o'48 
Ag 5s-5p 2S —*p? 3382 +090 +0°73 +0°20 
Sn 5p*-5 pos 'D —'p° 3262 +102 +0°71 -0°6 
1S -'p? 4524 to'1g —o'16 -O'9 
Pb 6p*—6p7s 7? 3639 +141 +4103 | —0'%4 
1D -'p? 3572 +1°I3. +0°§7 —o'l 
Bi 6p*-6p*7s 2—°—2P 3397 +o'6 torr —o°2 



































Of the two light atoms studied, Mg measurements are high and Al low by 
comparison with calculations. Doubts have been expressed (1) about the 
reported concentration of Al and these results do not allow of any generalization. 

University of London Observatory, 

Mill Hill Park, 
London, N.W.7 : 
1957 July. 
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Summary 


In applications of the two-streams theory, it is tacitly assumed that the 
velocity dispersions of the two drifts are the same. This paper describes 
a method whereby, assuming the positions of the drift apices and the values 
of h,W, and h,W, found from an analysis of the proper motions of a group of 
stars, the observed radial velocities of the group can yield a determination 
of h, and h, explicitly, together with the proportions of stars in the two 
streams. ‘The method is applied to stars of spectral types Ko-M, using 
Tannahill’s analysis for the GC proper motions and the radial velocities of 
Wilson’s GCRV. ‘The parameters of best fit give (for the GCRV stars) 56°5 
per cent of stars in Drift I, and essentially equal values for h, and h,. Thus 
the inclusion of different dispersions for the two drifts as unknowns does not 
improve the agreement of the two-streams theory with observation. 





1. Introduction.—It is well known (1) that there are two general models for 
correlating the observed motions of stars. The ellipsoidal theory supposes that 
the stars show a greater mobility in one particular direction than in others, so 
that, after the removal of a component of the observed velocities attributed to 
motion of the Sun, the distribution-function of velocity vectors of the stars 
has an ellipsoidal form. ‘The other model, the two-streams theory, supposes 
that the stars can be divided into two intermingling groups, designated Drift I 
and Drift II. For each drift, the velocities of the stars are supposed to have a 
Maxwellian distribution when measured relative to some centre of motion. The 
complete specification of a two-streams model would then consist of a statement 
of the velocities of the centres of motion of the two drifts relative to the Sun, the 
proportions of stars in the two drifts, and the velocity dispersion in each drift. 

Primarily because it can be more easily related mathematically to the 
dynamical properties of the Galaxy, the ellipsoidal theory has been preferred. 
It can be shown that, in the special case of equal numbers of stars in each drift 
and identical drift dispersions, relations exist between the constants of thé two 
types of model (2). However, from the phenomenological viewpoint, the greater 
flexibility of the two-streams theory may be deemed an advantage if the 
additional disposable constants enable a better correlation with the observed 
motions to be obtained. 

In practical applications, it is usual to include in the drift analysis the ratio 
of the number of stars in the two drifts as an unknown, and frequently the best 
fit with observations is obtained with a non-equal partition of stars between the 
two drifts. However, little, if any, practical use has been made of the other 
degree of freedom inherent in the two-streams theory, namely the possibility 
of different velocity dispersions for the two drifts. ‘The reason for this is not 


‘hard to find—it is simply that, from proper motions alone, the only measures 
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that can be obtained for the velocities of the drifts relative to the Sun are of the 
form h,W, and h,W,, where W, and W, are the drift velocities, and h, and h, 
are moduli of the velocity dispersions. Thus, from proper motions alone, it is 
necessary to make some assumption about the relative values of the two velocity 
dispersions in order to d&ive such quantities as the solar motion and the positions 
of the vertices of star-streaming. ‘The simplest assumption, and one that is 
always made in practice, is that h,=h,. The use of radial velocities in drift 
analyses seems to have been largely confined to showing that the radial velocities 
do not conflict with the two-streams hypothesis. However, a combination of 
proper motion and radial velocity statistics can (in principle) lead to an 
independent determination of the proportion of stars in the two drifts and the 
values of each of the moduli h, and h,. It is the purpose of this paper to describe 
a method whereby such a determination can be made, and then apply it to a 
selected class of star. In particular, the paper seeks to discover whether, for the 
particular class of star chosen, the inclusion of two different moduli as unknowns 
improves the correlation of the two-streams model with the observed motions. 

2. Theory.—It is assumed that the particular group of stars under examination 
has been subjected to a conventional two-streams analysis of proper motions, 








APEX II APEX I 
(A,,0,) (A,,D,) 
Nae «58 
VERTEX 
(Ay, Do) 
SOLAR APEX 
Fic. 1. Fic. 2. 


yielding the coordinates of the apices of the two drifts (A,, D,) and (A,, D,), 
as in Figs. 1 and 2. Consider now a small area of sky containing N stars, the 
angular distances of the centre of the area from the two drift apices being A, 
and A,, as in Fig. 2. The motions of the stars within this area are to be 
represented by a combination of two drifts; suppose that a fraction « of the 
stars can be supposed to be members of Drift I. The radial component of the 
velocity of Drift I relative to the Sun will be W,cosd,. For completeness, it 
is necessary to take into account two small terms in the radial velocities, due 
to galactic rotation and the so-called ‘“‘ K-effect”. The first term has the form 
(rA) sin 2(G— Gy) cos*g, where r is the mean distance of the stars, A is Oort’s 
constant, G, g are the galactic longitude and latitude of the centre of the area 
respectively, and G, is the galactic longitude of the centre of the Galaxy. This 
term may be considered to be constant over a given small area. The K-effect 
is represented by a constant radial velocity K. Then a typical observed radial 
velocity p of a star within the chosen area, and a member of Drift I, will be 


p=p,+ W, cosa, + K+(rA) sin 2(G—G,) cos*g 
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where p, is the random radial velocity, or, writing K + (rA) sin 2(G — G,) cos? g = Z, 
By definition of a drift, the random radial velocities p, will have a Maxwellian 


distribution, so that the number of stars of Drift 1 with observed radial velocity 
between p and (p+dp) will be given by 


Nh, d 
n,(p) dp = — : x P exp[—h,(p — W, cosd, —Z)?], 





or, on writing x=/,(p— W, cosa, —Z), 
N 
n,(p) dp = —= e~" dx. 
V 7 


Similarly, for the stars of Drift II, on writing y=h,(p — W,cosA,—Z), 
_ —a)N ore 
n,{p)dp= “— 2 e-'dy. 
Hence the total number of stars within the area with observed radial velocity 
between p and (p+ dp) is given by 





N ; : 
n(p) dp= war dx+(1—a)e~” dy]. (1) 
The mean radial velocity for the area, p, is given by 
p=a(W, cosa, +Z)+(1—«)(W,cos ry + Z) 
or p=«W, cosa, +(1—«)W,cosdA, + K + (rA) sin 2(G — Gy) cos? g. (2) 


If the sky be divided into areas, for each of which the angles A,, A,, G and g 
are known (G, being assumed), and p found from the observed radial velocities 
within the area*, then for each area we have an equation of condition of the 
form (2) which can be solved by least squares to give values for «W,, (1—a«)W,, 
K and (rA). The results of such a solution for the radial velocities considered 
in this paper are given in Section 4, Table II. 

The proper motion analysis will have given values for «, h,W, and h,W,. 
Thus, in principle, it would be possible to solve uniquely for «, h,, ha, W, and W, 
without further assumption. However, the value of « is likely to be poorly 
determined from the proper motions. In addition, in the application subse- 
quently to be made, the proper motion analysis refers to stars in the GC, while 
the radial velocity data are obtained from Wilson’s General Catalogue of Radial 
Velocities (GCRV) (3). Even for a given range of spectral type, the selection 
of stars in the two catalogues will be rather different, and a comparison of various 
drift analyses suggests that different values of « might be appropriate for different 
selections. (‘This point will be discussed in Sections 4 and 6.) For these reasons, 
it was decided to treat only h,W, and h,W, as known from the proper motions, 
and to use the radial velocity data themselves to determine the values of «, h, 
and h,. The method of solution adopted will now be described. 

Let n be the number of stars within an area, containing a total of N stars, 
whose observed radial velocity p<p,, where p, is the mean radial velocity for the 


* In the actual analysis, stars with p numerically greater than 100 km/sec were omitted 
(see Section 3). 
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area as calculated from the solution of equations (2). Then, if x,, y, are the 
values of x, y corresponding to p,, we have, from equation (1)*, 


=f" Pine (1—a)N f" ev dy 





Oe Va va 
ON, 3 (“ee Gao ( 2 ("erg ) 
=F (14+ Z[e ds) + 2 It ¢ ry 
or 
2n 
Waite erf (x,) + (1 —«) erf (y,). (3) 
Now 
x,=h,(p,— W, cos A, — Z) 
=h,[(a—1)W,cosdA,+(1-—«)W,cosd,] from (2) 
= —h,(1-—a)A 
where 
A= W, cos A, — W.. cos Ag. 
Similarly, 
Ve =+ hyo. 
Since erf (— x)= —erf (x), equation (3) becomes 
(: i wy) = aert Ce ~a)A]—(1—a) exf[AyeA). (4) 


Now, from the solution of equations (2), we find «W,/{(1—«)W,} =a, say, 
while from the proper motion analysis we find h, W,/h,W, = 6, say. Thus we have 
a ahs 


a known constant. Writing p=h,(1—«)A, equation (4) becomes 


(: ~ x) =aerf(p)—(1—«) erf(Cp) 





or 
2n 


: (:- *) + erf (Cp) 


 erf(p) + erf(Cp) (5) 


Equation (5) holds for positive and negative values of A. 

An equation of the form of (5) will hold for each area whose angles A, and A, 
are known, the quantity (1—2mn/N) being an observed quantity. Equations (5) 
can be solved for « and A, in the following way. Making an initial assumption 
as to the value of a, say «=a», the values of W, and W, can be found from the 
least-squares solution of equations (2). Hence A can be calculated for each area, 
and h, can be found from the value of h,W, of the proper motion solution. The 
values of p and Cp can thus be found for each area, and hence each equation of 
the form of (5) yields a value for «. A mean value «, can then be found from 
all the equations (5), and this new value used to calculate new values for W,, W, 
and f,, and hence for p and Cp for each area. The process can be continued 
until a stable value of « is obtained. Finally, the values of W,, W,, h, and h, 
can be calculated using this stable value of «, and hence also the coordinates of 





* See footnote on previous page. 
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the vertex of star streaming, and the direction and magnitude of the solar motion. 
The results of the application of this method to the observed radial velocities 
used in this paper are given in Section 5, Table III. 

3. Choice of stars.—The proper motion analysis with which the present 
analysis of radial velocities is to be compared is that of Tannahill (4) for the 
GC stars. Clearly, it is desirable to choose a definite spectral type range, as the 
stream constants show some variation from one spectral type to another. The 
range chosen was Ko-—M, since this is the range with the greatest number of stars. 
both in ‘Tannahill’s analysis and the radial velocity catalogue. In extracting 
data from the GCRV, the following categories of star within the chosen range of 
spectral type were omitted: (i) components of double stars; (ii) stars whose 
radial velocity is indicated as variable ; and (iii) stars whose spectra show emission 
lines. ‘The last category was considered unreliable because the presence of 
emission lines may indicate the presence of a gaseous envelope to the star, so 
that any determined radial velocity may be distorted by motions within the 
envelope. Finally, during the analysis, it was found that a small fraction of the 
stars had very large radial velocities which had an undue influence upon the mean 
radial velocity of their areas. For this reason, all stars (a total of 37) with observed 
radial velocities numerically greater than 100 km/sec were also omitted. There 
remained 3713 stars for analysis, compared with 10666 for Tannahill’s proper 
motion analysis of the same spectral type range. 

The relevant proper motion data for the chosen spectral type range are 
tabulated by Tannahill (4). His value of 1-02 for the ratio of the number of stars 
in Drift I to the number in Drift II corresponds to «=0-505. 

4. Solution for aW,, (1—a«)W,.—To solve equations (2), it is necessary to 
assume the positions of the drift apices. ‘These can be taken from Tannahill’s 
analysis. However, as the GC includes more stars of large magnitude than the 
GCRV, some evidence is required that the positions of the apices do not show 
any systematic change with stellar brightness. Such evidence for spectral 
range Ko—M is provided by the tabulation of Ewart (5, Table 1), where the 
catalogues considered, in order of increasing upper magnitude limit, are 
(i) the GC, (ii) the Cape Catalogue 1936 (6) and (iii) the Cape Catalogue 1939 (7). 

Reference to this table also shows that the variation between one catalogue 
and another as regards the positions of the drift apices is larger than the probable 
errors given by Tannahill. In the subsequent calculations, some saving of labour 
is achieved by assuming that the two apices differ in right ascension by 180°. 
For this reason, the positions of the apices for the present work were assumed. 
to be 

A,= 94°, D,=—-11°, 
A,=274°, D,= —74°. 


These positions lie only a little outside the range of probable error given by 
Tannahill. With these assumed positions, the angle 6=180°+ D,+D,=95° 
(see Fig. 2). Graphs were prepared giving A,, A, and % against (A, D), the 
equatorial coordinates of a typical star (see Fig. 2). For each star, the angles 
were read to the nearest degree. The convention was adopted that varied 
from — 180° to + 180°. 

In collecting the stars into areas, it is necessary to come to a compromise. 
The larger the area, the greater the variation of A, and A, across it. On the other 
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hand, the larger the area the greater the number of stars within the area, an 
advantage when the crude examination of the radial velocity statistics within 
an area is made by means of equation (5). It was considered that, for the 
application of equation (5), a minimum of about 30 stars was required, and the 
sizes of the chosen areas were determined by this consideration. As can be 
seen from Table I, only one area (No. 21) falls short of this minimum requirement. 

In forming the mean radial velocity p for a given area, some consideration 
must be given to the individual weights to be applied to each determined radial 


TaBLe I 


Observational data 


Area Range Range yp N Mean Mean Pp n 
No. of A, of A, A, As km/sec 
I O- 15 All = 77 10°6 96°6 23°01 37 
2 15- 35 <110 70 26°1 go'o 18-64 37 
3 15- 35 <110 63 26°3 89°8 ~18°51 31 
4 15- 35 110-131 73 25°8 117°7 “10°24 39 
5 I5— 35 110-131 — 50 28°3 120°2 10°50 =. 29 
6 35- 55 <100 67 46°3 75°5 15°50 37 
7 35- 55 <100 42 46°8 71°O + 24°88 17 
8 35- 55 100-125 70 451 114°0 + 47°12 31 
9 35- 55 100-125 ~ 73 43°2 116-2 +1703 2 
10 35- 55 125-146 + 60 46°6 137°! 10°50 25 
II 35- 55 125-146 - 102 46°1 135'1 6°25 53 
*12 55- 75 < go 63 65:2 59°1 g18 38 
*13 55- 75 < go -- 44 65°0 65°3 12°43 18 
14 55- 75 go-120 66 64°6 1051 + 7°86 35 
15 3S” 75 go-120 — 39 64°9 107°I ~ 7°94 17 
16 55- 75 120-145 r 54 65:2 133°3 + 1°42 26 
17 55- 75 120-145 — 64 65°6 133°7 + 2°56 28 
18 5S- 75 >145 58 65°5 155°6 086 24 
19 53-75 145 — 63 64°9 154°4 0:96 28 
*20 75- 90 < 80 43 82-6 48-1 6:22 29 
*21 75- 90 < 80 — 26 82:2 49°3 6°28 15 
*22 75- 9° 80-120 + 51 82-7 100°7 6-92 19 
23 75- 90 80-120 - 61 82-4 106°1 3°61 32 
24 75- 9° 120-145 t 33 82°7 133°6 2°01 16 
25 75- 90 120-145 — 47 83°4 132°0 7°20 23 
26 75- 90 > 145 : 59 83-0 161°6 1°05 2 
27 75- 90 >145 — 78 81°6 158°8 — 10°63 46 
28 QO-105 < 80 46 98°1 59°8 -~ 1°62 2 
* 90-105 < 80 -- 59 97°5 55°9 + 2°62 29 
a 90-105 80-110 57 98-9 96°5 0°39 2 
*31 QO-105 80-110 = 47 97°3 96°9 5°44 9.23 
*32 QO-105 110-140 } 42 96°5 122'0 sree 16 
33 QO-T105 110-140 = 70 97°8 122°6 -11°77 37 
34 QO-105 140-155 - 40 99°9 149°2 wd 17 
35. QO-105 140-155 ~ 56 97°7 148°0 16°73, 32 
36 90-305 “155 + 40 95°6 163°5 —11°73. 20 
37 QO-105 >155 ~— 42 97°3 161°8 — 14°36 21 
38 105-120 < 80 38 112°6 50°3 — 2°91 19 
39 105-120 < 80 _ 33 112°6 53°1 5°46 10 
40 105-120 80-110 + 37 112°0 94°0 — 4°42 14 
41 105-120 80-110 — 51 113°5 96°5 — 12°08 26 


*42 105-120 110-140 + 52 i11°7 125°7 —16°30 = 30 
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TaBLeE I (continued) 
Observational data 


Area Range Range ys N Mean Mean Pp n 
No. of A, of Az A Az km/sec 

*43 105-120 110-140 - gl 112°3 127°5 —12°04 42 
*44 105-120 >140 42 113°! 142°3 —19°41 27 
*45 105-120 >140 7° 111°8 147°0 16°42 33 
46 120-135 < 80 t 65 128-4 61°0 — 9°43 30 
47 120-135 < 80 5° 127°0 59°0 — $862 25 
48 120-135 80-110 t 60 128-4 981 12°61 29 
49 120-135 80-110 - 65 130°2 96°2 —16°41 35 
50 120-135 110-130 36 125°3 77 19°41 17 
51 120-135 110-130 71 127°3 120°2 —17°90 38 
"xo 120-135 >130 40 125°3 136°6 — 20°92 16 
*53 120-135 >130 64 125°3 136°3 12°98 = 27 
54 135-150 < 80 36 142°4 57°5 ~11'94 = 20 
55 135-150 < 80 ~ 52 143°6 58:3 —- 890 29 
56 135-150 80-110 48 144°1 93°9 — 53°73 24 
57 135-150 80-110 54 143°4 97°9 — 13°95 24 
58 135-150 “110 : 70 142°4 117°5 20°83 36 
59 135-150 >110 84 141°9 117°5 -17°68 39 
60 150-170 < 80 32 158°6 69°4 — 14°38 19 
61 150-170 < 80 - 58 159°3 69°3 — 19°67 34 
62 150-170 80-100 61 161-0 go'l 14°58 34 
63 150-170 80-100 93 161°3 go'l -13°06 4! 
64 150-170 > 100 55 156°3 105°7 20°67 22 
65 150-170 > 100 _ 54 1551 104°7 18-98 28 
66 170-180 All + 56 172°3 85°1 —18'49 29 


velocity. From the discussion in the introduction to the GCRV, it is clear 
that little precision can be attached to the relative quality-criteria, and it was 
decided to give each star unit weight. An exception might have been made for 
determinations labelled ‘‘e”’, as these are clearly much poorer than the other 
determinations, but as the proportion of quality “‘e”’ determinations is so small, 
it was considered not worth while making even this distinction between one 
star and another. For the least-squares solution of equations (2) each area was 
given a weight approximately proportional to the number of stars included 
in the area. The values of A, and A, appropriate to each area were taken to be 
the means of the corresponding angles for the individual stars, each star being 
given unit weight. The sign of y% was used to distinguish between any two areas 
with similar values of A, and A,._ ‘The values of G and g were read from the Lund 
Tables (8), G, being taken to be 325°. ‘The result of the least-squares solution 
of equations (2) is given in Table II. 


TaBLe II 
aW, 20'2 +0°5 km/sec 
(1—a)W,= + 8-7+0°7 km/sec 
K=-+ 0°3+0°4 km/sec 


(rA)=-+ 1:9+0°6 km/sec 
A K-term which is small and positive is in good agreement with the accepted 
values of o-o for type Ko and +04 for type Mo(g). Since giant and dwarf stars 
have not been separated in the analysis, no estimate of r can be made to compare 
the value of A with accepted values, but with the value of A quoted by Allen (9), 
a distance r~100 pcs is found. ‘This would seem to be a very reasonable figure, 
and to give confidence in the solution. 
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If the value of « from Tannahill’s analysis is supposed also to apply to the 
GCRV stars, then a determination of W,, W,., h, and h, is possible. In Fig. 3 
is plotted the value of the ratio h,/h, against «. The shaded area shows the range 
of this relationship allowing for variations in «W,, (1—«)W2, h,W, and h,W, 
within the ranges given by their probable errors. The single cross represents 
the position of the point corresponding to Tannahill’s value of «, with h,=h’,. 
This figure shows that the data have to be stretched in order to give a consistent 
solution with equal dispersions for the two drifts. 

The choice thus lies between the following assumptions: (i) that the value 
of « applies to both the GC and the GCRYV stars, giving a dispersion ratio 
h,/hy~o-79; (ii) adopting « from Tannahill’s analysis of the GC, and supposing 
his values of h,W, and h,W, are not applicable to the GCRV stars; or (iii) taking 
h,W, and h,W, from the GC as applicable to the GCRV stars, and determining a, 
h, and h, from equations (5). There is little doubt that the third assumption is 
the most realistic one, as it is well known that the proportions of stars in the two 
drifts vary with magnitude range when all spectral types are included (compare, 
for example, the results of (10) and (4)). Ewart’s Table II (5) also suggests 
that for individual spectral type ranges, « varies with magnitude range. In the 
next section, assumption (iii) is thus explicitly made. 


a 
1.0 
0.879 
~ 2 “YS 
06— SRRSSw—ws 


0.4— 
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The residuals p —p, (where p, is the mean radial velocity of an area calculated 
from the results of Table II) were plotted against galactic longitude, but no 
obvious correlation was noted. A similar result holds for the plot against galactic 
latitude. 

5. Solution for « and h,.—The method of solving equations (5) by successive 
approximations has been described in Section 2. The only additional point 
to be mentioned here concerns the relative weighting of the areas. From 
equation (5) it is seen that, for a given area, a change 8g in g-=(1—2n/N) 
produces a change 5q/{erf(p) + erf(Cp)} in «. Thus each area should be weighted 
by a factor proportional to | erf(p) + erf (Cp) |= ¢, the value of 5g being supposed 
to be independent of «. At first sight, this argument might appear invalid 
because p and Cp for a given area depend upon. But this is not so. The purpose 
of the analysis is to fit the observed radial velocities as closely as possible to a 
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two-streams model. Ideally, for every area the value of « would be the same, 
and W, and W, determined uniquely from the values of «W, and (1—a)W, 
given in Table II. In practice, we have to allow for a variation of « from area 
to area. However, we do not consider the possibility of W, and W, varying from 
area to area. The problem is thus to find the value to choose for «, say «=a,, 
to give W, and W, from Table II, such that «, is the best possible mean value 
of the a’s for the individual areas. Under these circumstances, once a value 
of a, has been assumed to define the drift velocities, the quantity ¢ is determined 
for each area, and cannot be considered to vary with the local value of « for the 
area. ‘The process of successive approximations ensures that the best mean 
value of the a’s is adopted for «,. The final weighting of an area was thus taken 
to be proportional to N¢. 

During an iteration, it was found that certain areas gave impossible values 
of « (i.e. negative values or values greater than unity). Such areas were always 
found to be those for which A was numerically small. Rather than arbitrarily 
ignoring impossible values, it was decided to ignore, in the calculation of the 
mean value of « for a given iteration, those areas with |A|<1o0 km/sec. (‘This 
figure should be compared with the supposed average probable error of between 
1 and 2km,sec for an individual radial velocity.) ‘The areas so ignored in the 
final iteration are indicated by an asterisk in ‘Table I. (It should be pointed out 
that, owing to the changes of the assumed W, and W, from one iteration to the 
next, the actual areas ignored are not quite the same from one iteration to the 
next.) 

The results of the successive iterations are :— 

%=0'500 (assumed) 

%=0'557 (0°543) 

%=0°565 (0°559) 

%=0'565 (0°559). 
The figures in brackets give the means found if the areas are unweighted by the 
factor £. The fully-weighted means are to be preferred. ‘The third iteration 
was therefore considered to give the best representation of the observations, 
and the full numerical results of this iteration are given in Table III. This table 
also includes the position of the vertex of star-streaming in equatorial coordinates 





Tas_e III 
Solutions 
ny a 
a=0'565 (+0°083) = ea =1°30 
W,=35°8(1) km/sec h, =0°0387 
W2= 20°0(0) km/sec hy =0°0385 
W,=42°4 km/sec 0. =hW,= 1°64 (»- me) 
Solar motion 
U,=21°3 ~km/sec Ayp= 274 
Do= +35 


Vertex of star streaming 
l,= 342° a,= 274° 
b= —2°8 8,=—17 
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(a,, 5,) and galactic coordinates (/,, 6,), the relative drift velocity W3, the solar 
motion U, and the equatorial coordinates of the solar apex (Ay, Dp), all computed 
in the usual way. The values of «, and Ay necessarily arise from the assumption 
that the drift apices differ by 180° in right ascension. 

The deviations («—a,) for individual areas from the final iteration were 
represented on the galactic longitude-latitude plane, to see if any systematic 
trends could be noticed. Two effects of possible significance were noted, namely 
(i) a tendency for positive values of («—a 3) to be associated with galactic 
longitudes in the range go° to 200° (approximately), i.e. a range centred on the 
longitude of the galactic anti-centre at /= 145°, and (ii) a tendency for negative 
values of (a—a,) to be associated with galactic latitudes within the range 
+20° to —20° (approximately). Neither of these effects is very pronounced, 
but, if they are interpreted as real effects, they suggest that the proportion of 
stars in Drift II is greatest in regions of the sky where the star density is greatest. 
This result bears some similarity to the conclusions of Hough and Halm (11). 

6. Discussion—A comment should be made about the significance of the 
rather large probable error for « in Table III. In normal statistical work, it is 
assumed that the observational data represent a sample of a potentially very 
large number of observations, whose statistical properties are sought. The 
probable errors then represent the precision with which the properties of the full 
set have been determined from the sample. In this work, however, we are not 
taking a sample from a much larger number of possible observations; we are 
dealing with effectively all the observations that can be available (within a given 
magnitude range) and are trying to find the parameters of best fit of the adopted 
model to the finite number of observations. A probable error, in this context, 
measures the degree to which the adopted model fits the data, and not the 
precision with which the parameters of best fit have been determined. In the 
present instance, no doubt the probable error of the calculated « could have 
been reduced by taking somewhat larger areas, but the additional labour would 
not be worth while, for the result would not be rendered physically more 
significant. 

On the other hand, Table III does answer, quite unambiguously, the question 
posed at the end of Section 1. For the particular class of star examined here, 
the inclusion of two different dispersions for the two drifts does not confer any 
additional advantage on the two-streams theory, since the parameters of best 
fit are such that the two dispersions are effectively equal. Indeed, the agreement 
between h, and h, in Table III is very much better than could have been expected. 
Since the dispersions are effectively equal, the agreement between the derived 
quantities in ‘Table II] and the corresponding quantities in Tannahill’s paper 
1s a necessary consequence of having assumed drift apices as given by Tannahill. 

‘Table III also shows that the best fit of the model with observation is given 
by a value of « which is greater than 0-5. The disagreement between the value 
of « found here and that given by Tannahill is no doubt a consequence of the 
fact that the selection of stars in the GCRV is different from that in the GC. 
Cc onsidering all spectral types, there is a general tendency for « to decrease towards 
0°5 as fainter stars are included (10, 4). Such a systematic variation is not 


obvious from the limited data available for discrete spectral type ranges (5), 


but the value of « obtained here is well within the range of determined values 
from other investigations. 
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RADIAL VELOCITY OBSERVATIONS 
OF SOUTHERN VISUAL BINARIES 


S. Archer 


(Communicated by the Radcliffe Observer) 


(Received 1957 October 14) 


Summary 


Radial velocities have been determined for the components of the following 
five systems : 


ADS 6914, « Cen, ADS 10417, Mlb 4 and y Cr A. 
Only for the first two pairs a significant velocity difference, due to orbital 
motion, was found. The known fact that ADS 10417 C (=HD 155849) is 


optical to the bright pair was confirmed by the radial velocity. The known 
membership of HD 156026 was confirmed by the radial velocity. 





Introduction and observations.—It is well known that observations of radial 
velocities of the components of visual binaries may be useful to fix the sign of 
the inclination, for a rough determination of the mass-ratio and in favourable 
circumstances may lead to a valuable result for the parallax (1, 2). 

The stars treated in the present article were selected from a list kindly supplied 
by Dr W. H. van den Bos. All but two of the spectra were obtained by the writer 
with the 2-prism Cassegrain spectrograph attached to the 74-inch Radcliffe 
reflector. A camera giving a dispersion of 29 A/mm at Hy has been used 
throughout. The spectra of each of the components of a pair were photographed 
on the same plate and in quick succession so that possible systematic errors 
should affect the results for both stars equally and cancel in the difference of 
radial velocity. At least four pairs of plates were obtained per system. 

In all instances of good seeing it was possible to observe the components 
separately—the distances ranging from 2” to 4”, the slit-width being 07-3. The 
orientation of the slit was east-west for all systems except Mlb 4 for which the 
spectrograph was rotated so that the slit pointed north-south. - 

The observations are given in Table I. 

Table II gives the mean results for each star together with the mean difference 
between the components at the mean epoch of observation. 

The mean difference from five spectra with five standard velocities was 
found to be +1-1 km/sec. This difference was, however, not applied in the 
results given in this paper. 

ADS 6914 (8208): The visual orbit of this star was subject to a discussion 
by W. P. Hirst (3) and by Dr W. H. van den Bos (4). They concluded that 
no reliable orbit can be found at present and that more observations are needed 
to settle the question of a possible third body. ‘Though radial velocity observations 
are difficult because of the small separation it is hoped that further plates of 
this system will be obtained during the next few years. The 1957 plates were 


kindly made by Dr A. D. Thackeray and Dr M. W. Feast after the writer’s. 


departure from Pretoria. 
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Star 


Dec. 1950 
ADS 6914 
(North) 
gh 37™-o 
—22° 30’ 


ADS 6914 
(South) 

gh 37™-0 

—22° 30° 


aCenA 


145 360-6 
—6o0° 38’ 


a Cen B 
142 36-6 
—60° 38’ 


ADS 10417 
(North) 


— 26° 32° 


ADS 10417 
(South) 


17h 12.3 
—26° 32’ 


ADS 10417C 

=HD 155849 

172 y2M-1 
—26° 29’ 


ADS 10417 

HD 156026 

172 173-1 
—26° 29’ 


Mlb 4 
(brighter) 
17 15-4 
—34° 57’ 
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my 


5°4 
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o'3 


5°3 


5°3 


aT 


6°7 


U.T. 


1956 V 
Vv 


1957 | 
Vv 


1956 V 
V 


1957 I 
Vv 


1956 V 
Vv 
Vv 
VI 
VI 


1956 IV 


TABLE I 


14°693 
29°681 

4°972 
12°692 


14°705 
29°691 

4°988 
12°706 


14°905 
29°796 
29°852 
17°79! 
22°766 


30°871 
14°888 
29°807 
29°840 
17°804 
22°777 


1°094 
14°981 
29°994 
17°878 
22°846 
28-861 


17089 
14°974 
29°989 
17°884 
22°854 
28°854 


14°994 
1°880 


30°026 
17°9OI 
22°875 
28-876 


15°029 
29°065 
17°976 


R.V. 
km/sec 
+37°4 
+37°0 
+38°7 
+45°1 


+43°0 
+45°8 
+45°1 
+46°5 


—20°2 
—21°4 
—21°! 
—23°3 
—21°3 


+ o°7 
+ o6 
+ 08 
— 2°97 


+ 4°6 
+. 29 
— 02 
+ 1°6 





Remarks 
\ Seeing poor at times. 


Seeing poor at times. 


Optical companion to 


ADS 10417. 


700” distant from 
ADS 10417 c.p.m. and 
c.r.v. physical com- 


panion. 
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Star 
R.A. 
Dec. 1950 


Mlb 4 
(fainter) 
170 ysM-4 
—34° 57’ 


yCrA 
(North) 
19 903-1 
— 37° 08’ 


yCrA 
(South) 
198 03™-1 
— 37° 08’ 


My 


eo 


S. Archer 


Taste I (cont.) 


U.T. 


1956 V 
VI 
VI 
VI 


1956 V 
Vv 


VI 
VI 


1956 V 
, V 


VI 
VI 


29°936 
17°946 
22°929 
28-919 


15069 
30°045 
18-028 
22°970 


15°081 
30°056 
18-013 
22°955 


R.V. 
km/sec 
+ 377 
+ 13 
+ 1°6 
£. 3°0 


—$4°1 
—52°8 
—5$5°2 
$s" 


—S§I°r 
—52°8 
—52°4 
33'S 


Remarks 


Vol. 117 


The spectral classifications of brighter and fainter component (MKK) were 
determined by the writer as F8 V and F6 V respectively. 
a Cen: This double was observed for radial velocity before by Wright, Lunt 
and Wesselink. The star is of interest as being the only one at present for which 
radial velocities have contributed to accurate knowledge of the parallax. The 
main purpose for observing this pair was to provide a good value for the radial 


velocity at this epoch when Av is diminishing rapidly. 


ADS 10417 (36 Oph): No significant change in the velocities of either of 
the bright stars has occurred since the early Lick observations of 1908 (5). The 
radial velocity result for HD 155849 corroborates the known result from proper 
motion, that it is optical to the bright pair. In ADS this star is called C. 


ADS 6914 
brighter 


ADS 6914 


fainter 


aCenA 
a Cen B 


ADS 10417 N 
ADS 10417S 
HD 155849 
HD 156026 


Mlb 4 (br) 
Mlb 4 (ft) 


yCrAN 
yCrAS 


Taste II 
Mean values of results given in Table I. 
Epoch km/sec 
1956.39 +37°2 
1957-19 +41°9 
1956.39 +44°9 
1957-19 +45°8 
1956.43 —21°5 
41 —27°2 
1956.42 + o-2 
42 — 05 
43 —17°0 
.46 — ov! 
1956.44 + 16 
.46 2° 
1956.43 —54°3 
43 —52°4 


br-ft 


Av km/sec 
A-B +5°7 
N-S +0o°7 

—o'8 
—I ‘9 
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The star HD 156026 is known to possess common proper motion with the 
bright pair. ‘The radial velocity result shown in Table I confirms the star as 
a physical companion. The star shows H and K fairly strong in emission. 
In the General Catalogue of Radial Velocities HD 156026 is listed as identical 
with ADS 10417C, but this does not conform to the notation of the double star 
observers. 

Mlb 4: The velocity (+2-7 km/sec) of the brighter star agrees well with the 
Mt Wilson (6) value ( + 3-4 km/sec) obtained in 1920 which is about one revolution 
period ago in the visual orbit. At this phase Av is small. The nodal passages 
are predicted to occur in 1973 and 1981. ‘The fainter component shows Cat1H 
in emission. 

yCrA: Av is small at present. It is not possible to compare the results of 
the present paper with earlier velocities (5, 6) as previously only the combined 
light of the two stars was measured. The former blended velocity is in agreement 
with the present mean result. 

Acknowledgments.—1 wish to thank Dr A. D. Thackeray for useful advice 
during this work and especially Dr A. J. Wesselink for detailed comments on 
the manuscript, Dr W. H. van den Bos for supplying details of the systems 
studied and Dr A. Beer for allowing me the use of a measuring machine at 
the Cambridge Observatory. 1 should also like to thank the Radcliffe Trustees 
for the award of a Studentship for the period of my work at the Radcliffe 
Observatory. 
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A NEW TYPE OF PENCIL BEAM AERIAL FOR 
RADIO ASTRONOMY 


J. H. Blythe 


(Communicated by F. Graham Smith) 
(Received 1957 July 19) 


Summary 


Aerials of large aperture giving narrow pencil beams are difficult to realize 
at low frequencies. In the present paper a new method of synthesizing a large 
aperture is described. ‘This system, known as the Moving-T aerial, employs a 
long narrow aperture and a small moving aerial. When phase-switch records 
obtained at a number of different spacings are combined, the same result 
is obtained as would have been given by a very large conventional pencil 
beaminstrument. The equivalent aperture is bounded in one direction by half 
the length of the long aperture, and in the other by the greatest spacing used. 
The method is considered in some detail, and the limitations are discussed. 
A pioneer realization is described, which has provided a map of about half the 
sky at 38 Mc/s with a 2°-2 beam. The results will be given in a subsequent 
paper. 





1. Introduction.—It is widely appreciated that the study of the complex 
distribution of radio frequency emission in the Galaxy requires pencil beam surveys 
on a wide range of frequencies with beamwidths of the order of 1°. Among 
many features of interest is the narrow band 2°-3° wide along the galactic equator, 
partly due to thermal emission from Ht! regions having temperatures near 
10 000 deg. K, and partly due to non-thermal sources (1, 2). Radio and optical 
observations of individual clouds suggest that at low frequencies the total optical 
depth due to H11 regions is sufficient to produce appreciable absorption. It is 
therefore important to compare the detailed distribution at low and high 
frequencies. 

At low frequencies, taken to be less than 60 Mc/s, surveys have been carried 
out by Moxon (3) (40 Mc/s, 35° x 70° beam, along the galactic equator), Shain 
and Higgins (4) (18-3 Mc/s 17° x 17° beam, from declinations —12° to — 52°), 
and Higgins and Shain (5) (9-15 Mc/s, 29° x 29° beam centred on declination 
— 32°). 

A description is given in this paper of a new means of obtaining a pencil 
beam, which has been used to carry out a survey of more than half the celestial 
sphere with a 2°-2 beam at 38 Mc/s. The results will be given in a subsequent 
paper. 

2. The problem of achieving sufficient resolving power at low frequencies.— 
The beamwidth Aé of an aerial is related to its maximum dimension d by the 
expression 


n 
Aé = 7 
which shows that to achieve a beamwidth of 1° at 38 Mc/s the aerial aperture 


must extend for 1500 ft. This rules out a steerable paraboloidal reflector, and 
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even a flat phased array would represent a considerable undertaking. It is 
however possible to use small aerials to give the same resolving power as large 
structures, although the sensitivity obtained may not be as great. A more general 
account of these methods will be given by Hewish and Ryle (6). 

The simplest example is the variable spacing interferometer. When the 
records obtained successively at all spacings up to a certain value are added, 
taking proper account of phase, the result is the same as that which would be given 
by an aerial with aperture as wide as the greatest spacing used. The arrangement 
has been widely used to study the one dimensional distribution of brightness 
across the Sun (7). In a similar way, Ryle and Scheuer measured the profile 
of the narrow band of emission along the galactic equator using an interferometer 
axis such that the high resolving power was arranged to be perpendicular to the 
galactic equator (1). 

A high resolving power in two dimensions requires aerials with physical 
extent in two dimensions and by analogy with these earlier experiments we 
can see that such aerials may be simulated by the addition of records in which 
both the spacing and the direction of the interferometer axis are varied (8). 
If for example we wish to simulate a uniformly excited square aperture, we can 
do it in the manner of Fig. 1, where A and B are small aerials, the two components 
of an interferometer, which occupy in turn all the possible pairs of positions in 
the square. 





Fic. 1.—Aperture synthesis. The small squares represent aerial positions. 


it may be seen from Fig. 1 that many combinations of spacing and direction 
occur more than once (e.g. AB, CD), and that the frequency of occurrence 
decreases linearly with increasing spacing (thus EF is unique). Clearly AB 
and CD give similar records, and the only advantage gained by using both is that 
of increasing the sensitivity. If the signal to noise ratio is adequate it is therefore 
unnecessary to use more than one record from each similar arrangement of the 
two aerials; the record amplitude obtained with each arrangement must then 
be given a weight proportional to the number of times the arrangement occurs in 
the desired aperture. The deterioration in the signal to noise ratio of the record 
depends on the square root of the number of times the arrangement occurs. 

The present system represents one application of this principle. A long 
aperture aligned in the E-W direction is used in conjunction with a small aerial 
which is moved to positions ona N-S line through the centre of the long aperture, 
as shown in Fig. 2(a). In this case all the possible interferometer arrangements 
of Fig. 1 with the same N-S spacing contribute simultaneously, and it is necessary 
to adjust the excitation along the long aperture to provide the appropriate weights 
for the various arrangements. For a uniformly excited square aperture we 
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have seen that the frequency of occurrence decreases linearly with spacing, so 
that a “‘triangular grading” of excitation is required. It is similarly necessary 
to decrease the sensitivity of the receiver linearly as the N—S spacing is increased. 

When the records are added in the way described in the next section, we have 
contributions from all the spacings which are present in the area outlined in 
Fig. 2(c), which therefore gives the area of a square aperture providing a pencil 
beam of the same size. 

This system may be compared to the Mills Cross (g), as shown in Fig. 2 (6), 
in which two perpendicular long apertures are connected to a phase switch 
receiver. In this case all the possible interferometer spacings and directions 
contribute simultaneously, and it is necessary to adjust the excitation along the 
arms of the cross to provide the appropriate weights for the various arrangements. 
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Fic. 2.—{a) Moving-T aerial; (b) Mills Cross aerial; (c) Conventional aperture giving equivalent 
resolution. 

3. The method of analysis.—We must now consider how the information 
obtained from each separate interferometer record is combined to give the 
brightness distribution across the sky as seen by a pencil beam instrument. 
Let us examine the way in which the receiver output depends on the true bright- 
ness distribution and the spacing. Ryle has shown (10) that the output of 
a phase-switch receiver due to a signal from a given direction is proportional to 
the geometric mean of the gains in this direction of the two aerials, multiplied 
by an interference factor whose periodicity depends on the spacing d. We use 
the concept of the equivalent area A of each aerial, related to the gain G by 

7A 
c=". 

The gain is a function of direction, and we take modified polar coordinates as 
shown in Fig. 3. If A,(8,¢) and A,(6, ¢) are the equivalent areas of the long 
aperture and the small aerial respectively, and d is the spacing which is in the 
N-S direction, then the output recorded at any given time will be proportional to 





C(d) =const | | T(6, ¢)V A,(9, ¢)A,(8, ¢) cos (eee) sin6 ddd. (1) 


Now A,(9, ¢) and A,(8, $) vary in the same way with declination, that is with ¢, 
so that it is convenient to put 


K A,(6, $) = 4,(0)A(¢) 


K, 406, #)= Al) A(9) (2) 
where A(o) = A;(o) = A,(o) = 1. 
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Both A,(6) and sin@ may be taken as constant over the few degrees for which 
A,(9) is appreciable. Integrating with respect to 0 we find 
and sing 


C(d)=K | A()T(#) cos es dd (3) 


where 7(¢) is the “‘ smoothed” value of T(6, ¢), (12), and K isa constant depending 
on the beamwidth and forward gain of the interferometer. On making the 
substitution 

S=sing 
the equation can be written 





C(d)=K | A'(S)T"(S) cos Ce =) dS (4) 


where 


A'(S)= A(g)/cos$ 
T'(S)=T(¢). 


This equation shows that C(d) gives one Fourier component of A’(S)7’(S). 
The corresponding sine term is found from the output S(d) with an extra quarter 
wave inserted in the path to one aerial. The discussion in Section 2 implied 
that it is necessary to take records at spacings of an integral number of aerial widths. 


Zenith 


-8 x 


90-68 | N 

















Fic. 3.—Coordinate system. 


In the present case the aerials are narrow, and the polar diagram extends from 
the northern to the southern horizon. ‘This means that A(S) extends over the 
range —1<S<1, and only spacings d=mnA/2 (mn integral) are needed. The 
record amplitudes are synthesized according to the expression 


T(S)A(S)=Taaot & >. (ca) cos ee ) + S(d)sin (“)} : & 


d=A/2 





The zero spacing term 7.) corresponds to the physically unrealizable super- 
position of the two aerials. However, it may be found from the total power 
output of the long aperture. Since the long aperture used alone contains greater 
spacings than have been used (see Fig. 2 (a)) it is necessaryin principle to re-arrange 
the excitation so that the same spacings are present in the same proportions; 
half the aperture would then be used with uniform excitation. In practice, the 
slowly varying nature of the record obtained showed that it was sufficiently 
accurate to ‘‘smooth” the record obtained without altering the excitation. 
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This synthesis uses all the values of receiver output at a given sidereal time 
for all spacings, and it gives the brightness distribution along a line of constant 
right ascension. A map of the sky is built up by repeating the synthesis at a 
sufficient number of sidereal times, corresponding to the E-W resolving power of 
the long aperture. 

The constant K in equation (5) may be found by measuring the amplitude and 
width of the receiver output from the intense point sources of known flux density 
in Cygnus and Cassiopeia. Observations by Adgie and Smith (11) have shown 
that the broad distribution of brightness is nearly the same at 38 Mc/s and 80 Mc/s. 
The primary declination polar diagram A’(S), very broad in this experiment, was 
therefore checked by comparison of the broad features with the survey at 80 Mc/s 
carried out by Baldwin (12). 

4. Theoretical limitations of the system.—F or an object to be observable with 
a radio telescope it must produce a signal large enough to overcome a variety of 
obscuring factors which include the inherent noise fluctuations, the signal due 
to other sources within a beamwidth, and the signal due to strong sources outside 
the main beam. It was shown in Section 2 that the moving-T aerial will produce 
a map on which the various features do not appear with as large a signal to noise 
ratio as they would on a map produced in the same time by the filled aperture with 
identical polar diagram. For the present structure the reduction is about five 
times. 

To calculate the size on the map of the spurious signals due to thermal 
fluctuations, let 7,, be the r.m.s. value on each record, expressed as an equivalent 
input temperature. Then if N spacings are used, and the sensitivity factor of 
equation (5) is 1/K, the r.m.s. level of spurious signals on the final map will be 
T,V N/2x1/K. 

Bright sources may produce signals in the wrong place, both through the side 
lobes of the long array and owing to errors in measurement when they appear 
in the main beam of the array. When a bright source is in a side lobe its zenith 
angle ¢ may take a range of values, and consequently it may produce spurious 
features on a locus one beamwidth wide through the star passing round the pole. 
Near the star this locus lies along a line of constant declination. Using the intense 
source in Cassiopeia the greatest side lobe level was shown to be 3} per cent of the 
main beam response, and the general level within 30° of the forward direction was 
about 1 per cent. 

Errors in gain and phase of the recording sensitivity produce signals which 
are not genuine, but which appear as sources with the same time of transit as 
a bright star. A convenient parameter for estimating the effect of amplitude and 
phase errors is 

A? = (Ad)? + (AA)? 
where A¢ and AA are the r.m.s. phase and fractional gain errors in the recording 
sensitivity. Errors of this type give side-lobes at random declinations, having 
r.m.s. value A x VN. 

__ This result enables us to predict only the kind of response we expect from 
single bright radio stars away from the main beam. This is of importance near 
very few sources, and the major limitation over most of the sky is the production 
of false features by errors in the measurement of typical signals arising from the 
superposition of the effects of many sources. Since the sources which give the 
largest outputs only do so when the aerial spacing is small, apart from a few bright 
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discrete sources, the largest errors correspondingly appear on the map only as 
large scale features. Thus we may conveniently divide the amplitude-spacing 
curve into two parts: the small spacings of high intensity and the larger spacings 
with considerably smaller intensity. The spurious features on the map will 
have r.m.s. values 7',, and T,,. related to the receiver outputs at various spacings 
and the gain and phase errors by 

I 


Ta=R { 2A? > [C(d)?+ say} 


d= N(4/2) | 
T= RWS [Cld*+ Stay} 
d =24 

The first group will give rise to spurious features of large size of surface 
brightness up to about 27,,,, and the second group will give spurious “‘ discrete 
sources” up to 27°». 

5. Aerial design.—Full wave dipoles were used, loaded with capacity to increase 
the phase velocity to 2c, so that they resonated in a half wave mode. A reflecting 
screen A/4 wide was placed at a distance A/10 below each dipole. The moving 
aerial consisted of a single dipole, and 48 were used for the array. It has been 
indicated above that a triangular grading of excitation is necessary to simulate 
the beam of a uniformly excited square aperture. It is possible to obtain a 
somewhat narrower beam at no extra expense in side lobe level by using certain 
other ways of grading the excitation. That used in practice was approximately 
Gaussian, truncated when it had fallen to a small value. The feeding system, 
which enabled the current in each dipole to be controlled independently, is shown 
in Fig. 4. 
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Fic. 4.—(a) Branching system used to feed the six groups of dipoles; (b) Feeding arrangements 


within a group of eight dipoles. The lengths of the twin feeder lines are shown; the trans- 
formers (T) lead to a balance-to-unbalance transformer (B). 





The dipoles were fed in groups of four, and each dipole was connected to 
a twin feeder by a quarter-wave whose characteristic impedance was made variable 
to control the current relative to the other dipoles in the same group. ‘Two groups 
of four dipoles were brought to one balance-to-unbalance transformer, the relative 
impedances and the matching into the transformer being controlled by means of 
another pair of quarter-wave transformers. A branching feeder system was 
used to connect the resulting six coaxial outputs to the receiver through equal 
lengths of cable. The power division and the matching at the junctions were 
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This synthesis uses all the values of receiver output at a given sidereal time 
for all spacings, and it gives the brightness distribution along a line of constant 
right ascension. A map of the sky is built up by repeating the synthesis at a 
sufficient number of sidereal times, corresponding to the E—W resolving power of 
the long aperture. 

The constant K in equation (5) may be found by measuring the amplitude and 
width of the receiver output from the intense point sources of known flux density 
in Cygnus and Cassiopeia. Observations by Adgie and Smith (11) have shown 
that the broad distribution of brightness is nearly the same at 38 Mc/s and 80 Mc/s. 
The primary declination polar diagram A’(S), very broad in this experiment, was 
therefore checked by comparison of the broad features with the survey at 80 Mc/s 
carried out by Baldwin (12). 

4. Theoretical limitations of the system.—F or an object to be observable with 
a radio telescope it must produce a signal large enough to overcome a variety of 
obscuring factors which include the inherent noise fluctuations, the signal due 
to other sources within a beamwidth, and the signal due to strong sources outside 
the main beam. It was shown in Section 2 that the moving-T aerial will produce 
a map on which the various features do not appear with as large a signal to noise 
ratio as they would on a map produced in the same time by the filled aperture with 
identical polar diagram. For the present structure the reduction is about five 
times. 

To calculate the size on the map of the spurious signals due to thermal 
fluctuations, let 7, be the r.m.s. value on each record, expressed as an equivalent 
input temperature. Then if N spacings are used, and the sensitivity factor of 
equation (5) is 1/K, the r.m.s. level of spurious signals on the final map will be 
T,,V N/2x1/K. 

Bright sources may produce signals in the wrong place, both through the side 
lobes of the long array and owing to errors in measurement when they appear 
in the main beam of the array. When a bright source is in a side lobe its zenith 
angle ¢ may take a range of values, and consequently it may produce spurious 
features on a locus one beamwidth wide through the star passing round the pole. 
Near the star this locus lies along a line of constant declination. Using the intense 
source in Cassiopeia the greatest side lobe level was shown to be 3} per cent of the 
main beam response, and the general level within 30° of the forward direction was 
about 1 per cent. ; 

Errors in gain and phase of the recording sensitivity produce signals which 
are not genuine, but which appear as sources with the same time of transit as 
a bright star. A convenient parameter for estimating the effect of amplitude and 
phase errors is 

A? =(Ag)?+ (AA)? 
where A¢ and AA are the r.m.s. phase and fractional gain errors in the recording 
sensitivity. Errors of this type give side-lobes at random declinations, having 
r.m.s. value A x VN. 

This result enables us to predict only the kind of response we expect from 
single bright radio stars away from the main beam. This is of importance near 
very few sources, and the major limitation over most of the sky is the production 
of false features by errors in the measurement of typical signals arising from the 
superposition of the effects of many sources. Since the sources which give the 
largest outputs only do so when the aerial spacing is small, apart from a few bright 
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discrete sources, the largest errors correspondingly appear on the map only as 
large scale features. Thus we may conveniently divide the amplitude-spacing 
curve into two parts: the small spacings of high intensity and the larger spacings 
with considerably smaller intensity. The spurious features on the map will 
have r.m.s. values 7,,, and 7’, related to the receiver outputs at various spacings 
and the gain and phase errors by 

I 


Lag = 2 2 |X? 
KAS (Cla + s(ay)} 


== 4 


T11= 


I d= N(i/2) 1/2 
Tre= eth? > [C(dP+S(d)]}> . 
K d =r 


The first group will give rise to spurious features of large size of surface 
brightness up to about 27,,,, and the second group will give spurious “ discrete 
sources” up to 27°. 

5. Aerial design.—Full wave dipoles were used, loaded with capacity to increase 
the phase velocity to 2c, so that they resonated in a half wave mode. A reflecting 
screen A/4 wide was placed at a distance A/10 below each dipole. The moving 
aerial consisted of a single dipole, and 48 were used for the array. It has been 
indicated above that a triangular grading of excitation is necessary to simulate 
the beam of a uniformly excited square aperture. It is possible to obtain a 
somewhat narrower beam at no extra expense in side lobe level by using certain 
other ways of grading the excitation. That used in practice was approximately 
Gaussian, truncated when it had fallen to a small value. The feeding system, 
which enabled the current in each dipole to be controlled independently, is shown 
in Fig. 4. 
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Fic. 4.—(a) Branching system used to feed the six groups of dipoles; (b) Feeding arrangements 
within a group of eight dipoles. The lengths of the twin feeder lines are shown; the trans- 
formers (T) lead to a balance-to-unbalance transformer (B). 


The dipoles were fed in groups of four, and each dipole was connected to 
a twin feeder by a quarter-wave whose characteristic impedance was made variable 
to control the current relative to the other dipoles in the same group. ‘Two groups 
of four dipoles were brought to one balance-to-unbalance transformer, the relative 
impedances and the matching into the transformer being controlled by means of 
another pair of quarter-wave transformers. A branching feeder system was 
used to connect the resulting six coaxial outputs to the receiver through equal 
lengths of cable. The power division and the matching at the junctions were 
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controlled by stubs. The E-W polar diagram was checked by observation of 
the bright source in Cassiopeia with a phase-switch receiver, using the array 
and a dipole as an interferometer with a large E-W spacing to remove the 
background. The beamwidth of the complete system was 2°-16 to half intensity. 

Survey. (1) The method of observation.—We have seen that the information 
required is one total power record giving 7'jW» in (5), and the output of the phase- 
switch receiver (a) with equal cables, (b) with an extra quarter-wave in one cable, 
for all spacings d=n(A/z) throughout the whole 24 hours. Thus ideally the 
movable aerial is first placed A/2 from the centre of the array in the N-S direction, 
and the output C(A/z) recorded for 24 hours. ‘Then an extra quarter-wave is 
inserted in one cable, and the output S(A/2) recorded for the next 24 hours. 
Then the movable aerial is placed at a distance A from the array, and the outputs 
recorded similarly for the next 48 hours at this spacing. ‘The procedure is then 
repeated for all d=n(A/2) up to the largest spacing N(A/2), the total observing 
time being 2N days. 

In practice observations were made initially for one night at each spacing. 
Both sine and cosine records were obtained simultaneously, using a device con- 
structed by J. H. Thomson which introduced a quarter-wave into one side of the 
interferometer at intervals of approximately 2 minutes, giving adequate sampling 
of both C and S records. This device employed coaxial switches driven by a 
synchronous motor. Calibrations were made weekly to ensure that the receiver 
sensitivity stayed constant, but unfortunately there proved to be serious variations 
in the amplitude and phase of the signal from the movable aerial due to the effect 
of moisture on some of the twin-feeders for this aerial. The gain and phase of 
the entire system were checked daily on the record fromthe bright stars. Although 
the observations were made near sunspot minimum, man-made interference 
made daytime records useless, so that it was necessary to repeat each spacing 
several months later. Most spacings were also repeated at an intermediate date. 

(2) The use of electronic computation.—Before describing the manner in which 
the computing was performed, it is interesting to consider the size of the task. 
At a given right ascension, observations were made at 38 spacings (up to 184A) 
giving 76 numbers which must be transformed according to equation (4). The 
record was read every 4 minutes, that is a total of 360times. This gives altogether 
about 2 x 10° operations without considering the insertion of sensitivities, grading 
or polar diagram. Through the kind cooperation of the Director of the Cambridge 
University Mathematical Laboratory the computing was handled by the 
EDSAC-I, a large electronic computer (13). A programme was devised which 
performed the entire computation in 15 hours. The information provided by 
the output of the computer was the declination of certain contours, rather than the 
temperatures at certain declinations. This eased the drawing of the map, and 
took no extra time to compute. 

Input to the EDSAC is by means of punched tape, so that a major task was 
the preparation of error-free data tapes. One difficulty was that owing to lack 
of storage space in the EDSAC it was necessary for the intensities for all spacings 
at a given right ascension to be arranged consecutively on the tape so that the 
calculation could be carried out for that R.A. before the next set of intensities 
was read in. This involved rearranging the data because each record gives the 
intensity for a given spacing at all right ascensions. This was done by writing 
down the intensities from the records in rows and punching the tapes from the 
numbers read down in columns. 
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(3) Summary of results—Observations have been made at all right ascensions, 
and are generally complete up to a spacing of 18}A. A map has been plotted 
from declination 70° to — 20°, ignoring the times of transit of the intense sources 
in Cygnus and Cassiopeia, where the side-lobe confusion is great. The map is 
presented in asubsequent paper. The total area of sky is just over half the celestial 
sphere. 

The most important source of error lay in the fluctuation of gain and phase 
of the moving aerial, which had the effect of giving large side-lobes in declination. 
The value of A? was estimated as 0-044 from the scatter of the amplitude and phase 
of the signals from the bright stars. This means that the r.m.s. level of side-lobes 
due to this is 4 per cent. A grading of receiver sensitivity with spacing was 
adopted which gave in the absence of errors a first side-lobe of 4 per cent the others 
being smaller. At the transits of bright sources we can therefore expect side-lobes 
over a wide range of declination of the order of 4 per cent of the source intensity. 
The more general spurious features depend on 7,,, and 7,2, which for typical 
regions away from the galactic plane have the value 700 deg. K and 500 deg. K. 
Noise fluctuations and random drifts have also been included in the latter. The 
background temperature in these regions is around 10000 deg. K. In spite of 
practical difficulties in this first realization of the moving-T aerial, a satisfactory 
survey with a beamwidth of 2°-2 has been carried out over a large region of sky, 
and many interesting features have been revealed. 
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Summary 
A pioneer tealization of a new aerial principle has produced a contour 
map of the sky at 38 Mc/s with a resolving power of about 2}°. The list 
of sources includes a number with appreciable angular diameters which have 
not previously been recorded. Some results for H 11 regions are also given. It 
is suggested that sources Nos. 18 and 24 represent regions where gas, dust, 
magnetic fields and relativistic particles occur together. 





1. Introduction.—A survey of the galactic radio emission at a frequency of 
38 Mc/s has been carried out using a new form of pencil beam “ radio telescope’’. 
High resolution was achieved by means of a new interferometric technique which 
has been described in a previous paper (1), referred to here as Paper I. The 
results are presented in the form of a contour map (Fig. 1) together with a list of 
reliably observed sources and features (Table 1). Some details of the observations 
are given in Section 2; the remaining sections comment briefly on the main 
features revealed by the survey. 


2. Details of the observations.—T he observations cover most of the sky north of 
declination — 20°. Some areas have been omitted where the results are untrust- 
worthy due to the proximity of intense sources, so that the total area covered by 
the map in Fig. 1 is about half the celestial sphere. The width of the beam to half 
intensity points is 2°-2 in the east-west direction, and in the north-south direction 
it varies according to the zenith angle from 2°-3 at declination 52° to 7°-4 at 
declination —20°. The observations were made from June to December, 1955. 

The temperature plotted in Fig. 1 is the synthesized aerial temperature 
T’(sin 6) obtained by means of equation (5) of Paper I (where 0 is the zenith angle). 
It is equal to the aerial temperature that would be measured with a conventional 
pencil-beam instrument having the above beamwidth, and will be referred to as 
aerial temperature, or 7.,. 

The uncertainties in the map of the sky require careful consideration with the 
new type of aerial. Reference should be made to Paper I for details, but in the 
present map we may expect the following type of error :— 


(i) At any point there may be a random error in the background temperature 
of about 200 deg. K. 


(ii) Individual intense sources produce spurious features at positions with the 
same time of transit (approximately along hour circles through the sources) and 
also on a locus through the star passing round the pole (near the source this locus 
lies roughly on a line of constant declination). These spurious features will 
appear with about 4 per cent of the intensity of the source itself. 
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TABLE | 
Position (1950) Aerial Size Flux density 
temperature 
No. R.A. Dec. T, (deg. K). w.m~? (c/s)—! x 1078 
h m 
I 00 o3 —16 3000 s XF 350 
2 oo 33 -1 2000 P 150 2Cso 
3 00 35 —I10 5000 24° x24 480 An extension joins 
to No. 1. 
4 ©O 40 40 2000 3° x3 390 IAUooN4A 
5 CO 53 26} 2000 P 150 2C72 . 
6 co 54 —1} 1500 P 110 2C73 Superposed on 
extended feature 
7 or of —133 3500 5° xP 600 
Ol 05 13} 2000 P 150 2Cg94 Superposed on 
larger source 
9 oO: 1 36 1700 2° x24 200 2Cg99 
10 o1 21 —13 1500 P 110 2C122 
sn “ee 3to5 ™~jz000 21° x6° ™~3000 
os 655 ' 
12 o1 34 20% zooo=—s«r$° x 13° 150 
tt Fu? 7k to 15 —~2000 11° x6 ~ goo 
“a 
14 03 00 154 2000 xs 210 
15 03 15 414 5000 - x? 365 IAU 
o3N4A 
t Me * een $3 ~anee Structure along gal. 


03-20 equator, 
/=105° + 112° 


03 Po Along gal. equator, 
04 ~=«10 sadhana aanied [=116° -, 120° 


18 a to = gto 19 ™4000 rxs° oo “ ; possibly 
9 «©6308 104 2000 P 150 2C349 Superposed on 
larger source 
20 04 14 384 3000 P 220 2C379 
21 04 33 29} 6500 P 80 _— 
J va 5 4 o4N3A 
‘ 1 IAU 
22 04 56 46 5500 3 15 440 o4NgA 
23 05 of 38 5000 P 370 2C440 
24 O§ 40 I —~3z000 10° X5 ™1300 BWSS-D 
25 o6 o2 15 2500 23° x2} 320 Along gal. equator 
1= 160° —. 164° 
26 05 #651 26 2500 P 180 
27 06 03 —8} ~8000 24° xP ™~ 300 2C520 
28 « to °° 15 to20 ™~2500 Belt A 
08 — 00 
29 06 14 22} 9500 P 730 ime 
30 06 24 —6} 3000 P 220 2C553 
31 06 42 —6} —~4000 Along gal. equator 
1= 184° — 187° 
2 07 12 —15} —~2000 l= 194° + 197° 
33 07: «+20 «-—8 5000 exe 520 
IAU 
34 08 07 484 2800 P 200 o8N4A 


35 o8 17 4} 2000 2° X4 270 
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Fic. 1.—Galactic radio emission at 38 Mc/s (each un 
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Position (1950) Pe. Size Flux density 
No. R.A. Dec. w.ni-? (c/s)-! x 10% 
wil IAU . 
36 09 #814 —13 19000 P 1400 ogS1A S,09-1 
37 o9-~=s«18 454 1700 a 130 2C805 
38 o9 ©6640 14 2000 - x 2$ 300 
39 09 «43 —12$ = 4500 2h x 8° goo 
49 09 43 8 1700 P 130 2C843 
41 09 «859~—l28 2000 P 150 2C855 
42 10 05 19 2000 3° «xX 4° 400 ; 
43 10 2-8 4500 ae 1700 S 10-1 
44 10 47 164 2000 s* xan’ 270 S 10+1 
45 1 12 — 63 2500 ” <9 180 
46 12 25 — 2} 14000 a° xP 1100 
47 12 36 — 24 —~Jooo 10° x11° ~5100 BWSS-H; S 12-0 
48 12 32 —124 g* Seg" 
IAU 
49 12 57 464 2000 P 150 12N4B 
50 12 §7 —17 —~6000 < =<? 630 2C1073 
gr 13 © 6 30000 2g x 24° 350 
52 13 03 — 2h 2500 2° xP 200 
IAU 
53 14 10 5243 2000 P 150 14NsA 
13 Og, 15 width 
54 16 00 20 6° to 20° Belt A 
17 53° 4 length 
> 90" 
55 14 43 — 4 sooo 7° x 7° 1700 
56 a os to be 3000 Belt B 
57 15 OF 26 ™~3000 P ~220 2C1259 
IAU : 
58 15 12 74 4000 P 300 rsNrA | 
59 15 13 —12 7000 3° x6" 1300 S 15-1 ; 
60 15 38 13 ~1500 4° xP 180 2C1315 \ Structure of | 
es 2 “£6 8 2000 2° xP 180 Belt A { 
62 16 o1 1} 3000 P 220 2C1359 j 
63 16 oF — 6 2000 P x¢s° 490 { 
64 16 24 8 2000 Oe? 220 2C1390 : 
65 16 27 392 2500 P 180 2C1402 
66 16 35 15 4000 Structure of Belt A 
IAU 
67 16 48 5 21000 P 1600 16NoA 
68 17 0% 11} 3000 2$° xP 300 2C1454 
69 17 16 —1 8000 P 600 2C1473 
7Oo 17 52 —8 20000 34° x 4° 3200 
71 18 24 —113 50000 - xe 6000 
72 18 30 9} 8000 24° xP goo 
73 18 54 14 6500 «14° xP 520 
74 #18 58 54 28000 3° x2° 3000 2C1607 
20 10, —974h i 
7 we “ee 124 ~10000 Belt A. Includes S 21-1 / 
76 20 30 — 7} 4000 6° xP 620 Structure on Belt A 
77 20 42 63 2500 2°. xP 240 2C1716 
78 20 44 504 8000 Pr xX3° 800 2C1725 
79 20 47 30 4500 1 x1° 380 
80 21 20 52 2500 Along gal. equator 
1=61° 
8: 21 24 —8 3000 24° xP 330 Structure of Belt A 
82 22 20 17 4000 Complicated object. 
The intense sources [AU Nos o5N2A, 12N1A, 19N4A, 23N5A which were used as 
calibration sources have been omitted from this list. 
The sizes are given East-West diameter first, then North-South diameter. P indicates. 
there is no detectable broadening of the polar diagram. 
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(iii) The individual “side-lobes” of numerous sources of relatively small 
intensity combine to give a random disturbance of the true distribution. Some 
features of small angular diameter with aerial temperature up to 500 deg. K, and 
some larger features up to 700 deg. K, may derive from this experimental error. 

It will be seen that with the adopted contour interval of 2000 deg. K appre- 
ciable errors can occur in the map due to effects (ii) and (iii). Contours which are 
uncertain because of the effect of the bright parts of the Galaxy or intense discrete 
sources are shown dashed. Some sources of low surface brightness do not show 
clearly on the map, but Table I includes any such sources which appear to be 
reliable. 


3. Galactic Structure 


3-1 Halo distribution.—Westerhout and Oort (2) have drawn attention to the 
considerable amount of radiation from high galactic latitudes. It has been shown 
by Baldwin (3) that most of this radiation originates in a spherical halo round the 
Galaxy. Inspection of Fig. 1 in relation to the model given by Baldwin (3) has 
shown that this component accounts for most of the radiation at 38 Mc/s from high 
galactic latitudes. 


3.2. Galactic plane.—In examining the radiation from regions within about 
15° of the plane we may bear in mind that at higher frequencies it has been asso- 
ciated with several different types of source. ‘Thus Westerhout and Oort (2) 
suggested a source distributed in the same way as mass in the Galaxy. Ionized 
hydrogen is observed optically within about 2° of the galactic equator, and may 
be expected to produce a narrow band of emission at high frequencies, and 
absorption at low frequencies. Ryle and Scheuer (4) observed such a band of 
emission at 80 Mc/s and 210 Mc/s, but the position may be complicated by the 
presence of non-thermal sources concentrated near or towards the plane (5). 























250,000 _ 
deg K - Tt r 
200,600 q 
/=0° 
150,000 
150,000 deg K 4 
100,000 100,000 
50,000 50,000 
—_ i ™ 1 
-20° -10° -20° -10° 0° 
50,006 | L : 
deg K /=60 
40,000 L 
30000 30,000 + 
degK ™ 
20,000 pers 
== _ 
10,000 + 10,000 + 
4 b 
— iL i j L 1 1 | 
-20° -10° °° 10° 20° -20° -10° 0° 10° 20° 


Fic. 2.— Sections across the galactic equator. The dashed curve is the profile expected from the halo 
distribution and a source distributed like mass in the galaxy. 
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The absence of a band of absorption in the sections across the galactic equator 
in Fig. 2 confirms the presence of non-thermal sources, for the total absorption 
through the Galaxy would be considerable if a large part of the band of emission 
at higher frequencies were emitted by H11 regions. 

Profiles expected if the radiation originated in a halo distribution and a source 
distributed like mass are also shown in Fig. 2, based on the curves given by 
Westerhout and Oort (2). Consideration of Fig. 1 shows that the observed 
distribution is very irregular. The profiles for the mass distribution have there- 
fore been adjusted, both by applying different intensity scale factors and by 
different translations in latitude. We note that radiation associated with the 
galactic plane does extend to latitude 15° at most longitudes. 

Since the band of emission is highly concentrated towards the galactic plane, 
and the radiation comes from distant parts of the galaxy, it may be used to locate 
the galactic pole. On the basis of a survey at 600 Mc/s, Piddington and Trent 
(6) have derived coordinates b=89°-1, /~330° relative to the Lund pole, with a 
solar distance of 56 parsecs. The position of the band varies irregularly about 
that expected on the basis of the above values, so that the present survey lends 
support to these coordinates. 

3-3 Hitregions.—H 11 regions are importantin the interpretation of the structure 
along the galactic equator, and a search was made for individual bright regions 
which might be observed at 38 Mc/s by their absorption of the background 
radiation. Fifteen regions observed optically are given in Table II: these lie in 
the areacovered by the present survey and have diameters greater than}°. Smaller 
regions do not fill a sufficient fraction of the beam to produce a detectable signal. 


Taste II 
Observations of H 11 regions 
T, wii T, T, 


R.A. Dec. Diam. : ; T EM 
units of units of 

h m . 1000 deg. K 1000 deg. K 
02 29 61 3° o+2 17 <0°3 < 1000 
03 «48 36 2°°4 x 0° 8 o+1 II ae 
o5 si rs 
oe a nn 
05 «(33 II 6°°8 2 13 21 > 3000 
06 833 ¥ 1°°2 o+ 16 Sine — 
05 34 — 6 i x7" 4 o+1 15 <o'2 < 600 
07 05 —II xe 2 33 0°09 300 
16 36 —10 10° 7 22 06 1800 
18 13 —12 .“*s 
18 14 —17 °° 8 
18 16 —12 1°°4 22 180 0*4 1300 
18 16 —14 or 
18 18 —16 °°'9 
22 18 55 1°°s o+1 19 <0%73 <1000 


Optical data were compiled from : 


Bok, Bester and Wade, Daedalus, 86, 9, 1955. 
Gum, Mem. R.A.S., 67, 155, 1955. 

Sharpless, Ap. F., 118, 362, 1953. 

Sharpless and Osterbrock, Ap. F., 115, 89, 1952. 
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Five regions of low brightness appear in the survey: ten of the fifteen H 1 
regions are contained in these regions, some containing more than one H 11 region 
owing to the limitation of resolving power. The remaining five regions were not 
detected. 

Details of the regions appear in Table II. T, is the aerial temperature 
observed in the direction of the region, and 7, is the aerial temperature observed 
in local directions. If we assume that 7, is the temperature that would be 
observed in the absence of the region, then (7;— 7) is the reduction of aerial 
temperature caused by the region. If we make the further assumption that each 
region is a uniform sphere with angular diameter given by the optical results, then 
the fraction of the beam which is filled can be calculated. Hence the total 
absorption or optical depth of the region may be found if its temperature is known ; 
in all cases the theoretically expected value of 10000 deg. K was taken. By this 
means limits may be placed on the emission measure, EM, which is proportional 
to the optical depth. The emission measure is defined by 


EM= J N22 ds 


where the distance s is measured in parsecs, and the electron density N, is the 
number per cc. Values of emission measure derived in this way are given in 
Table II. 


4. Localized sources.—F or those sources in Table I which have been observed 
previously the relevant survey numbers are shown. ‘The area covered by the 
survey of radio stars of Shakeshaft, Ryle, Baldwin, Elsemore and Thompson (7, 
referred to here as 2C) includes the area of the present survey, and stars down to 
a low sensitivity limit were detected. Comparison of the two surveys provides 
a check on the expected sensitivity limit. In view of the effect (iii) described in 
Section 2, there will be spurious stars of size about 500 deg. K, so that stars of 
aerial temperature less than 1700 deg. K have been ignored except for three at 
1500 deg. K which agree with other surveys. An aerial temperature of 1700 deg. K 
with the present aerial corresponds to a flux of 75 x 10-** mks units at 81-5 Mc/s, 
assuming an average spectral index of —o-7. There are 4o stars in the 2C list, 
with flux densities greater than this value, which lie in the area covered by the 
present survey. Of these 30 have been detected, and 4 lie in areas of the sky 
where the present sensitivity limit should be increased due to the presence of 
stronger signals, or lie in complicated regions where identification is difficult. 
The agreement is therefore 80 per cent. 

A number of the sources in Table I give records which are appreciably wider 
than the natural beamwidth. This broadening enables an estimate to be made of 
the source diameter at half intensity, if it is assumed to have a particular shape ; 
values assuming a Gaussian source shape are shown in Table I. For diameters 
less than 3° the broadening produced 1s small and the values are correspondingly 
inaccurate. 

A short list of broad sources was given by Bolton, Westfold, Stanley and 
Slee (8). Four lie in the area covered by the present survey and all appear in 
Table I. From the results of a pencil beam survey Shain and Higgins (9) 
recorded eight sources in the area of the present survey, of which seven agree. 
The other, S12-1, may possibly have been a blend of sources in a complex region 
which is now resolved, 
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Since individual galaxies have been identified as radio sources it is of interest 
to see if any of the broad sources originate in clusters of galaxies. Shane and 
Wirtanen (10) have given a contour map of the number of galaxies per unit area. 
This shows clustering, but the coincidences between the radio sources and the 
clusters are not significant. 

The present survey has revealed new structure in the extended feature denoted 
as Belt A in Table I (e.g. Nos. 18, 28, 54 and 75). Baldwin (11) has pointed out 
that this belt extends roughly round a great circle, and the present survey alsoshows 
an enhancement of radiation where this belt cuts the galactic equator at /= 162° 
(No. 25). 





| a eS ey kiloparsecs 
Galactic centre 
Fic. 3.—Comparison of sources of 38 Mc/s close to the galactic plane and the distribution of H1 
emission. Radiants represent the directions of sources: irregularities in H1 distribution have 
been marked with crosses. 





A number of the large diameter sources in Table I (e.g. Nos. 16, 17, 25, 31, 32; 
71, 74 and 8o) lie closely along the galactic plane, and possibly represent large 
scale irregularities in emission rather than localized sources. Radiants representing 
these sources have been plotted in Fig. 3, together with contours of the distribu- 
tion of the H1 emission in the galactic plane as derived by Van de Hulst, Muller, 
Oort (12). ‘The H1 contours show structure in the two large spiral arms, and an 
attempt has been made to distinguish the larger irregularities, which have been 
marked with across. There may be a relation between sources Nos. 17, 31, 32 
and 80 and such irregularities at distances of 2, 2, 3 and 3 kiloparsecs respectively. 
The only other source favourably placed for comparison, No. 16, does not coincide 
with any similar irregularity. However, such a relation can only be regarded as 
tentative at this stage. 

There also appears to be some correlation between the 38 Mc/s radiation and 
H1 radiation from sources considerably away from the galactic plane. Lilley (13) 
has given profiles of peak H1 intensity at /= 147° and also at b= — 15°, and found a 
correlation with the extinction of extra galactic nebulae. As may be seen from 
Fig. 4, sources Nos. 18 and 24 also appear on Lilley’s profiles, and both lie in or 
near extensions of the zone of avoidance. For No. 24 the shape at 38 Mc/s is in 
good agreement with the shape of the extension. Since neither gas nor dust can 
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produce the radiation at 38 Mc/s, we are led to the hypothesis that sources Nos. 18 
and 24 represent regions where gas, dust, magnetic fields and relativistic particles 
occur together, as has been conjectured to be the case in spiral arms. 
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Fic. 4.—Comparison of distribution at 38 Mc/s with peak intensity at 1420 Mc/s and with 
mean total extinction. (Hubble, Mt. W. Contr., 21, 139, 1933). 
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Summary 


Continuum emission from the Cygnus-X radio source has been studied at 
seven frequencies in the range 92 Mc/s to 1425 Mc/s. Its spectrum is that 
expected from an H 11 region having a peak emission measure of 6,000 and a 
constant flux density of 47 x 10-4 w/m?/c/s. Evidence is presented which 
indicates that although its structure is irregular it contains no intense regions 
of small angular diameter. 

Measurements taken at the position of Cygnus-X in the 21 cm line of 
neutral hydrogen show an anomalous “ absorption’’ at the frequency character- 
istic of the second spiral arm in Cygnus. The results can be explained 
adequately in terms of the conversion of neutral hydrogen at the position of the 
source to the ionized state. If Cygnus-X is in the second spiral arm at a 
distance of perhaps 6 Kpc as the “ absorption ”’ effect suggests, its mass is 
107 Mo. There is a concentration of O- and B- stars in the Cygnus-X region, 
some of which are contained in a cluster at about 2 Kpc. If these stars were 
within a spiral arm they could produce a large fraction of the ionization 
required. 





Introduction.—All surveys of extra-terrestrial radio emission have indicated an 
excess of emission near the intense source Cygnus-A (19N4A). Bolton and 
Westfold (1) showed that it was an extended source centred approximately on 
Dec. = 40° and R.A.=20"30™, and suggested that it might be the excess of 
galactic background emission associated with the galactic spiral arm extending 
away from the Sun in this direction. However, Piddington and Minnett (2) 
who named it Cygnus-X, concluded from observations of its shape and spectrum 
that it was probably an Hu region. Scheyer and Ryle (3) extended this con- 
clusion by suggesting that the Cygnus-X region was the integrated H 11 emission 
from a spiral arm. However, Baade and Minkowski (4) were unable to make an 
optical identification of the source. 

Recent observations have been made of the radio and visual emission from 
the Cygnus region. They include (2) high resolution surveys covering a wide 
frequency range of the continuum emission from the source ; (4) neutral hydrogen 
absorption measurements which give added information about the source; and 
(c) recent studies of early-type stars lying in this part of the Galaxy. This paper 
describes the results of seven high resolution surveys of the continuum emission 
from Cygnus-X and relates them to 21cm hydrogen line observations of the 
region. A more definite statement about the nature of Cygnus-X is now possible. 

2. Radio surveys of the Cygnus-X region in continuum emission.—Three high 
resolution surveys of the Cynus-X region made at Jodrell Bank are described. 
They were made in the continuum emission of the radio frequency spectrum as 
distinct from the line emission at 1420 Mc/s. 
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(a) 1425 Mc/s.—A total power switched receiver (§) with a bandwidth of 
1:1 Mc/s was used in conjunction with a 30 ft paraboloid having half-power beam- 
widths of 1°-45 and 1°75 in the H and E planes respectively. A facsimile of a 
tracing taken at Dec.=40°-o using an output time constant of 30 seconds is 
shown in Fig. 1. The r.m.s. noise ripple on the chart (+2 °K) is greater than 


aroo™ 2040" 20"20" oo Wao” 
RIGHT ASCENSION 





Fic. 1.—Facsimile of a total power tracing at Dec.=40°-o at 1425 Mc/s. 
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Fic. 2.—Contour maps of the Cygnus region at (a) 1425 Mc/s, (b) 158 Mc/s, (c) 92 Me/s. 


expected from mixer noise alone and is attributed to short term instability in- 


herent in the switching process. The brightness temperature calibration was 
made by the method described in Section 4. 
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In order to construct a map of this region transits were taken every 0°*5 in 
declination in the range Dec. =37° and 45°. At least two transits were obtained 
for each declination setting and the derived mean transits were considered 
accurate to +2°K. ‘The resultant map of the region containing Cygnus-A and 
Cygnus-X is shown in Fig. 2(a). 

(6) 158 Mc/s.—At this frequency a more detailed examination * was made of 
the recordings already taken by Hanbury Brown and Hazard (6) in 1950-51 using 
the 218 ft transit radio telescope. The aerial beamwidth to half power is 2° in 
right ascension and 2°-2 in declination. ‘The Cygnus region was surveyed from 
Dec. =37° to 53°. At least two transits were taken at fixed declinations at 
intervals of less than 1°. The brightness temperatures are those derived in 
Section 4. The recordings were corrected for the side lobes of the intense 
Cygnus-A source between Dec.=41° and Dec.=37° (the greatest tilt of the 
aerial mast at that time). These side lobes of Cygnus-A were symmetrically 
placed in right ascension about the main beam and could also be distinguished 
since they exhibited ionospheric scintillation along with the main beam. No 
attempt was made to derive an absolute zero level for this survey, so the zero was 
taken as the minimum level observed within 50° of the galactic plane. The 
contours of the region are plotted in Fig. 2(d). 

(c) 92 Mc/s.—The Cygnus region was surveyed at 92 Mc/s in collaboration 
with Mr D. Walsh using the 218 ft transit radio telescope. ‘The aerial beam- 
width to half power was 3°:2 in declination and 3°:o0 in right ascension. ‘Two 
transits of the region were taken at fixed declinations at intervals of 2°:5 in the 
range Dec. = 32° to 48°. The brightness temperature scale is that derived in 
Section 4. Again, the zero level is the lowest observed within 50° of the galactic 
plane. ‘The final map is shown in Fig. 2(c) and indicates that there are no side- 
lobes of Cygnus-A greater than 2 per cent. 

3. Observations of the Cygnus region at the hydrogen-line frequency. (a) 
Equipment.—The double comparison spectrometer previously described (§) was 
used in conjunction with a 30ft paraboloid. ‘The spectrometer bandwidth was 
reduced to 18 Kc/s. A time constant of 20 seconds was used when taking spectra 
and this was increased to 1-5 minutes when taking drift curves through the region 
at a particular frequency. The brightness temperature scale was standardized 
against the peak temperature at /=50° and b=o’, taken as 100 °K. 

(b) Spectra.—Itis necessary to obtain the absorption spectrum of the Cygnus-X 
radio source in order to determine its distance. This involves deriving two spectra 
characteristic of Cygnus-X; the first is the spectrum actually observed in the 
direction of the source, 7'(v), and the second is the spectrum of the background 
neutral hydrogen emission, 7',(v), expected if the source were not present. The 
difference between 7',(v) and 7',(v) is the difference or absorption spectrum T>)(v). 
The determination of 7',(v) involves an extrapolation of the spectra taken in the 
vicinity of Cygnus-X. 

Since the source is elongated across the galactic plane, sets of spectra were 
taken along a line parallel to the galactic plane, i.e. across the source. ‘This 
procedure has the added advantage that the hydrogen-line spectra change less 
sharply with change in distance along such a line and interpolation across the 
source is more trustworthy. An attempt was made to obtain 7(v) at Dec. = 40°°5 
and R.A.=20"28™, near the centre of Cygnus-X. Spectra were accordingly 
* Dr Hazard kindly assisted with the interpretation of these results. 
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taken at the following points along a line: Dec.=36°, R.A.=20" 19; 
Dec. = 38°'5, R.A.=20%24™; Dec.=40°%5, R.A.=20" 28m; Dec.=42°'5, 
R.A.=20"32™; and Dec.=45°, R.A.=20"38™. Smoothed curves through 
one such set of spectra are shown in Fig. 3. The spectra show the existence 


Dec. R.A. 
cna. 2 se o.: 
chin 2. ae (x 






a 40.5 20°28" 
----- 42°5 20°32" 
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~40 
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Fic. 3.—A series of hydrogen-line spectra taken across the Cygnus-X source 
(Dec.=40°'5, R.A.=20h 28™) 

of 3 spiral arms in the direction of Cygnus-X. In addition there is a 
frequency drift with declination which is different for each spiral arm. It can be 
seen that the spectrum at the position of the source, T'(v), is clearly reduced at 
the frequency of the first and second spiral arms. The spectrum 7;(v) was 
derived by interpolating, both with respect to position and frequency, between 
spectra on opposite sides of the source. 7(v) is then obtained by subtraction 
of T,(v) from 7;(v). 

Four sets of spectra were obtained at these coordinates with the spectrometer 
and then averaged to give the full line curves of T,(v) and 7p(v) shown in Fig. 4. 
The standard deviation of any point on these curves for 7',(v) and Tp(v) is + 3° K. 
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Fic. 4.—Averaged plots of the expected spectrum Tp») and the difference spectrum T p(v) at the 
position of Cygnus-X. 
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The difference spectrum 7>p(v) shows (a) clear absorption in the nearby spiral 
arm, (6) a fall in temperature in the second spiral arm greater than the observed 
brightness temperature of Cygnus-X, and (c) no absorption greater than the 
limits of error in the outer arm. The Cygnus-X source must accordingly lie 
nearer to the Sun than the outer arm and it exhibits a value of T')(v) for the second 
spiral arm which is anomalous (see Section 5). 

(c) Drift curves.—In order to investigate the difference spectrum of Cygnus-X 
more fully, and in particular to confirm the anomalous value of 7p(v) in the 
second spiral arm, this spectrum was derived from a different set of observations. 
Accordingly, drift curves were taken at the declination of Cygnus-X and at adjacent 
declinations with the spectrometer set at given frequencies. These curves were 
normally taken at Dec. = 36°, 38°-5, 40°-5, 42°°5,45° and 47°. Since the frequency 
of a given spectral component changes with declination the frequency setting at 
each declination was always chosen to be similarly placed relative to the spiral 
arm under investigation, e.g., in order to obtain Tp(v) for Cygnus-X at the mid- 
frequency of the second spiral arm, the drift curves were taken at the mid-frequency 
of the second spiral arm at each declination. One such set of drift curves has been 
smoothed and is shown in Fig. 5(a). The contours of Tp(v) over the Cygnus-X 
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Fic. 5 (a).—A set of drift curves taken at the centre frequency of the second spiral arm in Cygnus. 
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Fic. 5 (b).—Interpolation of Tp(v) (broken line) across the Cygnus-X region at this frequency. 
T(») is shown full-line. 
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region for a given frequency are derived by interpolation across the region along 
lines of galactic longitude. Fig. 5() illustrates the interpolation procedure 
applied to the drift curves of Fig. 5(a) at four different galactic latitudes which 
pass through Dec. = 40°:5 at the right ascensions given; the interpolated bright- 
ness temperature, 7';,(v), is plotted as a broken line and the observed value of 
T(v) is plotted full-line. The derived contour map of 7')(v) at the mid-frequency 
of the second spiral arm is plotted in Fig. 5(c). It can be seen that the contours 
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Fic. 5 (c).—The derived contour map of Tp(v) at this frequency. 


of T>(v) are similar in shape* to those of the continuum emission from Cygnus-X 
plotted in Fig. 2(a) taken with the same aerial at a nearby frequency in the con- 
tinuum. 

This procedure was repeated at 20 frequencies in the neutral hydrogen 
spectrum of the Cygnus-X region. The map of 7>(v) for each frequency was 
similar in shape to that in Fig. 5 (c). The values of 7p(v) and 7,(v) at Dec.=40°°5, 
R.A. = 20" 28™ for each map are plotted as crosses on Fig. 4 in order to give a 
comparison with the results of the previous method. The scatter in the crosses 
gives a measure of the errors involved in interpolation; the difference between 
the results of the two methods lies within the errors of each. Moreover, the large 
value of T)(v) at the frequency of the second spiral arm is confirmed and its value 
at the centre frequency appears to be 34°+2°K. The corresponding value of 
T,(v) is 64 °K. 

(d) The distribution of neutral hydrogen in the spiral arms in Cygnus.—The 
distribution in right ascension of the brightness temperature of neutral hydrogen 
in Cygnus is given by the interpolated plots of 7,(v) already derived for 
Dec.=40°5. The optical depth, 7(v), of the neutral hydrogen at each point in 
the distribution may be obtained from the value of 7,(v) at that point using 


T;(v) = T,(1-e-™) (1) 


where 7, is the kinetic temperature and is assumedt to be 125°K. The distri- 
bution of neutral hydrogen in each spiral arm may be expressed in terms of the 


* The relatively smaller values of Tp (v) at the position of Cygnus-A are expected (see Section 5) 


since Cygnus-A is 5° from the “‘radio”’ galactic plane where the optical thickness of the neutral 
hydrogen is about half that at the position of Cygnus-X. 


+ This assumption may not be justified for the outer spiral arms (see Section 5). 
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number of neutral} hydrogen atoms percm? in the line of sight (=f N,.dl) at 
each point on the distribution. It can be shown (7) that 
{ Ny. dl= 4°62 x 10% 75 Tx (2) 
where 7, is the optical depth at the centre frequency of an arm and 7 is the half 
width of the + spectrum of the spiral arm. 
The distribution of { Nj, .dl across each of the three spiral arms is given in 
Fig. 6 along with the brightness temperature distribution across Cygnus-X at 


@# ARM i 
w © ARMS 263 





12424 


fN,dl IN UNITS OF 107! cwr2 


1425 Mc/s BRIGHTNESS 
TEMPERATURE OF GyGNnus-x (°K) 








° 





RIGHT ASCENSION 


Fic. 6.—The distribution of f Ng. dl across the three spiral arms in Cygnus. The brightness 
distribution across Cygnus-X at 1425 Mc/s is also shown. 


1425 Mc/s. It can be seen that the neutral hydrogen is more broadly distributed 
in right ascension than is the emission from Cygnus-X. In addition, all four 
distributions show the same sharp rise and slow fall with increasing right ascension. 
The shape and position of each distribution is summarized in Table I. 


TaBLe I 
Distance Position of peak f N,,.diat peak Width between 
(kpc) intensity (R.A.) (cm~?) half values 
First spiral arm o-4 200 29M. 17 x 107 11°°5 
Second 5°8 20h 23m 4°5 X 107! 9°°0 
Third 10° 20) 24m 4°3 x 107! 9°'0 
Cygnus-X 20 25m 5°°5 


4. An analysis of the radio continuum measurements in Cygnus. (a) A deter- 
mination of the galactic background spectrum and distribution.—Before an estimate 
can be made of the radio emission from Cygnus-X it is necessary to determine the 
contribution due to the galactic background. Seven surveys of the Cygnus 
region have been made with beamwidths of 3° or less; 92, 158 and 1425 Mc/s 
(Jodrell Bank), 242 Mc/s (8), 400 Mc/s (9), 600 Mc/s (10) and goo Mc/s (12). 








670 R. D. Davies Vol. 117 


These were used to determine the distribution of background emission at 
Dec. = 40°'5, the declination of Cygnus-X. An examination of the seven radio 
maps showed that the background distribution changed only slowly with galactic 
longitude through the position of Cygnus-X and was nearly the same on either 
side of it. The distribution at Dec. = 40°5 was accordingly derived by smooth 
interpolation from either side of the source. The observed radio galactic plane 
(defined as the line of peak intensity) in this region is at +0°*5 galactic latitude 
(defined by the Lund pole). The brightness temperature, T;,’, at each frequency 
is obtained by a comparison with the intensity of the Cygnus-A source and is given 
by the equation 


ny  HayS ~=—P(b-g) 

Va= Gh * P(Cyg-A) (3) 
where yz represents the fraction of the energy falling into the main lobe compared 
with that falling into the area containing the major side lobes (+ 4 beam widths 
from the main lobe), a, is the physical area of the aerial, S is the flux density of 
Cygnus-A derived from a recent study by Adgie and Whitfield, kindly supplied 
in advance of publication, k is Boltzmann’s constant, and P(b-g) and P(Cyg-A) are 


the chart deflections of the background and Cygnus-A respectively. yz and ay are 
given in Table II. 





TasBLe II 
Frequency ao 
(Mc/s) (m*) " 
92 3485 I‘o 
158 3485 0°85 
242 327 1‘o 
400 490 I‘o 
600 94°3 1'0 
goo 44°4 Io 
1425 65°9 I'o 


The resultant brightness temperature distribution is plotted for each frequency 
in Fig. 7. Furthermore, the peak brightness temperature of each distribution 
is plotted against frequency on a log-log scale in Fig. 8; the best-fitting spectral 
law is Ty’ xc y-?24+'16, 

(6) The brightness temperature distribution across Cygnus-X.—The background 
component of galactic radiation derived above can now be removed from the 
observed intensity contours in the region of Cygnus-X. Then, after the removal 
of the effect of Cygnus-A, the emission from Cygnus-X itself remains. However, 
there is some uncertainty in plotting the faint preceding edge of Cygnus-X owing 
to the proximity of the intense Cygnus-A source. Maps of Cygnus-X at the 
seven frequencies are plotted in Fig. 9 along with the aerial polar diagram used at 
each frequency. The shape of Cygnus-X is seen to be similar at all frequencies, 
apart from the effect of aerial smoothing. A comparison of detailed structure 
not readily seen in the maps is provided in Fig. 10, which plots the brightness 
temperature distribution across the source in right ascension at Dec. = 40°°5. 

The source covers 12° in right ascension, with a width of 5°-5 between half 
intensity points. It shows broad maxima at R.A. =20" 24™, 20" 36™ and 20" 48™ 
on the highest resolution surveys. The extension of the source in declination is 
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Fic. 7.—The brightness temperature distribution of the galactic background emission at Dec. = 40°°5 
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at 7 different frequencies. 
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Fic. 8.—The spectrum of the peak of the galactic background radiation in Cygnus. 
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FIG. y.—Brightness temperature maps of Cygnus-X at 7 frequencies. The aerial polar diagrams 
used at each frequency are also shown. 
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Fic. 10.—Brightness temperature distribution across Cygnus-X at Dec.=40°:5 for 7 frequencies. 


less than in right ascension and may be obtained, after correction for aerial 
smoothing from the relation * 

A® = 0? — B? (4) 
where A, O and B are the half power widths of the actual and observed distri- 
butions across the source and the aerial beam respectively. The half power 
width of the resultant declination distribution averages 2°-5 over the frequency 
range; there is some indication that the source width may vary from 3° at the low 
frequencies to 2° at the high frequencies, although such an effect is close to the 
limits of error. 

(c) The spectrum of Cygnus-X.—The flux density of the total Cygnus-X 
complex can be compared directly with the flux density of the Cygnus-A source 





Tasie III 
en Flux of Cygnus-X Standard flux density of Estimated flux density of 
— . Cygnus-A Cygnus-X 
(Mc/s) Flux of Cygnus-A (x 10% w/m?/c/s) (x 10 w/m?/c/s) 
92 0°42 118 49 
158 0°67 87 58 
242 0°97 62 60 
400 1‘2 43 52 
600 1°6 29 47 
goo 1°4 20 30 
1425 26 13 37 


* This relation is correct for Gaussian distributions. 
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by integration of the intensity contours of the two sources. The relative flux 
densities are given in Table III for the seven frequencies. ‘The best estimate of 
the flux density of Cygnus-X is then obtained, applying the standardized flux 
density spectrum for Cygnus-A used above. 

Table III suggests that the flux density of Cygnus-X is some 20 per cent 
greater at low than at high frequencies; however, this difference is scarcely 
significant since uncertainties in specifying the background distribution at the 
lower frequencies could lead to such errors. ‘The mean flux density of the 
Cygnus-X complex is 47 x 10-4 w/m?/c/s. 

A second and independent method of specifying the spectrum of Cygnus-X 
is to measure the brightness temperature of a given region of the source at 
different frequencies. ‘The brightness temperature of the most intense part of 
the source (Dec. = 40°, R.A. =20"25™) at each frequency is given in the second 
column of Table IV. 


TasBLe IV 
Frequency Observed Aerial dilution Derived T (1000 Me/s) 
(Mc/s) temperature re temperature (x) 
(«) (xk) 

92 2200 1°3 2900 34 
158 1400 1°2 1700 42 
242 140 3°0 pod os 
400 210 1'2 250 39 
600 53 I'5 81 29 
goo 23 1°8 43 35 

1425 16 1°2 19 38 
3200 4 1° 4 40 


Included in the table is a value of 4° at 3,200 Mc/s; the value given by Haddock, 
Mayer, and Sloanaker (12) was converted to a brightness temperature and then 
reduced by 1 °K to allow for the background contribution at this frequency (see 
Fig. 8). In order to derive the actual brightness temperature of the source at 
each frequency, the reduction in temperature by the aerial beamwidths is first 
calculated. Column 3 of Table IV gives for each frequency this dilution factor 
which was obtained by a method of graphical integration; Column 4 gives the 
derived brightness temperature; each value is reduced to that expected at 1000 
Mc/s assuming a 7','oc v-* spectrum * and is given as T(1000 Mc/s) in the table. 
T(1000 Mc/s) has a mean value of 35 °K and shows no significant variation with 
frequency. 

(d) The fine-structure of Cygnus-X.—Information relating to the fine structure 
within the Cygnus-X complex was derived from interferometer observations 
taken at Jodrell Bank. A swept-lobe interferometer operating at 158 Mc/s with 
a spacing of 56A showed no feature less than 1° diameter with an intensity greater 
than 15 per cent of the peak value of Cygnus-X +. 

These results were confirmed by a search for scintillations on the g2 and 
158 Mc/s records. At the time of large fluctuations on Cygnus-A no fluctuations 


* If this spectrum arises from free-free emission from ionized hydrogen at 104 °K, the finite 


optical depth at the lowest frequency would cause an apparent reduction in the value of 7(1000 Mc/s) 
by approximately 15 per cent. 


+ Data kindly supplied by Dr H. P. Palmer. 
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greater than 15 per cent of Cygnus-X could be detected. Fluctuations would | 
occur for a source having angular dimensions much less than the angular extent 
of the small dimension of an ionospheric irregularity; this is taken as 0°-4. 


5. Interpretation of results. 


(i) The radio results.—Measurements at the hydrogen-line frequency show 
a fall in the observed brightness temperature at the peak of Cygnus-X of 16°, 34° 
and o°K in the first, second and third spiral arms respectively. These results 
may be compared with the fall in temperature, A7(v), expected when the con- 
tinuous radiation of a source of brightness temperature 7',’ passes through a 
spiral arm of optical depth 7(v) at a frequency v 


AT(v)=Ts(1-e™) = aoe 5) 


The maximum value of AT(v) is Ts’ (which is 16°K) and occurs when the spiral 
arm is optically thick. Accordingly, for the first spiral arm where 7;,(v) is 117° 
at the peak frequency and assuming 7’; to be 125 °K, A7(v) is expected to be 15 °K 
at this frequency, if Cygnus-X lies behind this arm. ‘The observed value is 
consistent with this picture, and moreover the distribution of AT(v) and 7,(v) 
with frequency are similar, as expected. 

On the other hand, the observed value of A7(v) in the second spiral arm is 
twice the value expected on the above picture. This result cannot then be 
interpreted in terms of an absorption effect alone. However, all the available 
data may be explained satisfactorily by supposing that hydrogen in the neutral 
state has been removed from the position of Cygnus-X. Arguments will now be 
presented which suggest strongly that a large fraction of the neutral hydrogen in 
the second spiral arm has been ionized to produce the Cygnus-X complex of H11 
regions. 


(a) The contours of the deficiency of neutral hydrogen emission, T;p(v), and 
consequently of neutral hydrogen itself (if Av(v)<1), are similar to the emission 
contours of the Cygnus-X complex. Hence the two phenomena are related. 
Such an effect is expected if the original neutral hydrogen atoms have been 
ionized and are now radiating by the process of free-free emission. Since the 
neutral hydrogen emission from unit cross-sectional area is proportional to 
f N,.dl and the free-free emission is proportional to f N,*.d/ there may be 
differences in the two contours if N,,, and consequently N,, are not constant over 
the whole extent of the Cygnus-X source. Ny will begin to fall near the outer 
edges of a spiral arm, which, in the case of the second spiral arm in Cygnus, is 
beyond the R.A. limits of the Cygnus-X source. Consequently the contours of 
T(v), (or more accurately Az(v)) are likely;to be similar to the contours of the 
brightness temperature (or more accurately the optical depth) of the free—-free 
emission from the associated source. 


(6) The spectrum of the Cygnus-X complex is that of an H11 region since 
its flux density is constant within the limits of error, having the value 
4°7 x 10° w/m?2/c/s. Moreover, for an H11 region the optical depth, 7'(v), is 
proportional to v-* so the brightness temperature of the source is given by 

Ts'(v)=T(1-€"™) 
=T,.7'(v) for 7'(v) <o°4 say. (6) 








5% 
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For the frequency range used in this investigation 7’(v) < 0-4 if we assume the 

electron temperature, 7, is 10'°K, and consequently 

T,'(v).v?=constant. ~ (7) 
As is seen from the last column of Table IV this quantity has the average value 
3°5 x 10! °K sec for the peak of the Cygnus-X complex. Its constant value 
over the frequency range confirms the H11 region nature of the source. 

(c) The density and depth of the assumed H11 region and of the neutral 
hydrogen deficiency are compatible and lead to acceptable values of these 
quantities. This can be shown as follows with the simplifying assumptions that 
the hydrogen density is the same throughout the volume of Cygnus-X. 

The optical depth of an H11 region with electron density N, and depth / is 


given by (13) 
7'(v) =0°53 T3? v?= N21 x 10° (8) 
where / is measured in parsecs. 
Whence the emission measure 
E.M.=N21= 1-89 x 10°. T!*(T'(v)v*) (9) 

The emission measure for the peak of the Cygnus-X source has the value 
6,600 for T,=10'°K; it will not vary markedly with T,. This relation between 
N, and / will be compared with a similar relation for the neutral hydrogen. 

The line integral of the neutral hydrogen removed from the second spiral arm 
at the peak of Cygnus-X can be estimated from the expected brightness tempera- 
ture, T,(v), of 64°K and the observed brightness temperature, T,(v), of 30 °K. 
Now, 


T,(v) = Tx(1 -€-"™) (10) 
where T, is the kinetic temperature of the neutral hydrogen. Then, if 7,(v) 


is reduced to 7,(v) by the removal of neutral hydrogen of optical depth Ar(v), it 
can be shown that 


Tx —Ts(v) 
Tx-T3(v) (1) 
Furthermore the line integral of the neutral hydrogen which has been removed 
will be given by (see equation (2)), 

Ng -!=Ar(v). Tn x 1°54 x 107° (12) 


Ar(v) =In 


where / is measured in parsecs; 7 is 10km/sec for T:p(v) in the second spiral arm. 
Accordingly the estimate of Nj, ./ is dependent upori the kinetic temperature 
of neutral hydrogen in the Cygnus-X region. The value of 7’, found for the inner 
parts of the Galaxy is 125°K (14) but there is some evidence of strong local 
variations (15). In the case of the two outer spiral arms in Cygnus the brightness 
temperature nowhere reaches 80 °K; this represents a possible lower limit of Tx. 
Table V illustrates the variation of Ny ./ with 7’, in the range 80° to 125° K. 


TABLE V 
Tx Nu. N,=Ng l 
(°x) Ar (v) (cm? parsecs) (cm-) (parsecs) 
80 1°14 1410 4°6 310 
100 0°66 1020 6°5 160 


125 0°44 850 78 110 
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Now if the ionized hydrogen is derived from the neutral hydrogen removed 
from the second spiral then N, and Ny may be equated and equations (g) and (12) 
will yield the estimates of N,, Ny and / given in Table V for different values of Ty. 
Thus the suggestion that the Hi region can be identified with the neutral 
hydrogen deficiency leads to no contradictions, and it will be shown that it results 
in quite acceptable values of N,, Ny and I. 

The depth of the region / can be compared with the transverse dimensions 
of the Cygnus-X complex at the distance of the second spiral arm (6kpc); the 
half-intensity dimensions are 270 parsecs by 600 parsecs. The neutral hydrogen 
density is perhaps a factor of two greater than is normally found in the centre of 
spiral arms (14), but it is not excessive. Moreover the number of neutral hydro- 
gen atoms per cm? in a line through Cygnus-X is about 3 x 107! compared with 
2 x 10” for the Cygnus region and 1 — 2 x 10” elsewhere in the galactic plane (14). 
The mass of the Cygnus-X complex ranges from 8—11 x 10® solar masses for 
the range of 7, in Table V. It would correspond to 1-9 — 2-7 x 108 solar masses 
if at a reduced distance of 3 kpc and may be compared with the galactic mass of 
1-6 x 10! solar masses (16). aa 

(d) The brightness distribution across the Cygnus-X complex is that expected 
from an H 1! region rather than an aggregate of “point’’ sources. The results 
of Section 4(d) showed that only a small amount of the radiation of the Cygnus-X 
complex could arise from intense regions of the order of 10’—20’ of arc in diameter. 
The source would appear to be irregular in shape and to show only smooth 
variations in intensity across its area. 

(e) Further evidence for the H 11 nature of Cygnus-X comes from the 22 Mc/s 
survey by Burke (17) which shows that “this area of the sky is less bright than a 
suitable average over the remaining sky”. Such an effect is expected if the 
kinetic temperature of the H11 region is less than the brightness temperature of 
the emission from behind the source. An extrapolation of the spectrum in Fig. 8 
indicates that the brightness temperature of the total background component at 
22 Mc/s will be of the order of 25 000 °K compared with a likely kinetic tempera- 
ture of the Cygnus-X H11 region of 10000°K. ‘Thus the Cygnus-X region 
would absorb the high temperature background radiation and appear cooler than 
its surroundings at 22 Mc/s. 

(ii) Optical studies of the Cygnus-X region.—This region of Cygnus is rich 
with H 11 regions interspaced with dust clouds. However a direct identification 
between H 11 regions and parts of the Cygnus-X complex will be difficult because 
of the heavy obscuration amounting in places to 7 magnitudes (18). Nevertheless 
recent studies (19, 20) of early type stars in Cygnus reveal a likely origin of the 
intense ionization required. A catalogue of the types of O-star given by Hiltner 
and Johnson (20) within the boundary of Cygnus-X is given in Table VI. These 
stars lie beyond the Cygnus Rift at 700 parsecs (21). The value of the quantity 


TasB.e VI 
Ss Number of stars fN.?.aV ~ fN,?.dV 
tar tyPe within Cygnus-X complex for one star for all stars 
05 5 3°4 X10% 17°0 X 10% 
06 6 16 x10% 9°6 x 10% 
07 8 0°81 x 10% 6°5 x 10% 
08 8 0°35 X 10% 2°8 x 10% 
o9 13 o'12 X 10% 1°6X 10% 


47 
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J N,2.dV integrated over the volume V ionized by each star type was estimated 
from the data of Strémgren (21) and is given in the third column. 

The total value of { N,2.dV for all these stars is 3-7 x 10° cm™*. An estimate 
of this quantity can be derived from the flux density measurements of Cygnus-X in 
the continuum. It can be shown that for free-free emission from an H 11 region 


ce : 1/2 ‘ 

.S§ I 
o3s4h ih 
where c= velocity of light 

T,=electron temperature, taken as 10*°K 
k=Boltzmann’s Constant 

D = distance of the source 

S=flux density of the source in c.g.s. units. 


For the source at 3kpc and 6kpce the values of fN,?.dV are 7-1 x 10% and 
2°8 x 10%cm-* respectively. Hence the stars already classified can produce a 
substantial fraction of the ionization required. Most of these stars are found 
within the preceding half of Cygnus-X; the following half is much more heavily 
obscured. i 

An early-type star cluster, 4° in diameter and called VI Cygni, has been 
discovered within the boundary of Cygnus-X at Dec.=41° and R.A. =20"29™ 
(18, 23,24.) This cluster, provisionally placed at 2 kpc, shows absorption up to 7 
magnitudes and colour excesses up to 3 on the B—V index (25). 20 members have 
been classified as O-type and another 65 fainter members are shown, in low dis- 
persion spectral surveys (26), to have similar spectra. The 14 stars classified 
would give an emission measure (27) of 2800* if N, were taken as 4 (see Table V). 
The radio surveys indeed show a peak in this position and the emission measure 
appears to be 4000. This difference in emission could readily be accounted for by 
the unclassified stars of the cluster. 

6. Conclusion.—The above analysis shows that Cygnus-X is not the integrated 
effect of radio sources following the Cygnus spiral arms as defined by the neutral 
hydrogen measurements. The shape, extent and spectrum are inconsistent with 
such an hypothesis. The continuum measurements indicate that the source is an 
H 11 complex ; this suggestion is confirmed by a marked absence of neutral hydro- 
gen within the contours of the source. ‘This neutral hydrogen effect occurs in the 
second spiral arm at an apparent distance of 6kpc. Optical observations show a 
concentration of O stars in the section of the Galaxy between /= 40° and 48° and 
the stars already catalogued could provide a large amount of the observed H 1 
emission if they were within the hydrogen of a spiral arm. ‘The VI Cygni O-star 
cluster, }° in diameter, is most likely responsible for the observed emission 
maximum at Dec. = 41° and R.A.=20"29™. However, the measured distances of 
the radio and optical phenomena do not agree closely ; this does not invalidate the 
above results since neither measurement is capable of great accuracy in this 
region. It is concluded that the Cygnus-X source can probably be identified 
with the H11 regions likely to be produced by the observed concentration of 
O-type stars. 
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: * By virtue of the heavy absorption this would be reduced in the optical range to about 30 which 
is on the limits of detection, even in a clear part of the sky (28). 
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THE SPECTRA OF RADIO STARS 
G. R. Whitfield 


(Communicated by M. Ryle) 
(Received 1957 September 23)* 


Summary 


By considering all available measurements of the flux density of Cassiopeia 
A it is concluded that its spectrum obeys a simple power law. Using the fluxes 
given by this law and the relative intensities of other stars given in all available 
surveys, the spectra of 31 radio stars are derived. 

In the majority of cases it is only possible to assume a simple power law, 

flux density oc (frequency)?, 

and to obtain the spectral index, x, for each source. It is found that the mean 
spectrum is steeper for extra-galactic sources than for galactic ones, while that 
of the unidentified sources is steeper still. 





1. Introduction.—The determination of the spectra of radio stars is of 
considerable importance, both in relation to theories of the origin of the 
emission (1, 2, 3, 4) and to the possibility of distinguishing between different 
classes of source. The establishment of accurate absolute measurements of 
flux density over a range of frequencies has proved difficult, as can be seen by 
the spread of the values obtained for the four most intense sources by different 
observers (Figs. 1 and 2). 

The two principal uncertainties in such measurements are the absolute gain 
of the aerial system and the calibration of the receiver sensitivity. It is generally 
only possible to compute with accuracy the gain of small arrays of dipoles, or 
horn aerials; the former have been used for absolute measurements at low 
frequencies (5, 6), and the latter at higher frequencies (7). 

Earlier measurements of the absolute sensitivity of the receiver have been 
based on the shot noise in a saturated diode, but measurements have shown (8) 
that even at low frequencies serious errors can arise due to impedance trans- 
formation in the electrode structure unless special coaxial diodes are used. 

Recent observations (5, 6, 7) have used the thermal noise from a resistance 
as a standard source, and these measurements are believed to be considerably 
more accurate. 

Although absolute measurements are difficult, it is relatively easy to compare 
the fluxes of two sources, provided that the aerial can be directed at each source 
in turn without disturbing effects due to ground reflections, and that the linearity 
of the receiver response can be established. If, therefore, the spectrum of one 
of the most intense sources were known from accurate absolute measurements, 
it would be possible to use observations of the relative intensities of other sources, 
made with uncalibrated systems, to establish the spectra of a large number of 
sources. ‘This method was first suggested to me by Adgie. 

In the present paper, the evidence provided by many different observers 
has been used in this way to derive the intensities of a number of sources relative 

* Received in original form 1957 July 22. 
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to that of Cassiopeia A. By giving special weight to those observations which 
have used thermal calibrations and detailed computation of the aerial gain it 
is concluded that the spectrum of this intense source does not differ significantly 
from a straight line except at the lowest frequencies. The intensities at other 
frequencies can then be interpolated and used to derive the absolute fluxes of the 
other sources. 

_2. The spectrum of the intense source in Cassiopeia.—All known observations 
(with the exception of a few early measurements) of the flux density of 
Cassiopeia A are plotted on a logarithmic scale in Fig. 1; the slope of such a 
curve will be referred to as the spectral index x. The errors indicated represent 
the observer’s estimate of the accuracy of his observation, where this is given; 
but a few observers do not quote errors, and that shown is the author’s estimate. 
Those measurements that included accurate determination of the aerial gain 
and thermal calibration of the receiver are indicated by the large circles. 
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1000 10,000 
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Fic. 1.—The spectrum of the intense source in Cassiopeia. Observations which included thermal 
calibrations and accurate calculations of the aerial gain are indicated by the larger circles. 


The large decrease in flux density at low frequencies is almost certainly due 
to absorption in H 1 regions (9), but it can be seen that in the range 30-10 000 Mc/s 
the most accurate observations lie on a straight line, and can be represented by 
the relation 

flux density « (frequency)*, 


where the spectral index, x, has a value of —o-8o. 

In 1955 it was suggested (10, 11) that the relatively small values of flux density 
found in the 400-500 Mc/s region represented a true departure from this simple 
relationship, and further measurements in this region are needed. It is felt 
however that the present observations do not support this suggestion, as the 
adoption of a linear relationship for Cassiopeia A also leads to simple spectral 
curves for the sources in Cygnus, Taurus and Virgo (Fig. 3). If there were a 
reduction in the 400-500 Mc/s region it would imply that a similar irregularity 
also occurred at these frequencies in the spectra of these and other sources. 
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Fic. 2.—The spectra of the intense sources in Cygnus, Virgo and Taurus. 
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Fic. 3.—The spectra of the intense sources in Cygnus, Virgo and Taurus. 
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In this paper it will therefore be assumed that, apart from the low-frequency 
reduction which is consistent with absorption in an H11 region, the spectrum of 
Cassiopeia A is a straight line. 

TasBie | 


Details of surveys used 


Symbol Observer and reference Aerial Frequency Correction 
system Mc/s factor used 
A Adgie and Smith, 1956 (5). 20-30 A 38 0°07 
interferometers 81°5 0°97 
210 0°85 
500 1‘20 
B Burke and Franklin, 1957 (12). 1}° pencil beam 22 ° 
BL Blythe (private communication). 3° pencil beam = 38 0°97 
BSS Bolton, Stanley and Slee, 1954 (13). 50A 100 0°99 
interferometer 
BWSS Bolton, Westfold, Stanley and Slee, 1954 Multispacing 100 o°4 
(14). interferometer 
Cc Conway, 1957 (15). 470A 500 
interferometer 
D Denisse, Le Roux and Steinberg, 1955 24° pencil beam goo o'91 


(16). Revised values from personal 
communication of Denisse to Seeger.t 


G Grebenkemper, McClain and Hagen, 8° pencil beam 1935 0°56 
1957 (17). 
HB Hanbury Brown and Hazard, 1953 (18). 2° pencil beam 158 1°40 
Hanbury Brown, Palmer and Thompson, Multispacing 158 1°40 
1954 (19). interferometer 
H Hey and Hughes, 1954 (20). 20A 22°6 0°53 
interferometer 
HMH Hagen, McClain and Hepburn, 1954 (21). 1° pencil beam 1420 0°83 
HMS Haddock, Mayer and Sloanaker, 1954 4° pencil beam 3200 0°79 
(22, 23). 
HM Haddock and McCullough, from personal 10’ pencil beam 9500 1°O 
communication of Haddock to Seeger.t 
J Jennison, (9 and private communication 11° 24° beam = 125 1°24 
from Lovell). 
K Kraus, Ko and Matt, 1954 (24). 17° X1°'2° beam 250 1°58 
L Lamden and Lovell, 1956 (9). 10°—15' 16°5 >3°5 
pencil beam 19 1°65 
22°6 1°09 
30 I'll 
M Mills, 1952 (25). go A 101 o'9! 
& interferometer 
MS Mills and Slee, 1957 (26). 49’ pencil beam = 85°5 0°75 
(Mills Cross) 
MSS McGee, Slee and Stanley, 1955 (27). 2° pencil beam 400 5°57 
OB O’Brien, 1954 (28). 50A 38 1°45 
interferometer 
P Piddington and Trent, 1956 (29). 3° pencil beam 600 0°76 
PM Piddington and Minnett, 1951 (30) and 3° pencil beam 1210 0°67 
1952 (31). 
R Razin and Pletchkov, 1957 (32). Not known 300 4°8 
32’ pencil beam 9400 1°06 


* These observations give only relative intensities. 
+ Note added in proof.—These measurements have now been published (8). 
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TaBLe I (continued) 
Details of surveys used 


Symbol Observer and Reference Aerial Frequency Correction 
system Mc/s __ factor used 
RB Rishbeth, 1957 (33). 49° x 83’ 85°5 0°75 
pencil beam (part 
of Mills Cross) 
Ss Seeger, 1956 (7). 20A 400 1°09 
interferometer 
Seeger, Westerhout and van de Hulst, 2° pencil beam 400 1'09 
1956 (34). 
SC Seeger, Westerhout and Conway 2° pencil beam 400 109 
(private communication). 
SH Shain and Higgins, 1954 (35). 17 pencil beam 18 0°29 
SS Stanley and Slee, 1950, (36). 25-100 A 40 0°67 
interferometers 60 0°74 
85 0°96 
100 0°94 
160 1°54 
WwW Whitfield, (preliminary results, 60 A 38 * 
unpublished). interferometer 
WE Wells, 1956 (6). 10-20 A 12°5 0'50 
interferometers  15°5 0°73 
18°5 * 
26°5 0°43 
WH Westerhout, 1956 (37). 1°"9 X2°'8 1360 0°74 
pencil beam 
2C Shakeshaft, Ryle, Baldwin, Elsmore and 157A and 14A_ 81°55 0°97 
Thomson, 1955 (38). interferometers 
3C Shakeshaft, Edge and Archer. 307A and 27A_ 159°5 Me 
(Preliminary results ; private interferometers 


communication.) 


* These observations give only relative intensities. 


3- The observational data.—The sources examined were selected originally 
from the second Cambridge survey (38) and all sources in this list having a flux 
density greater than 100 x 10-** wm~*(c/s)-! were examined. A number of other 
intense sources in the southern hemisphere were also included. Since the 
sources were all selected from low-frequency surveys thermal sources have been 
largely excluded. 

Of the 54 sources selected it has been possible to derive reasonably good 
spectra for 32 sources. Fifteen have not been observed over a sufficiently wide 
range of frequencies and 7 require more detailed examination. The latter 
class includes Cygnus X, the source at the galactic centre, and other sources 
which appear to have large angular diameter and fine structure and which are 
partially resolved in some of the observations. 

All known surveys have been used, with the exception of some early ones 
which have been superseded by more recent work. The complete list of surveys 
used is given in Table I, which, in addition to published work, includes 
unpublished Cambridge work at 38 Mc/s by O’Brien, Blythe and the present 
author: preliminary observations at 159 Mc/s by Shakeshaft, Edge and Archer 
with the large Cambridge interferometer (39): recent observations at 400 Mc/s 
by Seeger, Westerhout and Conway (private communication) ; and observations 
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at 9500 Mc/s by Haddock and McCullough (given in personal communication 
from Haddock to Seeger). 

In each survey a correction factor is derived which brings the measured flux 
of Cassiopeia A to the assumed straight line of Fig. 1, and this factor is used 
to correct the measured flux of each other source. In some of the surveys 
Cassiopeia A was not included, and the correction was detived from the 
observation of some other intense source, usually Cygnus A. 
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Spectra from ratios with Cassiopeia A. 


4. Results.—The spectra of 31 sources are shown in Figs. 3-8. Fig. 3 shows 
the spectra of the sources in Cygnus, Virgo and Taurus. Figs. 4-8 show the 
spectra of the 28 weaker sources, arranged in order of decreasing spectral index. 
The curves in Figs. 4-8 have been displaced vertically, so the units of flux density 
are arbitrary. Where two or more points are too close together for clarity they 
have been averaged ; such points are labelled with the symbols of all the observers 
whose values are included in them. 

It can be seen that the observations of the more intense sources fit closely to 
smooth curves; this good agreement is felt to justify the scaling method adopted 
since it indicates that the discrepancies noted in the earlier discussion can be 
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explained in terms of simple instrumental errors. The precise shape of all the 
curves depends on that assumed for Cassiopeia A, which is based on a small 
number of measurements using thermal calibrations; it is clearly important 
to repeat these observations and to make similar measurements at different 
frequencies, but the correction factors assumed seem to be a good first 
approximation. 

It is interesting to note that, while the spectra of Taurus A and Virgo A can 
be well represented by straight lines, that of Cygnus A, besides the low-frequency 
cut-off ascribed to H 11 absorption, shows a marked change of spectral index near 
350 Mc/s. 

In a number of the other sources there is evidence from 38 Mc/s observations 
of a low-frequency cut-off, but this does not appear to occur predominantly for 
sources at low galactic latitudes; nor would such a large effect be expected at 
38 Mc/s except for sources very near the galactic plane. Even Cygnus A and 
Cassiopeia A, at galactic latitudes of 5° and — 2°, show no appreciable absorption 
above 30 Mc/s. 
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Fic. 9.—Histograms showing the distributions of spectral index for different classes of source. 


For most of the sources it is only possible to assign a mean slope of an assumed 
linear spectrum. The spectral indices of all the sources can, however, be 
determined with reasonable accuracy and these are listed in Table II. 

5. The spectral indices of different classes of source.—A number of the sources 
have been identified with optical objects and it is interesting to compare the 
spectral indices of different classes of object. Initially the identified sources can 
be divided into galactic nebulosities, including supernovae remnants, and abnormal 
extra-galactic nebulae. The Andromeda nebula, the only normal nebula included, 
has a spectral index of —0-6+0-2, in good agreement with the value of —o°5 
found for the Galaxy by Adgie and Smith (5). 

Histograms are plotted (Fig. 9) of the spectral indices of these two classes 
of object and also of the unidentified sources. It can be seen that there is a 








Vol. 117 


G. R. Whitfield 


688 





(bx) .c~ pure (Lb) £ x /S.9 


(L¥) +x £ 
«1 > 
«I 


«5 


+(€S) aD I 
(v1) ool 


(S¥) ,z71> 
(1S) ,vz 


(6b ‘gh ‘Lh) S.£x 5.5 
“o> 

(61) oF: I 

#1 > 


*&~ 


(g€) ,oz pur ,(St) ,F.z 


(zv) ,ob1 
#(M-a) <> 
#2 
S/O], COI = AvaU 


*($$) Atuo 
20INOS J9}JBWIBIP [[BUWIS JO SI UNIQDIdG «= *-gzZIS DON 1.0 99:0— o$L1 
(WS) JajaureIp 103 pajde1109 
ore saiisusjuy “(08 ‘ov) Lg JA ‘Axeye8 yewsouqy ZzI.o $l.o— ogli 
9-0 6.1- of 
9-0 g-1— ol 
S.o Z.2— og 
p'Axeye3 ajqnog 1.0 O.1- oot 
*(zS) Aqisojnqou onoeey) = z.0 g-0- ogS 
}-Aazernoad 
SNOIAQO OU YM 3nq ‘UuOTZISOd jOBxa UT Axeles AUIeY z.0 $g.0 $6 
*(1S) 1939WeIP 10J pd}da1I09 are SaIyIsuajUyT «“ Ebb DT z.0 SS.o— oo$ 
‘(oS 
‘ov) Lgor jo eAousadns jo yuBUUAI ‘eINgeNY quid) $0.0 Qz-o— ool! 
z.0 Z.1— S11 
(61) Aqtsopngou s0RTes) ~=—s«£.0 o.1— oti 
z.o Sg.o— oOgz 
$.0 g-o— $6 
*A[UO 901n0s 
JajJ@UIvIp |[eUIs Jo UNIQIadG ‘(gv ‘ov) SLZTIDQON = z-0 6.0— oft 
“(ob ‘Eb ‘zb) fl ‘ejnqou epowoipuy z.o 9-0—- 061 
» A®AINS DE 9y3 UO paafosez 
Ajjeizsed st pue (gz) S—N popuaixa st s01nos styy, = S.0 g-o— oL 
*(1¥ ‘ob) ZZS1 ‘eaoutadng scyeig oysAy, =z-0 L.o— oSz 
+ 1-(8/9)g—U1M,_O1 > 
SyIBUlOY 10119 8/DJA[ COT 38 
UMNUITXey, «=- XOpuy Ayisusq 
paeunsy emoedg xn] J 


saoiput poajoads fo Aamuungy 
I] a1av1, 


JajouBIq] sewrxoiddy 


Vv wo 
Ivol Dz 
VIIA 
ZIOI Dz 
SSg az 
fg oz 
9°8 De 
Vv PAH 


zl Dz 
Les gz 
Igv az 
V nel 
ob az 
Ltr az 
bor az 
6LE az 


g6z az 
9§ Dz 


o$ az 
bf az 


‘ON Oz 


pue oweyy 


vrsti 


VINZI 


Wig6o 
Wrsgo 


VNB 
VzN9° 


WZNSo 
Vrnso 
VENto 


vrNto 
vPNoo 


V9NOO 


‘ON 
‘A'v'l 





Ae AEE A TO 








The spectra of radio stars 


No. 6, 1957 


*UONBOTUNUIUIOD a38ALIg 


*(6S ‘zS) Aqyisojnqou onoepes 


Tv eredoissey 03 sejruns ‘AjIsopngqou Arejuaueyy 


*s/9J| OSE aaoge adojs sty} sey umy90dg 
*s/OJAI OSE 03 OF wiorZ adoys sty} sey 
*(gS ‘ZS) eurxeu om} yuM ‘arnjonI3s 

sey a0mn0g ‘(zS) sarees Burpryjoa 


umi330dg 
xo] duos 


*(g£) bogi ‘eaousedns s_19]da yy 

‘[]]2” AI9A 99188 Jou Op suOTjIsod 9y} 
inq f‘Axeyes rernsed & YIM UONvOyHUApI s[qQISsOg 
*(g£) sorxeyed moj jo 193sN]9 ‘gg1g DON 

§-Ajaytun sity} oyeur (, gf) suomeasasqo yu9001 


ing ‘(9S) ast] “Qy'] ut Lz¢S DON YUM peygnuapy 


syseUlsy 


*UOTJBOTUNUTWIOD SIVA 
“UOHBITUNUTWIOD 9JBALI 

“UOIBOTUNUIWIOD 33BALI 
‘sayory ‘[ °H ‘A ‘S pue e3pq “qd 


$0.0 0g.0— 00061 (gS ‘ZS) ,¥ 
S.o L.2- og (gf) ,4> 
9-0 g-1-—> 06 e I> 
b.0 S.o— ors (61) ,S.1 
£.0 I.I— Str * I~ 
1.0 z0.1— } 
go.0 99-0— } oogit (LS) ,of x ,z 
$.0 o.1— S11 Pa 2 
$.0 Sz.1— og i> 
Z.0 Sg.0— oot (gf) ,9 
£.0 6.0— 061 (gf) ,o7<— 
z.0 €.1— ogs #,59-2 
S.o g.I— OoI «I~ 
£.0 Z.1— 06 *,2:I> 
Z.0 g-0— oor (S¥) ,z1> 
¥ 1~(8/9)g—-UIMy, OF x 
IOI S/DJAI COT 38 
UMUTTXBIA, «= XOpuUy Aysueq s/XJAI COL =Iw8U 
powumsy eroedg xn] J JojourViq] ayeunxoiddy pue swe 


sastpur poujzads fo Aanuuuny 
(panusjuor) {] a1av], 


"ageysayeys *y ‘f 

‘PSMOUTA, “Y 
“"ssaTyMaq] “MA “CT 
‘ageysoyeys “yf 





o6gI az 

V 88) 
oLgI az 
SLL1 az 
SzL1 gz 
9891 Dz 


ELg1 Dz 

V BAD 
6981 Dz 
SgbI az 
ELv1 gz 
Shi Dz 
ztri az 

VPH 
ZOvI Dz 
6Sz1I az 


SLI1 92 


‘ON Oz 


SEI SOR RETIN 





VSNEz 


Wrnoér 


VON9! 
VINO! 


WSN?I 


‘ON 
‘NVI 





690 G. R. Whitfield Vol. 117 


difference between the galactic sources (mean x= —0°74) and the extra-galactic 
sources (mean x= — 1°05). The unidentified sources have a mean spectral index 
of —1-21 which is steeper than that of either class of identified source. These 
differences are small and may not be significant, but if real they are of considerable 
importance. 

I would like to thank all those who have provided unpublished information, 
and my colleagues for their advice and assistance. I am indebted to Trinity 
College, Cambridge, for a Coutts Trotter Studentship. 


Cavendish Laboratory, 
Cambridge: 
1957 July 20. 
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OBSERVATIONS OF DISCRETE RADIO-SOURCES 
AT A FREQUENCY OF 500 MC/S 


R. G. Conway 
(Received 1957 September 26) 


Summary 


New observations have been made of eleven discrete radio-sources, to 
determine their flux densities at 500 Mc/s relative to that of Cygnus A. A 
new technique is used that increases the available sensitivity by enabling 
very long integration times to be used. The flux densities thus obtained have 
been combined with the flux densities at lower frequencies by Whitfield (12) 
to derive a spectral index for each source. The mean spectral index of the 
eleven sources investigated here, none of which is yet identified, is significantly 
different from the spectral indices of those identified with galactic objects. 





1. Introduction.—A recent survey at 81-5 Mc/s by Shakeshaft, Ryle, Baldwin, 
Elsmore and Thomson (1) has revealed a large class of radio-sources the diameters 
of which are less than the limit of measurement, i.e., about 10 minutes of arc. 
These sources, which in the present paper will be termed “ radio-stars ”’, appear to 
be distributed isotropically, and only a few of them have been identified with 
objects observable at optical wave-lengths. It is not known whether or not the 
remainder are similar in type to any of the identified sources. 

One way of establishing such a relationship would be to compare the radio- 
frequency spectrum of a number of radio-stars with the spectrum of identified 
sources. ‘This comparison has not hitherto been possible, because most published 
observations of radio-stars have been limited to a single octave in frequency, 
80-160'Mc/s. Apart from the four brightest sources, which are all identified, 
sources observed at higher frequencies by Hagen, McClain and Hepburn (2), 
Haddock, Mayer and Sloanaker (3) and Piddington and Trent (4) belong to a 
different class, sometimes referred to as Class I, which have appreciable angular 
diameters, and show a concentration towards the galactic plane. A number of 
them are known to be galactic nebulae emitting thermal radiation. 

The observations described below were made at a frequency of 500 Mc/s with 
an interferometer of high resolving power. The observations were thus restricted 
to the small-diameter radio-stars, without confusing effects from the more intense 
extended sources. Eleven of the most intense small-diameter sources were 
selected from the list of Shakeshaft et al. (1); their positions are given in Table I. 
Observations were made of the relative intensities of the sources at 500 Mc/s, with 
Cygnus A as acomparison source. The flux densities found in this way are com- 
bined with the flux densities at lower frequencies to give a spectral index for each 
source that may be compared with the spectral indices of identified sources. 
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TABLE I 


List of sources 


Identification R.A. Declination Flux density at 81-5 Mc/s 

Number hm s . 10-** MKS 

2C 50 00 35 57 —O1 59 145 
379 04 15 12 +38 o1 ; 14I 
404 04 33 57 +29 32 400 [AU 04N3A 
440 05 O1 17 +37 58 175 
724 08 10 12 +48 35 124 o8N4A 
806 09 15 47. —I2 02 570 o9S1A (a) 
1175 14 09 49 +52 37 130 14N5A 
1259 15 02 48 +26 02 145 
1432 16 48 41 +05 06 goo 16NoA (6) 
1775 21 19 52 +2017 124 
1843 22 12 o1 —17 16 202 
CygA 19 57 45 4° 35 14 200 19N4A 


The data in this table are from the survey by Shakeshaft et a/. (1). 


Notes : (a) Hydra A, possibly identified with a faint double galaxy. 
(6) Hercules A, possibly identified with a faint galaxy. 


2. The observations.—The aerials used in this experiment were two 7:5 m 
paraboloidal mirrors, each with a beam width of about 6 degrees. The aerials 
were connected through preamplifiers to a phase-switching receiver, to form an 
interferometer with a spacing of 470 wave-lengths. Such a system gives an output 
that varies sinusoidally with time, with a period of about 30 seconds. A specially 
designed filter, modelled on one by McAdam (an account of whose work is in 
preparation), was placed after the receiver to detect the signals with greater 
signal/noise ratio than would be possible with a normal RC integrator. 

The filter consisted essentially of a reversing switch, followed by an RC inte- 
grator with a time-constant of one minute. ‘The switch could be set by a gear 
train to reverse with a frequency close to the calculated frequency of the star’s 
trace, i.e., about 2 c/min. If the switch frequency were exactly equal to the 
frequency of the unfiltered trace, the output would be unidirectional and propor- 
tional to the amplitude of the unfiltered trace. In practice, to avoid error due to a 
possible zero wander, the switch frequency was made to differ slightly from the 
signal frequency. The output is then quasi-sinusoidal, with a frequency equal to 
the difference frequency. 

The unfiltered signal may be written as 


A(t).cos2mfyt 


where A (t) is a time-varying function describing the reception pattern of each of 
the two aerials, and f, is the frequency of the signal. ‘The switching function is 
cos 2nf,t+higher harmonic terms. 
The higher terms arise because the switching function is a “‘ square-wave”. ‘The 
output signal is thus 
| A(t). cos2m(f,—fo).t+ higher harmonic terms. 


The subsequent integrator reduces the unwanted higher frequency terms and the 
noise level of the trace, but allows the low frequency signal to pass. 
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The output from this filter contains only half the available information. The 
remainder is obtained from a similar switching filter, synchronised with the first 
in phase-quadrature. ‘The smoothed output from the second filter may be written 
as 


4 A(t). sin2m(f,—fo)-t. 
The two outputs, “ sin ” and “ cos’, were recorded separately. 
Fig. 1 illustrates the increase in signal/noise ratio that can be achieved with 


the filter. The traces are records obtained from Virgo A, (a) before filtering, (b) 
from the “ sin ” output, and (c) from the “ cos ”’ output. 





Fic. 1.—Records of Virgo A before and after passing through the switching filter: 
(a) Unfiltered trace 
(b) “‘ Sin’’ trace 
(c) ‘‘ Cos’’ trace. 
The dotted line indicates the time of transit of the source across the meridian plane of the 
interferometer. The vertical markers on traces (b) and (c) are 30 minutes apart. 


An accurate estimate of the R.A. of each source may be made from the phase 
of the output traces. The results obtained in this way are being incorporated by 
Mr B. Elsmore in a paper at present in preparation. 

The filter is used for a similar purpose as the “ rotating-lobe interferometer ” 
described by Hanbury Brown, Palmer and Thompson (5). The signal/noise 
ratios obtainable with the two instruments are the same, but the present instrument 
is preferable when, as in the work desctibed here, the phase of the original signal is 
required. ‘ 

Each radio-star was observed daily near transit for a period of about 30 minutes. 
The observations were in most cases continued for about 14 days. The switch was 
started manually before each observing period and at the same sidereal time each 
day, in order that the phase of the difference frequency (f,—f)) should remain 
constant from day to day. The daily records could therefore be summed to in- 
crease the signal/noise ratio. ‘To do this the traces were sampled at intervals of one 
minute throughout the observing period and were summed numerically. The 
Fourier component of frequency (f, —f,), the amplitude of which is proportional to 
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the flux density of the radio-star, was extracted from the summed records. ‘The 
scale of flux density was determined from records of Cygnus A, which were treated 
in exactly the same way. 

In order to correct for any non-linearity of the system, the receiver and filter 
were calibrated together by replacing the aerials by a diode noise source, which was 
switched on and off, at a frequency of 2c/min, to simulate a radio-star. The 
strongest signal injected gave an output deflection equal to that of Cygnus A, and 
the weakest signal was 1/1ooth of this value. The calibration indicated a slight 
non-linearity, which has been corrected by increasing the apparent amplitude of 
Cygnus A by 1°30 + 0°30. 

The method of calibration allows for the variation of gain due to the extra noise 
power from the sources themselves, over and above the power from the pre- 
amplifiers (6). ‘The variation is in any case small, being less than 2 per cent. 

A correction must be applied for the reduction in amplitude due to the 1-minute 
integrator. ‘This correction depends upon the frequency (f,—f)); in every case it 
was less than 20 per cent. A second correction must be included to allow for the 
fact that the amplitude is an average over the observing period, and is less than the 
central value by a factor which depends upon the directivity of the aerials and the 
exact length of the observing period. The required correction factors were 
calculated from the known directivity pattern, A (t); they differed from one another 
by not more than 10 percent. If any of the sources have a finite angular diameter, 
a further correction will be needed ; this point is discussed below. 

The values of flux density at 500 Mc/s may be in error owing to the noise 
fluctuations of the trace and to the effects of confusion with nearby sources. ‘The 
combined effect of the filter and the subsequent numerical summation of the daily 
records should give the same noise level as an integration time equal to the total 
time of observation, which for most of the sources was about 7 hours. This 
estimate of the noise uncertainty was confirmed by the differences between the 
amplitude of the ‘‘ sin ” trace alone, and that of the “‘ cos ”’ trace alone (each trace 
leads to an independent estimate of the amplitude, with an error appropriate to 
one half the total observing time). 

The effect of confusion with nearby sources may be computed from the 
observations at 81-5 Mc/s (1), assuming that the nearby sources have the same 
spectra as those investigated here; this computation shows that the effects of 
confusion are less than the residual noise fluctuations and may therefore be ignored. 

3. The effect of a finite angular diameter.—The apparent amplitude of a source 
as measured by a 470A interferometer will be less than the amplitude measured at a 
small spacing if the angular diameter of the source exceeds a few minutes of arc. 
‘The amplitude of a radio-star relative to that of Cygnus A must therefore be 
corrected to allow for the angular extent of Cygnus A, which is some three 
minutes; a further correction will be needed if the radio-star itself has an 
appreciable diameter. 

Measurements of the angular extent of Cygnus A at 500 Mc/s (7) show that + 
(i.e., the factor by which the amplitude is reduced) for an east-west spacing of 470A 
iso-75 +005. Comparison of these measurements with those at other frequencies 
(8, 9) suggests that y is independent of frequency, at least within the range of 
frequency 125-500 Mc/s. It will be assumed that this result is true also of the 
other sources; measurements at different frequencies may therefore be used to 
estimate y at 500 Mc/s. 
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‘The angular diameters of the fainter sources at 159 Mc/s have been obtained by 
Mr S. Archer (personal communication), who directly measured the reduction in 
amplitude at a spacing of 308A, and, assuming the brightness distribution across 
each source to be Gaussian in shape, estimated the angular diameters. As the 


Tas.e II 
The measurements of angular diameter 
Source Diameter y (at 470A) 
50 <2’ 0°76 to 10 
379 3° t0'5 o55t0'1 
404 to 0°93 to 1°0 
440 3 0°76 to 1°0 
724* <12" 1‘o 
806t 1°6’ 0°85 
1175* <s2" raze) 
1259 <i1''5 0°87 to I'o 
1432 ast 0°65 +0°'2 
1775 <i1''s 0°87 to 1'0 
1843 <—s 0°93 to 1'0 


* Measured by Morris, Palmer and Thompson (10) 
+ Measured by Carter (11) 


extrapolation from 308A to 470A is practically independent of the shape of the 
source, however, the derived values of y at 470A do not depend upon the assumed 
shape. Archer’s results are listed in Table II, together with two measurements by 
Morris, Palmer and Thompson (10). y is appreciably different from unity in 
the case of sources 379 and 1432. ‘The uncorrected amplitudes of the remaining 
sources are not in error by more than 25 per cent, and in most cases by even less. 
These errors may be ignored. 


TasBLe III 
The results of the observations at 500 Mc/s 
Source Relative Relative Flux Density Spectral 
Amplitude Amplitude at 500 Mc/s Index 
81°5 Me/s 500 Mc/s 10° MKS 
50 10 5'ot2 18 —o's 
379 10 8:0+6 29 —o'8 
404 28 12°5+4 45 —0°85 
440 12 <2°7 <1o —1°2 
724 9 4'°9+2 18 —0'85 
806 40 I7*1+5 62 —1'o 
1175 9 7442 27 —0'6 
1259 10 4°'8+2 17 —1'2 
1432 63 18:0+6 65 —1°3 
1775 9 <319 <9 Pre | 
1843 14 <1'7 <F “19 
Cyg A 1000 1000 3600 - 


4. The results.—The results of the observations at 500 Mc/s are summarised in 
Table III, in column 3 of which the relative amplitudes of the sources are shown, 
with that of Cygnus A normalised to 1000 units. Corrections have been applied 
to allow for the non-linearity of the receiving system, and for the finite angular 
extent of Cygnus A. In addition, corrections have been applied to the amplitudes 
of sources 379 and 1432, to allow for the diameters of those sources. 
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For most of the sources, the error quoted in ‘Table III corresponds to the noise 
uncertainty. For five sources, however, Nos. 404, 724, 806, 1175 and 1432, this 
uncertainty is a small fraction, and the quoted error has been increased to include a 
possible error from the linearity calibration, described in Section 2, the accuracy of 
which is estimated at 25 per cent. 

The values of flux density in column 4 of Table III are dependent upon the 
value for Cygnus A, as corrected by Whitfield (12). So long as the absolute 
flux density of Cygnus A is still under discussion, therefore, the values are all 
provisional, but Whitfield’s arguments make it unlikely that the corrected value is 
in error by as much as the error of the present measurements, i.e., 25 per cent. 

‘The new results have been used by Whitfield to derive the mean spectral index 
for each source, defined as the index x in the relation S(f)ocf”, where S(f) is the 
flux density at frequency f. Whitfield’s values derived from all available observa- 
tions in the range 30-500 Mc/s are given in the final column of Table III. Froma 
study of a larger number of sources Whitfield (12) shows that the spectral indices 
for extra-galactic sources are lower than for those identified with galactic objects, 
and that unidentified sources such as the eleven investigated here form a class 
with a still lower spectral index. ‘The mean spectral index derived from the present 
measurements is — 1°10 + 0°10, with a mean error per point of 0-3. ‘l’o include the 
experimental error at 500 Mc/s, which is 25 per cent, the error in the mean should 
be raised to + 0°15 (in this calculation it is assumed that the errors introduced by 
the measurements at lower frequencies are unimportant). According to Whitfield 
the mean spectral indices for galactic and extra-galactic sources are — 0°74 and 
—1°05 respectively. ‘he mean spectral index of the sources investigated at 
500 Mc/s thus differs significantly from the value for sources identified with galactic 
objects. 
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ERRATA 


M.N., 116, No. 6, 1956: 
M. A. Ellison, A study of a blow-off prominence (1956 May 18). 
p. 626, last two lines, 


for 20 km/sec read to km/sec; 
for 600 km/sec read 300 km/sec. 


David S. Evans, The sense of rotation of NGC 253. 


Plates 9 and 10, facing pp. 660 and 661 respectively. 


‘The angular scale was not changed to conform to the reduction of the original photographs. 
The actual scale of the photographs reproduced is 15”:2/mm. 


M.N., 117, No. 3, 1957: 
Obituary Notice: John Evershed 
Pp. 254, penultimate paragraph, first and third lines respectively, 


for 1950 read 1953; 
for Nancy Acworth read Nancy Acworth Orr. 
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